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XXV. The Onset of Convection by Thermal I nstability in Spherical Shells 


By 8S. CHANDRASEKHAR, F.R.S. ° 
Yerkes Observatory, University of Chicago* 


[Received November 1, 1952 


SUMMARY 


In this paper the problem of the thermal instability of an incompressible 
sphere consisting of an inviscid core and a viscous mantle is considered, 
and it is shown that the pattern of convection which sets in, at marginal 
stability, in the mantle shifts to harmonics of the higher orders as the 
thickness of the mantle decreases. Thus, when the mantle extends to 
a depth of half the radius of the sphere, the harmonics of orders three and 
four set in about simultaneously, while the harmonic of order five follows 
very soon afterwards. The bearing of this result on the problem of 
convection in the earth’s mantle and of the interpretation of the earth’s 
topographic features is indicated. 


§1. [INTRODUCTION 


In a recent paper (Chandrasekhar 1952: this paper will be referred to 
hereafter as I) the problem of the thermal instability of a fluid sphere 
heated within has been treated. In this paper the related problem of 
the onset of convection in spherical shells will be treated. The interest 
in this latter problem arises, in part, in connection with Vening Meinesz’s 
(1951) interpretation of a feature of the earth’s topography, to which he 
drew attention at the .same time. The feature in question is this: 
Using Prey’s (1922) development in spherical harmonics of the earth’s 
topography, Vening Meinesz combined the coefficients of the (2/-+1) 
tesseral harmonics belonging to a given order / to obtain the root mean 
square amplitude with which a harmonic of order / is present in the 
development ; and he showed that, apart from the strong first-order term 
(corresponding to the near separation of the land and the sea hemispheres 
on the earth), the orders most predominantly present are three, four and 
five. From this predominance of the orders three, four and five Vening 
Meinesz concluded that phenomena in which these orders play an impor- 
tant part must have occurred at some time in the earth ; and he suggested 
that thermally induced convection might be the principal responsible 
phenomenon. On this suggestion the differentiation of the surface of 
the earth into a land and a sea hemisphere (i.e. into a hemisphere of high, 
respectively, low topographic elevation and corresponding differences 
in crustal composition) is assumed to have taken place at a time when 
the whole earth was a highly viscous fluid sphere of nearly uniform 
" *Communicated bythe Auth, 
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density, and the easiest pattern of convection to excite (and the last which 
might be expected to survive in a cooling earth) is a spherical harmonic 
of order one (cf. I). And the predominance of the orders three, four 
and five in the topography is traced to a subsequent stage in the develop- 
ment of the earth when it had reached a constitution not unlike the 
present with a core and a mantle of thickness of about 2900 km. The 
surmise here is that convection in the mantle at this later stage ‘ tore ’ 
the land hemisphere into something like its present topographic dis- 
tribution (cf. also Jeffreys 1952 in this connection). The validity of this 
surmise requires that in a mantle extending to approximately half the 
radius of the sphere the modes of convections most easily excited are 
those of orders three, four and five; we shall see that this is confirmed 
by the detailed calculations presented in this paper. 


§2. THE FORMULATION OF THE PROBLEM 


Consider an incompressible sphere of radius R in which a stationary 
temperature distribution of the form (I, eqn. (2)) 

T =B(R?—1°) oe P o eile ee ee 
can be maintained in virtue of a uniform distribution of heat sources’ 
imbedded in it. We shall suppose that though the sphere is of nearly 
uniform density it consists of an inviscid fluid core extending to a fraction 
y of the radius of the sphere and a viscous mantle overlying the core. 
We wish to determine, under these circumstances, the criterion for the 
onset of convection by thermal instability in the mantle. 

For the application to the earth’s mantle contemplated in §1 the 
model we are adopting is admittedly a very idealized one: for, in the 
earth, what principally distinguishes the mantle from the core is the 
existence of a relatively sharp discontinuity of density between them. 
In the model we have described we are not allowing for any such dis- 
continuity. On the other hand, by assuming that the core is inviscid and 
the mantle viscous, we are making the interface between them a free 
surface across which no tangential viscous stresses can act (cf. I, §3); 
and we may expect that this latter condition, appropriate for a free 
surface, will be actually fulfilled in the earth at the interface between the 
core and the mantle. Consequently it would not appear unreasonable 
to hope that the consideration of the model we have described would 
disclose the essential features of the convection patterns which may 
predominate in the earth’s mantle. 

Under the conditions we have stated the equations governing the 
problem are the same as in I. In particular, the criterion for the onset 
of convection in the mantle as a disturbance in the form of a spherical 
harmonic of order 1 will be determined by the conditions that the 
equation (cf. I, eqns. (41) to (43)) 


GP 2d Uri p ee 
ita ge a W=—lI+1)C,W eee 4) 


y2 
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admits a solution satisfying the boundary conditions 


_@W T@® 2d +1) 


ade adr “rar 


2 
| W=0at7r=land r=, . (3) 


where it may be recalled that 7 is now measured in units of R and 
(I, eqn. (42)) 
2 
C= 2BY pe. . . . . . . . . . (4) 


KV 


(For the meanings of the other symbols see I.) 


§3. THE VARIATIONAL SOLUTION OF THE PROBLEM 


It can be readily verified that the variational principle established in 
I, §5is applicable, with suitable modifications, to the problem on hand, and 
that it can be used to solve the characteristic value problem formulated 
in §2 in the following manner. 

Letting 


(5) 


Orn r 


@ 2d La ee, 


assume a form for Ff (containing one or more parameters A;) which 
vanishes at r=1 and r=7; then solve eqn. (5) for W and arrange that 
it satisfies the remaining boundary conditions 


_@w 
te 


Ww a OME Abt al ANC t= 9 os, ees. =" (6) 


and finally, evaluate C, according to the formula (cf. I, eqn. (54)) 


Lf, (PR fdrP UL FS dr 


C= (7) 
: L(l+-1) eC r2W Fdr 


and minimize it with respect to the parameters A,. In this way we shall 
obtain the best value of C, for the chosen form of F. 

The form of the equation to be solved for W suggests that we assume 
for F a solution of Bessel’s equation of order (J+) which vanishes at 
r=landr=y. Thus, writing 


FHS 14.4)2(47)+OI _ci4sjay(&); gen oes le pee) 


where a, 6 and « are constants, we shall arrange that F vanishes at 
r—landr=y. For this purpose we proceed as follows: Let 


@ 11/2,» (2)=I_cr4sjay (en) (2) —F 14 142(0%0)J—(2), one (2) 
where « is a constant unspecified for the present. Then, 


€ pre, rya(2r =F crayon %M)F resya(ar) —F resyo(om)F—(a41/2)(2") sir. (10) 
RZ 
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clearly vanishes at 7=7; it will also vanish at r=1 provided 
Trap 0q)d r4172(%) —S py 1j2o(%)I—cr4sya(%)=9.  . « . (11) 
It is known (cf. Gray and Mathews 1922, p. 82, theorem X) that 
eqn. (11) admits an infinite number of roots all of which are real and 
simple; and that if «; (j=1, 2, ...) are the distinct roots of the equation, 
the functions @), 5/2, p41/9(%j”") (=1, 2, ...) form an orthogonal set with 
the integral property 


1 
{h PE riryo, rerpal%q) ryryo, r41p2(%x") 7=N 5, p43729jn, + +» (12) 
where 6,, denotes the usual Kronecker symbol and 
2 (J*144/0(a;7) 
M,, see gage pegs tt. 8 
j, 1+1/2 ajar” J” 144/2(%;) ( ) 


For later use we may note here that the derivatives of @), 1/2, 41 jo(%j) 
at r=1 and r=y (which we shall denote by @’;,4)9(a;) and @’;,4/9(a,;n) 
respectively) are given by 


; d 2/ a, 
@nnld=[ gp @rene, onto”) | (1 eda 


and 


d 
© 141j9(%j) = Or @ psrjo, renal) | =(-1/= «oe satis (15) 
bh 
(Equations (12) to (15) can be derived by following the procedure given 
by Carslaw and Jaeger 1947, § 79, in obtaining similar relations involving 
the Bessel functions of order zero.) 
Also, since @;,1,,(z) represents a solution of Bessel’s equation of 
order y it must clearly satisfy the relations implied by the equations 


d v 
Fra 1 14/9, rp 1(2)}= +2" 1 170, (z) 
and Fie 


d 
ae {2-16 1119, »-1(2)}= Po @ cstee (z). 


Returning to the variational solution of the problem described in the 
beginning of this section, we shall assume for F the trial function 


il 
EO Fe 1s Siaays, 141/2(%)7), + iso ha) Rane 


where the A,’s are the variational parameters. With this choice, F 
vanishes at r=1 and r=y as required. Moreover, with F given by 
eqn. (17), eqn. (5) governing W can be explicitly solved and we have (Ciel 
eqns. (68) to (73)) of 


1 A, 


Waes Vr aad Cis, 412%") + By + Bor? + Bo U4 D4 By C-1), . (18) 
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where B,, B,, B, and B, are constants of integration to be determined 


by the Bomidary Borduons (6). These boundary conditions lead to 
the equations 


Bre b. oh 0. Ae ee ee 119) 
By + Bor ?+ By ae WE a0 eee tis .(20) 
U(d—1)(B, +B,)+ + 1)(0+2)(By+ Bs) =25;( A/a?) p41)2(a4) . (21) 
and. 
U1) By? Boy) + UENCE 2 Bolt Boy) 
= 2 972) (A 5/0052) @ 141)0(2-9n), - + + (22) 
where @’;,1/2(%;) and @’,,4/2(%) are defined as in eqns. (14) and (15). 


Combining eqns. (19) and (21), and similarly (20) and (22), we obtain 
ibrea rl 
B, + B,=—(B,+B;)=— Wr one GS isja(%j) » - (23) 


and 
By + By Y=—(B age ys ee!) 
yah) 
A ae a; 0; ss €' 141/2(%7) 5 
and solving these equations for B, and B, we have 
1 A ees : 
B,= (QF 118) 24% ap il @ 141/9(0)—1 3? reuelon | ae) 


-and 


; Be 
By= (21-+1)(y GG ~ hae 1+1/2(% le, rerelopn) | (26) 


Turning next to the eae of C,, we find that for / and W given 
by eqns. (17) and (18) the integrals occurring in (7) can be reduced in 
a manner quite analogous to that followed in the reduction of the 
corresponding integrals in I (cf. 1, eqns. (76) to (78)). Thus, on using the 


orthogonality relation (12), we now find that 


{{e(q) eye: ‘hdr 2A; BN as eet ye) 
n 


and 
Pe, D(A, aj )N;, r1/2 


EA, | @ ryayo, reajal% Br? B, pt Bote Bart +52) dr, 
(28) 

The remaining integrals on the right-hand side of (28) can be evaluated 
if proper use is made of eqns. (16), (23) and (24); and we find, after some 
lengthy reductions, that 
[iP Wrdra=ZiAplAN5, rp 

—2By, Zi (A, /a;? MN v4 172, 145)2(%)— resi ME UE 1+5/2(%77)] 

2B, A,|«,? NE tss/2, 1-aj2(% an =e he Ma 1-3/2(%5)) |- (29) 
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The resulting expression for C, is, therefore, 


21-41 


C= TT) 


(3A aN, yal | LENE (A, 259N5, co 


My A; 
— rr 2 3 pa syal(ay) 
(1—7?!+) j a8 14+1/2\%9 


A, / 
—7f 2" 1 1 j0(xjn) 1&5 ga Gis, 145/2(%1) 1°? E 14172, 145/2(%27)] 


2 An 
~ GFA) re: [@' p4.1j2(%,) 


A; 
ip) HIRE io(a aj) = 53 Gate, ie 3/2(% x) 


al 
—n HP io, -an(e)] | > (30) 


where we have substituted for B, and B, according to eqns. (25) and (26). 
By minimizing this last expression with respect to the A,;’s we shall obtain 


the ‘ 


Table 1. The Characteristic Numbers C,; for Various Values of / and 
l n=0:2 7 =0-4 7=0°5 7=0-6 7=0°8 
1 4-141 108 | 2-140 104 | 5-373104 | 1-73410° | 8-604 x 10° 
2 3°826 10% | 1-041 x 104: | 2-347 10* | 6-962x10*) 3-017X%108 
3 5:790 x 10 | 9-730 108 | 1:803x104 | 4602x104 | 1-626x 10° 
4 9-441 103 | 1200x104 | 1-889x104 | 3-938x10* | 1:075x 105 
5 1-486 x 10* | 1-649x104 | 2-175x104 | 3-939x104 | 8-070x105 
6 2-331 x 104 | 2-767x104 | 4:340x104 | 6-610x 105 
7 3°281 x10* | 3-633 104 | 5-072104 | 5:779x 105 
8 4-808 x 104 | 6-130 104 | 5-312 105 
9 6-340 x 104 | 7-545x104 | 5-079x 105 
10 8-276 x 104 | 9-358 104 | 5-012 105 
11 1-067 105 | 1-162x105 | 5-073 10° 
12 1-357 xX 108 | -1-439105 | 5-239x 105 
13 1-705 105 | 1:774x105 | 5-499 x 105 
14 2-115 10° | 2-172x105 | 5:848 x 105 
15 2-594 10° | 2-640x105 | 6-282 x 105 


best ’ value for C’, for the chosen form of Ff’. 


The simplest trial function for / of the chosen form is 


1 


f= 0 @ i+1)2, 1+1p2(%17")s (31) 


where «, is the first zero of the expression on the right-hand side of 
eqn. (11); @p44/9, r41/2(%,"7) has then no node between r=y and r=1. 
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For this choice of F there is no variational parameter with respect to 
which we have to minimize, and eqn. (30) directly gives 
C= ene 1+1/2 {2+ Vo Ny ity 

ae am [O14 1/o(o)—9f 92 @ 1 alan) IS p70, 145/2(%1) 

Cy so, 145/2(%17))]— ry [@" p41yo(oe1) 


7 PC spl @ i412, 1-3/2(%1)—9 9° 1a r-aelan)]} 


al 


(32) 
Fig. 1 
a SILC oe ee haa A 
Marites 3 SS, ah) Ca ee ee oe 
The criterion (C;) for the onset of convection in the mantle as a spherical 


ic di i i the mantle. 
harmonic disturbance of order / for various thicknesses of 
The different curves are distinguished by the values of the fraction (%) 
of the radius of the sphere occupied by the inviscid core. 


In table 1 are listed the values of C, obtained with the aid of formula (32) 
and the recent tabulation of Chandrasekhar and Donna Elbert (1953) 
of the roots of eqn. (11). The dependence of C; on / and y is further 


illustrated in fig, 1. 


— 
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Fig. 2 
l 
I 
10 
9g 
8 
7 
6 
Bae 
rae 
o 
; 
I 
02 0-4 0-6 0-3 ” 1-0 


The order of the spherical harmonic / of the convection pattern which appears 
in the mantle, at marginal stability, in its dependence on the fraction (7) 
of the radius of the sphere occupied by the inviscid core. 


Fig. 3 


Ll 2 3.4 5 eR ee oy et end a een meee 


Comparison of the root mean square amplitude (in km) of a spherical harmonic 
of order / in the topographic analysis of the earth’s surface (the full line 
discontinuous curve) after Vening Meinesz with 1/C; (the discontinuous 
curve passing through the open circles) in the case when the mantle extends 


to a depth of one-half of the radius of the sphere. (The curves have been 
adjusted to agree at 1—5.) 
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The values given by (32) could be improved by including a second term 
in F. However, in view of the fact that for the case n=0 the first 
approximation already gave values accurate to a fraction of a per cent 


(cf. I, table 1), it was not considered necessary to carry out this improve- 
ment. 


§4. Discussion of THE RESULTS aND APPLICATION To CONVECTION 
IN THE EartH’s MANTLE 


From fig. 1 it is at once apparent that, as the thickness of the mantle 
decreases, the pattern of convection which first appears at marginal 
stability shifts progressively to harmonics of the higher orders. This 
is shown more clearly in fig. 2, in which the order of the harmonic which is 
the first to be excited by thermal instability is plotted as a function of 
the fraction of the radius occupied by the inviscid core. From this 
figure, as well as from fig. 1, it is seen that for 7 ~ 0-5 the harmonics of 
orders three and four set in about simultaneously, while the harmonic 
of order five follows very soon afterwards. This is a confirmation (go 
far as it goes) of the suggestion of Vening Meinesz that the topographic 
distribution in the land hemisphere may have arisen as a consequence of 
convection in the earth’s mantle. A more quantitative test of this 
hypothesis is possible if, following Vening Meinesz, we compare the 
relative amplitudes with which the spherical harmonics of the various 
orders are present in the topographic analysis of the earth’s surface, with 
1/C, as a measure of the probability with which a mode of disturbance 
of order / will be present in fully developed convection. This comparison 
is made in fig. 3. It will be seen that the agreement between the two 
curves is very satisfactory; indeed, as satisfactory as that found by 
Vening Meinesz (1951) on the basis of an empirically founded theory of 
convection in spherical shells. 


In concluding this paper I wish to record my indebtedness to Miss 
Donna Elbert but for whose assistance with the numerical work this 
paper might not have been written. My thanks are also due to Professor 
G. P. Kuiper for drawing my attention to the geophysical interest 
attached to the problem considered in this paper. 
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ABSTRACT 


Eight hundred tracks of slow particles have been photographed during a 
cloud chamber experiment on penetrating showers. Measurements of tracks 
produced by ten heavy, negative, non-decaying particles are presented. 
The total observed time of flight of these particles is 1-6 10~* sec. 
After combining this time with the number of slow V-particles observed 
to decay in the cloud chamber, it is concluded that at least one type of 
heavy meson must have a mean life greater than 4x 10~° sec. 


§ 1. INTRODUCTION 


THe photographs discussed here were obtained in an experiment on 
penetrating showers at the Jungfraujoch (3580 metres above sea level). 
A description of the experimental arrangement has been given already 
(Astbury et al. 1952). The cloud chamber had an illuminated volume 
of 54549 cm? and was operated in a magnetic field of 5800 gauss. 
All the photographs considered below were selected from 10000 taken 
in the first 26 weeks of operation. They were obtained with selection 
systems which are illustrated in figs. 1 (2) and (6) of our previous paper. 

The photographs showing the decay of V+-particles were discussed 
in our previous paper and it was concluded there that the mean life of 
the particles was at least of the order of 10-® sec. More valuable 
information can be obtained from an estimate of the number of charged 
V-particles that have been observed to traverse the cloud chamber 
without decaying. The present paper shows how such an estimate 
has been made for the slow, heavily-ionizing particles. 

Mesons with masses of about 1200m, can be distinguished by their 
ionization from 7-mesons if their momenta are betWeen 3107 ev/e 
and 610% ev/c. However, visual estimates of the ionization density 
in a cloud chamber track are so unreliable that it is difficult to distinguish 
tracks due to particles of mass 1200m, from those due to protons. For 
this reason in the following analysis an attempt has been made to identify 
only negatively charged heavy mesons. 


§ 2. MEASUREMENTS AND RESULTS 


Tracks of 800 slow particles were selected for analysis; the criteria 
for selection were that the particles should be associated with showers 


se er a is ee 


* Communicated by Professor P. M. S. Blackett, F.R.S. 
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and that their tracks in the cloud chamber should be at least 10 em 
long. The curvatures of 600 tracks were measured by comparison with 
standard curves. The remaining 200 tracks, including those with 
smallest curvature, were measured with a travelling microscope. The 
errors in these measurements were always less than the probable spurious 
curvatures introduced by convection in the cloud chamber or by multiple 
scattering of the slow particle. 

The specific ionizations of the 800 tracks were estimated by a single 
observer. 200 tracks were then excluded from the analysis because 
their estimated ionizations were greater than 10. A further 180 tracks 
were excluded because either the technical quality of the photograph 
or the position of the track in the cloud chamber made a reliable estimate 


Fig. 1 


TRACK CURVATURE (m7?) 


The curvatures and estimated ionizations of 400 slow particles. The broken 
lines give the expected values for particles of mass 275m,, 1000m, and 
2000m,. 


of ionization impossible. Finally, tracks of twenty apparently negative 
particles were rejected because their direction of flight was not certain. 
Event CC 17 (Plate 11 (b)) showsan example of a particlein this last category. 

In fig. 1 the estimated ionization densities of the remaining 400 tracks’ 
are plotted against curvature. ‘To avoid confusion no probable errors 
have been indicated on the figure. The actual estimates of ionization 
were made by setting upper and lower limits which differed, in general, 
by a factor of nearly two. The mean value of these limits for each track 
has been used in plotting the figure. 
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Figure 1 also shows the theoretically predicted values of ionization and 
curvature for singly charged particles of mass 275m,, 1000m, and 2000m,. 

The measurements made on ten tracks apparently due to negatively 
charged heavy particles are presented in table 1. These tracks were 
all measured with a travelling microscope and their ionization densities 
were estimated independently by three observers. The errors in the 
momenta given in column 2 are probable errors including errors of 
measurement, distortion due to gas motion in the cloud chamber, and 
spurious curvature due to multiple coulomb scattering of the particles. 
The ionization values given are the extreme limits of the estimates 


Table 1. Negative Heavy Particles 


Ionization Mass Track length | Time of flight 
(L/I) CL) (cm) (10~® sec) 


Event | P(108ev/c) 


P31 
JM41 
JV76 
DN20 
KD79 
JC3 
JY82 
DP22 
CH19 
FK24 


>1100 20 
>1030 2: 
1230 to 3600 
5 to 10 1580 to 4600 
630 to 1400 
700 to 1900 
950 to 3400 
> 660 
910 to 2800 
> 1000 


KS DK ON RRR 
ASDwWwOwowrH aa 
HEHEHE HE HE HE HE He HEH 
CoOooe oof Oo SC So 
Gc HB Si de WD Go GO DD 
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made by the three observers. The limits given in column 4 to the mass 
of each particle have been calculated by combining the lower and upper 
limits of momentum and ionization from the previous columns. Column 5 
gives the length of each track visible in the cloud chamber, and column 6 
the observed time of flight of each particle. The times of flight (7’) 
are measured in the rest systems of the particles, and have been 
calculated from the relation L=(P/M)T, where P is the measured 
momentum and M the mass of the particle. In calculating the values 
of 7 the masses of all the particles in table 1 have been taken to be 
1200m,. The photographs of three events listed in the table are 
reproduced in Plate 10 (a and 6) and Plate 11 (a). 


§ 3. Discussion 


Bearing in mind the size of the errors discussed above, the main 
features of the results shown in fig. 1 can be readily explained. Most 
of the tracks measured were produced by heavy positive particles. 
The majority of these are certainly due to protons, but deuteron and 
a-particle tracks probably broaden the distribution of points on the 
figure. The number of heavy positive mesons is discussed in § 4. 


——— 
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The tracks of positive and negative light particles occur in roughly 
equal numbers and can all be interpreted as due to 7- or p-mesons. 
Since they are all associated with nuclear interactions these light mesons 
are probably 7-mesons, but this could not be deduced with certainty 
from our measurements. 

The points in fig. 1 apparently due to heavy negative particles are the 
most interesting feature of the analysis. We show below that few 
of these points can result from errors of measurements or interpretation ; 
the majority represent non-decaying heavy mesons. If this is so, the 
measurements in table 1 can be used to give a lower limit to the mean 
life of the particles. 

There are two possible causes of erroneous classification of tracks 
in this group. First, gross distortion in the cloud chamber could have 
resulted in a negative 7-meson appearing to have a very small curvature 
or in a proton track having its curvature reversed and appearing to be 
negative. Second, and more important, a positive particle travelling 
towards an interaction could have been mistaken for a negative particle 
coming from the interaction. No imaginable error in the ionization 
estimates could cause any light meson to have been included among the 
particles listed in table 1. 

The figures in table 1 allow the likelihood of gross distortion to be 
estimated. The probable errors in the momenta are statistical in the 
sense that they are derived from the known performance of the cloud 
chamber. In every case the measured curvature of the tracks is more 
than three times the probable error in curvature. Further, a distortion 
of more than five times this probable error would be required before a 
proton track of the same ionization would have the measured curvature. 
The distortion necessary to make a 7-meson appear heavy is still greater 
and this source of error need not be considered. Since the total number 
of tracks included in our analysis is 400, it is unlikely that even one 
of the group in table 1 has been misinterpreted as a result of distortions. 

The probability of a proton moving towards a nuclear interaction 
and being interpreted as a negative particle produced in the interaction 
can only be calculated approximately. This is because no satisfactory 
criterion of ‘ association’ can be laid down which is applicable to all 
events. A nuclear interaction near the cloud chamber has to be located 
by projecting backwards the visible tracks in the chamber. The precision 
with which the site of the interaction can be determined depends upon 
its distance from the chamber, the number of secondary fast particles 
that cross the sensitive volume and the angular spread of these particles. 
In addition, when a slow particle is emitted from the interaction, its 
probable scattering in the material above the cloud chamber may be 
several degrees and its association with the interaction may be obscured. 

A loose criterion of association can be used to give an upper limit 
to the number of spurious association among the ten events. For this 
purpose any track is assumed to come from a shower if it meets the 
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chamber wall (or the lead inside the chamber) at a point less than one 
centimetre from the shower axis. This condition is generally less 
stringent than those applied to the tracks in table 1 since it ignores the 
direction of the track. However, event KD79 (Plate 11 (a)) does not 
satisfy the condition—magnetic deflection alone would account for the 
spread of the particles in this event. 

Using this criterion and making the additional assumption that all 
heavily-ionizing particles are randomly orientated in space, the probability 
of a chance association between any single track and a given interaction 
is found to be ~0-:002. <A spurious association cannot arise, however, 
unless more than one nuclear interaction occurs near the chamber in 
the same event. If n is the average number of interactions producing 
tracks in the chamber in each event, the geometrical probability (0-002) 
must be multiplied by (n—1) to find the probability of a spurious 
association of any proton track with any nuclear interaction. The value 
of n is between 1-5 and 2 so that the final probability, with these assump- 
tions, for a proton being interpreted as a heavy negative particle is ~0-002. 

Among the original 800 tracks analysed there were 600 where this 
misinterpretation might have occurred. Multiplying this number by 
the calculated probability, we find the expected number of spurious 
events in table 1 to be rather greater than one. 

All the particles listed in table 1 are travelling downwards in the 
cloud chamber. A rather more refined geometrical argument based 
on the number of upward moving protons (84 in all) and the precision 
with which each of the ten nuclear interactions can be located leads to 
a result similar to that given above. 

From these considerations it is likely that at least seven of the 
particles in table 1 are negatively charged and heavy (mass> 700m,). 
In § 4 only the total time of flight of these particles is used ; this is taken 
to be 11x10-* sec to allow for errors in interpretation. It is worth 
emphasizing that the particles with the longest times of flight are those 
whose tracks are longest and most curved. Therefore they are the most 
accurately measured and the least liable to misinterpretation. Thus 
the total time of flight has probably been under-estimated, 


§4. THe Mean Lire or CHarcep V-Parricies 


From the discussion in §3 we conclude that negative, non-decaying, 
heavy mesons have been observed for a total time of ~11 x 10~° sec. 
In our earlier paper (Astbury et al. 1952) we reported measurements 
on two heavily-ionizing charged V-particles which decayed inside the 
cloud chamber. Table 2 reproduces the relevant measurements. 

If we assume that all the non-decaying particles discussed above are 
of the same type as the decaying V-particles an estimate of their mean 
life can be made. 

No positive heavy mesons have been identified directly because of 
the uncertainty attached to ionization estimates, It is known, 
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however, that decaying positive and negative V-particles occur in nearly 
equal numbers (Armenteros et al. 1952, Astbury et al. 1952). The best 
estimate of the time for which positive non-decaying particles have been 
observed is therefore assumed to be 11x 10~® see also. 


Table 2. Slow Charged V-Particle Decays 


Event Track length (cm) PIM Time of flight (10-° sec) 
JB127 13 0-47 0-9 
JK23 10 0-40 0:8 


These figures refer to the primary particle before decay. 
These events are events 2 and 3 in table 4 given in Astbury et al. (1952). 


Combining this figure with those in table 2 we find that slow particles 
of both signs have been observed for a total time of 24 10~° sec. 
During this time two particles have been observed to decay. This leads 
directly to an estimate of the mean life 7~12x10~° sec. The 
statistical error in this estimate is very large but, if the mean life were 
as short as 4X 10-® sec, the probability of our observing as few as two 
decays would be less than 10%. 

It is possible that the decaying and non-decaying particles are not 
all of the same type. Three distinct types of heavy charged meson 
have been identified in the photographic emulsion. These are the 
k-, x-, and 7-mesons, all of which have masses between 950m, and 
1500m,. These mesons decay according to the following schemes 
(Powell 1952): 
K+ > w*+2 neutrals, 
y+ > 7+-+1 neutral, 
ce >qgt+tattor. 

The ¢-meson is excluded from the discussion because its mass (500 
to 570m,) is outside the mass limits for any of the particles in table 1. 

The non-decaying heavy mesons could, from our measurements, 
be identified with any of the three types listed above. The decaying 
V-particles cannot, from their mode of decay, be 7-mesons but could 
be «-megons, y-Mesons or a mixture of the two. If there were no 
information about the relative abundance of the three different types 
of meson, we could only conclude that at least one type must have a 
mean life of the order of 10-8 sec to account for our non-decaying 
particles. However, if a large fraction of these particles were 7-mesons, 
it would imply that t-mesons had a long mean life (since none are 
observed to decay) and that they were abundant. if this were sO, they 
would be expected to decay at rest in the photographic emulsion relatively 


frequently. What evidence there is (Menon 1951) suggests that this is 
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not so; 7-meson decays appear to be much less frequent than those 
of other heavy mesons. It follows that most of our non-decaying particles 
are probably «- or y-mesons, and that at least one of these mesons must 
have a mean life greater than 4 10~° sec. 


$5. Tae Mean Lire oF 7-MEsSons 


For interest, the same method of analysis has been used to estimate 
the mean life of thez-meson. Thirty-three tracks due to slow, light mesons 
of both signs were identified among the original 800 tracks. The sum 
of the observed times of flight of these particles is 5-3 10-8 sec and 
three examples of 7->. decay can be identified. The assumption that 
all the light mesons are 7-mesons leads to a most probable mean life 
of 1:8x10-8 sec. This result is in good agreement with the accepted 
value of 2:5 10-8 sec. 


§ 6. CONCLUSIONS 


An analysis of 800 tracks due to slow particles has shown the presence 
of a small group due to heavy negative mesons which cross the cloud 
chamber without decaying. The total time of flight of these mesons 
is about 10-® sec. Assuming that similar positive particles are 
equally numerous, it can be concluded that at least one type of heavy 
meson has a mean life of the order of 10~-* sec. It is most probable 
that the mesons observed are the x- or y-mesons which have been 
identified in the photographic emulsion. Whether the charged V-particles 
are «- or y-mesons or a mixture of the two remains uncertain. 
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ABSTRACT 


Statistical methods of estimation of the mean life-times of unstable 
particles from cloud chamber photographs are discussed (i) for decaying 
particles only, (ii) when non-decaying particles can be recognized. 

Some discussion is also given of the effect of momentum errors, and 
of estimation and discrimination problems for mixtures of particles. 


§1. InrropucToRY REMARKS 


In view of recent interest (see, for example, references at end of paper) 
in estimating the mean life-times of unstable particles from cloud chamber 
photographs, some notes on purely statistical aspects of this estimation 
problem are recorded below. They were initially stimulated by conversa- 
tions with research physicists at Manchester and it is hoped will be of 
use to other workers confronted with this type of statistical problem. 

It is well-known that there are many difficulties (see Wilson and Butler 
1952) in isolating unambiguous and accurate observations from the 
photographic material, quite apart from the usually sparse total amount 
of statistical information on the true mean life-times contained in such 
observations even if extracted free from error. The main purpose of 
these notes is, however, to discuss the estimation problem after the 
observations have been extracted, although a few remarks have been 
included (§ 3) on the effect of errors in the assessment of the momenta of 
the particles. How far some of the more ‘ precise ’ methods of estimation 
(more precise in the sense of relying on fewer approximations) will be 
justified on particular data will of course depend partly on the latter’s 
overall reliability, but it will be shown that the more precise methods 
are not much more complicated (and sometimes simpler) than the more 
usual ones available in the case of larger samples of data, and are much 
less liable to mislead. Some brief comments on the possibility of 
discrimination and estimation in the case of mixtures of more than one 
type of particle are also added. 


§ 2. OnE TYPE OF PARTICLE, WITH NO RECORD OF UNDECAYED PARTICLES 


We consider first the case of m observations for one type of particle, 
the time to decay from the (effective) edge of the chamber being, after 
conversion from the corresponding measured length, denoted by t, for 
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the rth observation. The associated total time available for decay is 
denoted by 7',. The variation of 7’, from one observation to another 
slightly complicates the estimation problem ; otherwise this is similar 
to one discussed many years ago in connection with problems of artificial 
radioactivity (see, for example, Peierls 1935, Bartlett 1936). For 
completeness, however, the maximum likelihood equation of estimation 
is summarized, especially as its use in this context calls for special 
comment and extension. 
The life-time probability distribution for all the particles is 
n n e—tr/é 
Ti 


1—e-fr/é ? (0<t,<T,), : . 0 (1) 
r=1 r=1 


given that t,<T'., where @ is the unknown mean life-time, or equivalently, 
A=1/0 is the average rate of decay. It should be noticed that t, may 
be measured from any previously chosen zero, provided 7’, is also (thus 
permitting a choice of origin which avoids, as far as possible, edge effects 
in the chamber). From (1), we obtain the estimation equation 
oO. gee 
rar ol he Neel Fe U,Q,/P,)=0, Shs a (2) 


where L,= log f(t,), u,=t,/0, U,=T,/6, and P,=1—Q, is the chance 
of decay 1—e~?r in the time interval (0, 7',). Equation (2) may be 
solved by iteration or interpolation ; the latter method is recommended 
in conjunction with the further extensions below, and its use will be 
illustrated presently. 


The ‘ information function ’ 


nm OU \ ue UrQ,\ 1 $ 
b2(—Fe)=p 2 (Se) gE. @ 


r=1 =1 aye 


say, where H denotes statistical averaging, in this case over the 
probability distribution (1), is best known as providing the reciprocal of 
the asymptotic error variance in large samples of the estimate 6 from 
eqn. (2), but it is also important in providing the exact variance of 
XOL,/00 in (2), whose true mean value YH(AL,/é0) is moreover exactly 
zero. 

There are two limitations in the use of (3) to provide a standard error: 
for @: firstly, it is an approximation valid only for small standard errors 
and secondly, it depends on the unknown quantity 0. The dependence 
of (3) is mainly on 1/6* for large U,; log @ will have an asymptotic error 
variance obtained from the same expression (3) with the 1/6? removed 
and thus under such conditions the second limitation can be largely 
removed by working on the log@ scale. However, in most cases the- 
U, are not large, and even this change of scale may be inadequate 

An alternative method of indicating the statistical accuracy of fe to 
give a ‘ confidence interval’ for 6, which assigns statistical limits to 6° 
such that the chance of @ being outside these limits is equal to some. 
chosen value 0-05, 0-10, ete. It is hardly practicable in the present. 
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context to obtain these limits exactly; it is, however, quite possible to 
obtain much more valid limits than indicated by the crude standard 
error of 6. We may consider directly the standardized quantity 


== n ay n 
S(¢)= 2 -- —_ Pa, ek as 
(= 3 (u—14+ fe) || By, (4) 


which has zero mean and unit variance. The only approximation now 
involved is to treat S(#) as normal (Gaussian), so that from a graph of 
S(@) against 6 the confidence limits corresponding to any probability level 
may be read off by equating S(@) to the associated value ». Thus for 
P=0-05 (sum of both distribution ‘ tails’) »=-+1-96, for P=0-05 (one 
tail only, or 0-10 for both tails), .=1-645, ete. 

An interesting feature of S(@) is that its value as A= 1/0-> 0 is 


A—43P,) 1 f(ZE\  o (1 
carat tan (az) +0(g)- ee) 


Workers in this field have realized that unless Xt.<}X7, no finite 
estimate 6 is possible, but (5) provides a rapid method of making more 
‘precise’ statements about @ in other conditions for which 6= 00 is a 
possible value. From it we see that in the neighbourhood of 1/9=0 the 
use of S(@) (which is invariant to any change of scale for 6) is approximately 
equivalent to the estimate 
1 2A(t,—4tT,) 1 : 
(9 DPD wpm 9 f°) 
The use of the 1/0 scale is also most convenient with the more exact 
eqn. (4), as in addition to 1/0=0, values of S(@) may be computed for 
convenient values of 1/9 such as 1/6=0-1, 0-2, etc. To facilitate this 
computation a table of the quantities 
1—UQ/P, I,=1-—U°Q/P 
has been prepared (table 1). 

Finally, it is noted that the remaining approximation, that S(@) is 
normal, may be improved if desired (for further theoretical details, see 
Bartlett 1953). Thus the sampling distribution of S(@) always has a 
positive skewness, and if this is corrected, the resulting quantity will be 
more nearly normal. A general formula for the third moment H{(eL,/00)*} 


of 0L,/00 is ; a 
alg. cl Op 
on (Set) +358 (— Spe): ani et metas (7) 
whence in this case we find 


Cates. 


In place of S(0) in (4), we may as a refinement when calculating confidence 
limits for @, use the amended quantity 


Sei Seah aca eee eet (0) 
s2 


with standard error 


2 
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Table 1 


= 


SSS SSS OS SS 


0-7275 
0-7534 
0-7769 
0:7981 
0-8173 
0-8347 
08504 


Fe ee 


0:8647 
0-8775 
0-8892 
0-8997 
0-9093 
0-9179 
0-9257 
0-9328 
0-9392 
0-9450 


09502 
0-9550 
0-9592 
0:9631 
09666 
0-9698 
0:9727 
0:9753 
0-9776 
0:9798 


‘0 
“il 
2 
3 
“4 
5 
6 
oi 
8 
“9 
‘0 
all 
“2 
3 
“4 
5 
6 
cy 
8 
“9 
‘0 
ol 
2 
3 
“4 
‘D 
‘6 
otf 
8 
AY) 
‘0 
| 
“2 
3 
“4 
5 
6 
ooh 
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Table 1 (cont.) 


U P(=I,) |1—UQ/P I, [xs}1 [xs]s 
4-0 | 0-9817 0:9254 06959 0-738 1-744 
4-] 0-9834 0-9309 OV7119 0-779 1-763 
4-2 0-9850 0-9361 0-7274 0-820 1-781 
4:3 0:9864 0-9409 0-7422 0-861 1-798 
4-4 0:9877 0-9453 0-7564 0-901 1-813 
4:5 0-9889 0-9494 0-:7700 0-942 1-828 
4-6 0:9899 0-9533 0-7830 0-981 1-841 
4-7 0-9909 0-9569 0-7954 1-021 1-854 
4:8 0:9918 0-9602 0-8072 1-059 1-865 
4°9 0:9926 0-9632 0-8185 1-097 1-876 
5:0 0:9933 0-9661 0:8293 1-135 1-885 
5-1 0-9939 0:9687 0-8395 TeLfl 1-895 
5-2 0:9945 0-9712 08492 1-207 1-903 
5-3 0-9950 0:9734 0-8584 1-242 1:91l 1s 
5-4 0-9955 0:9755 0-8671 1-276 1-918 
5:5 0-9959 0-9774 0:8754 1-309 1-924 
5:6 0-:9963 0-9792 0-8832 1-341 1-930 
lath 0-9967 0-9809 0-8906 1-372 1-936 
5-8 | 0:99:70 0-9824 0-8975 1-402 1-941 
5:9 0-9973 0-9838 0:9041 1-431 1-946 
6-0 0-9975 0-9851 0-9103 1-459 1-950 
6-1 0:9978 0:9863 0-9162 1-486 1-954 
6-2 0-9980 60-9874 0-9217 1-512 1-958 
6-3 0-9982 09884. 0-9268 1-537 1-962 
6-4 0-9983 0-9893 0-:9317 1-562 1-965 
6-5 0-9985 0-9902 | 0-:9363 1-585 1-968 
6-6 09986 0-9910 0-94.06 1-607 19710 
6-7 0-9988 0-9917 0-9446 1-628 1-973 
6°8 0-9989 0-9924 0-9484 1-648 1-975 
6-9 0-9990 09930 0-9519 1-668 1074 
7-0 0-9991 0-9936 0-9552 1-686 1979 

mil 0-9992 0-9941 0:9583 1-704 1-981 
aio 0-9993 0-9946 0:9612 Laz) 1-982 
7:3 0-9993 0:9951 0:9640 1-736 1-984 
7:4 0-9994 0-9955 09665 1-752 1-985 
7-5 0-9994 0:9958 0-9689 1-766 1-986 
7°6 0-9995 0-9962 0-9711 1-780 1-988 
fiers 0:9995 0:9965 0-:9731 1-793 1-989 
428 0-9996 0-9968 0-9751 1-805 1-990 
7°9 0-9996 O-9971 0:9768 1-817 1-990 
8-0 0-9997 0:9973 0-9785 1-828 1-991 
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where the skewness coefficient y is given by 


y= 2 %(e)-/| 3 a,| ue LO 


To facilitate the computation of y, the values of the dimensionless 
quantity 

Px3—=2—U2Q(1+-@)/P? 
are also tabulated in table 1 (as [«3],, [x3], being defined in § 4). To 
the same order of approximation, since S=-p to the first approximation, 
the equation for the limits may be alternatively written 


S—=dy(p?—l= ee ne 
The expression (9) has the advantage of being available for any p; 
(11) has the advantage that if interpolation for S is approximately 
linear, e.g. on the 1/6 scale, the left-hand side of (11) will probably be 
rather more so than (9). 

To illustrate the use of these various methods, the data of Fretter et al. 
(1953) on neutral V-particles have been used. These authors give for 
their second type of particle V,° an estimated mean life-time of 
4-3 10-1 sec (based on eleven observations), their estimation being 
in effect based on the eqns. (2) and (3) above. We find, however, from 
the crude but simple method (6) that 


1 
q~ 0-21250-15 x 1019 sec-2 5 oy See 


indicating that a value 1/9=0 is quite consistent with these data (a 
conclusion not apparent from the standard form of maximum likelihood 
estimate used by the authors). A more exact use of (4), using the 1/0 
scale, shows that interpolation on this scale is nearly linear (see fig. 1). 
Weare therefore more justified in quoting the maximum likelihood estimate 
for 1/0, this being 0-256-L0-198. Since the value of J, —- tT Ve ae 
is only —1-39, a finite upper 0-05 confidence limit for @ is not possible ; 
for the lower 0-05 limit we obtain 


1 
H ~~ 0-587K 100. een ee 
6 0:05 


(compare the approximate value 0-582—0-256-+ 1-645 x 0-198). With the 
skewness correction (the use of both the alternative corrections (9) and 
(11) is shown in fig. 1), this is modified to 


1 
H AO0-621501020, as Ieee 
6 0-05 


It may be of interest to note that the crude method (8) gives for the 
data on V,° (23 observations, with @ estimated by Fretter et al. to be 
10-+7 X 10~?° sec) the estimate similar to (12) 


> 


1 
G7 00044. 0-086X 108 sect ees 


again compatible with 1/9=0. With such data it is quite evident that 
often only lower bounds to 6 can be given. 
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§3. THe Errect or Momentum Errors 


One of the main sources of actual error in the measurements t¢, 
and 7, is in the conversion factor from the corresponding length 
measurements d, and D,, say. As the conversion factor is the same 
for both measurements, it has been wondered whether a more accurate 
method of estimation might be obtainable by working as far as possible 
with the ratios d,/D,=t,/T,. This question is examined briefly here, 
with reference to the type of problem discussed in §2. The nature of 
the momentum errors will of course differ in different circumstances, for 
example, with charged particles in a magnetic field more direct determina- 
tions are possible than with uncharged particles ; we shall merely consider 


Fig. 1 


Example of S(9) plotted against 1/@ (eqn. (4) ; the lower continuous line 
corresponds to the skewness correction of eqn. (11), the dotted line to the 
alternative use of eqn. (9)). 


first-order effects when the error in 7’, is assumed to be of the fairly 
general form 
OTe Tete eT Ee Mase fe: Wey es. (16) 

where the first component is a bias or systematic error, and € in the 
second component is purely random with mean zero and variance oe 

In the eqn. (2), we see that each 0L,/00 is weighted equally in the 
complete equation. The use of t,/T,, in place of t, will still permit (in 
the absence of momentum errors) an asymptotically unbiased estimate 
of @ if we write alternatively, corresponding to weights 1 ales 


1% fu,—l = | 
= is ais ==(); ° ° e . ° (2 a) 
ve OF SEP: 
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We denote the expression to the left of eqn. (2) by A and the corresponding 
expression in (2a) by B. Then we find 


Ube sae ie 


Consider first the resulting bias i in the eee of @. Since 


E(8A)= a VEPs) ? 


r=1 


TART a Ud he AW I ? 


we find that the proportionate ae in Ms estimate of @ (or to the same 
order, of 1/0, etc.) is vXT'.°-1(L,),/2(1,), for 6 and v3T,2-(I, )/AUL Bh hed he 
for the estimate obtained nae (2 - If we put «=1, this gives the 
trivial but correct answer v in both cases, for all measurements are then 
wrongly scaled and the final estimate must also be. A more relevant 
assumption in many cases is that «=0. This corresponds to the 
approximate relations 
T,Cd,/Pyy Py XPy's 
where p, is the estimated momentum and 5p, its error, the second 
relation 5p,op,2 being probably representative at least of the errors in 
momentum determinations from the track curvatures of charged particles. 
In such situations we should have «=$8=0. This gives biases by the two 
methods of 
pel ee a Voile aly 
211), 3 rd ee) ; 
By Schwarz’s inequality, the first is always less than the second, so that 
eqn. (2a) appears in this respect at least to be inferior to (2). For 
example, if we consider the same range of 7', as in the sample of eleven 
observations referred to in § 2, we obtain biases 0-136v, 0-249v respectively. 
In the case of the random error component, the position is not quite 
so clear-cut. The total error variance of our estimate will be made up 
of two parts, one due to the smallness of the sample and the other to the 
random momentum errors. By the usual approximate methods of 


evaluating these errors we find for the overall error variances of our two 
estimates 


(17) 


2(I,),+07(5A) ade hi pact (OD) 
mPa N mae ee eipvircemperae cage o. SS? 

where 

o°(8A) woPZ ST, [(1,)(1,+1)],; 

o*(8B) wo PZ 2T 2P-4[ 7.2). 
When o,” is zero, the second method is of course inferior to the use of the 
maximum likelihood equation. For non-zero o,2, however, it is possible 
that in some cases the second overall error variance would be less than 
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the first, though this seems rather unlikely. For instance, for the range 
of 7’, in our previous example and B=0, we obtain as our two approximate 
error variances 

0-629(1+0-0375¢2),  0-629(1-83-+0:00935 o,2), 
where 10~1° sec is our unit time interval. The second error variance is 
only less than the first if o,>5}, which for 7, ranging from 1 to 20 
represents a very large random error. 


§4. OnE TyPE or PaRTICLE, witTH A RECORD OF UNDECAYED PARTICLES 


It may sometimes be possible to take into account tracks of charged 
particles which have not decayed. For such data we thus assume particles 
1...n as before, with n+1...N further undecayed particles. As the 
probability for the sth undecayed particles is Q,=e-7s’°, and the 
probability for the rth decaying particle with decay in the range 
t,, t,+dt, is e~'r’? dt,/0, we now have for the total data a probability p 
such that 


0 log p ia. Wi a 
ee a + | 2 tt ee | eee arri9} 
so that the maximum likelihood estimate now becomes 
ai l n N 
6=-| Xt4+ 2 al M, Pere eee, 
n Koa s=n+1 


It may at first sight appear a little surprising that this estimate is simpler 
than that for decayed particles only, but it will be seen that the more 
complete records permit a direct estimate of the mean life-time to be 
made (cf. Astbury et al. 1953). The information function corresponding 
to (19) and (20) is 


N 
Ce me es 21) 


say. 

D graph of P=I,, the information (apart from the factor 1/0") per 
particle, is plotted against U in fig. 2; it represents the fraction of the 
information 1/62 available per particle in a chamber of unlimited size. 
For comparison, the fraction PJ,, expected if undecayed particles cannot 
be recorded, is also shown. For small U, it will be seen that the ratio 
PI,/I,=1, is very small, so that there is then a large gain in information 
if the more complete records are obtainable. 


From (19), we have 
dlogp _® Cay 22 
6 AC 1) ae SORA a oe 


with mean zero and variance X(I,),. Hence a confidence interval may 
be obtained, on the assumption that 0 log p/00 1s approximately normal, 


from the equation 
A N 
n(5—1)/ 9 2 Coes a) ge (28) 
Ne ran 


[ 


———---—— ————-— — 
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As a further approximation we might substitute within the square root 
the estimate n, giving confidence limits for 0 of 


Si(izpia/n; . ee 
though this will not be advisable if n is too small. It is interesting to 
notice from (24) that here also 1/@ is a more useful scale to use than 0. 

If a correction for skewness is required, we find from the general 
formula (7) 
ut 2s 
Z (Ks)o= 7s 2 (2Pg—3Q Ua + - » - + 0) 
q=1 q=1 
and make use of a modified equation similar to (9) or (11). The 
dimensionless quantity 6°, is tabulated in the table (as CAA while 
negative for a small range of U, it is in general positive. Of course in 
this type of problem n is one of the random quantities, and, being discrete, 
will make precise approximation impracticable for very small n. 


Fig. 2 


The expected information for decaying particles, and the information for all 
particles, expressed as ratios of the information available as J’ > oo. 


If mixed data are available, some being of the type considered in § 2 and 
some of the present type, we merely add the eqns. (2) and (22), and the 
total information function will be the sum of the separate ones; 
similarly, if required, for the third moment. 


§5. Two TypsEs oF PARTICLE 
There is no difficulty in principle in extending these methods to two 
types of particle in unknown proportions, but a general investigation of 
accuracy divorced from actual data becomes laborious and it seems 
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preferable to proceed directly in any problem, using if necessary such 


approximations as 
_ & log p 2 log p 
0a0B —— 0a0B 


where p denotes the probability of the entire data. Thus for data of the 
type considered in § 4 we have 


p=ll {D4 Px+WaPp}s CARO See a O (26) 
where p,=1—q,, assumed constant for all particles, is the probability 
that a particle is of type a with decay constant «; f is the other decay 
constant. (Notice that p, can only be constant relative to the observa- 
tional conditions, and will depend on «, 6 and the distributional origin 
of the particles relative to the (effective) initial boundary of the chamber.) 
Px (pg) denotes the probability of the observation given that the particle 
is of type a (b) and is of the mixed type given in § 4 (i.e. it is an infinitesimal 
probability only if decay occurs, in which case the infinitesimal element 
dt, is dropped). JI denotes multiplication over all particles, the dis- 
tinguishing suffices for which have for convenience been omitted. The 
maximum likelihood equations for p,, « and f are 


Ologp _ 5 _Pa~Pp 


ee 


=i) 
Oda PaPatTCaPp 
0 log P 3 pals 0 log Pa =0, «ee ys (27) 
TEE PaPaotIaPp 9 
Glog p _ »__GaPs 9 108 Pp _y 
a PaPatdaPs OB 


where @ log p,/d«, Alog p,/@B are as in eqn. (19). The approximate 
standard error of any estimate is the square root of the corresponding 
diagonal term in the inverse of the matrix of terms Tee, Laas CC, With! 
respect to the three unknowns «, f and p, (see, for example, Kendall 1946, 
ch. 17). It is likely to be large in most cases, unless the data are very 
extensive. 
In the particular case where one decay constant, say P, becomes large, 
a simpler procedure becomes available. The chance of decay in the 
chamber of the second type of particle is now negligible, so that (unless 
Pq is nearly zero) the problem reduces to the estimation of a by the 
method of $1, and then the estimate of p, from a knowledge of the 
chance of decay P, for each of these particles, the ratio of decayed to 
undecayed particles being equated to 
N N 
Se ea ae (v. a Ort Naa) ihe Ae AV ee SS) 
q=1 C= 
In the case where only decayed particles are available, the form of 
eqns. (27) is modified, because the chance of such a decayed particle 
being of type a is no longer p,, but 


Da’ =DoF (a)lPaP(a)+daP(B)], - - + + + (9) 
depending in a variable way on « and 8. Equation (26) now holds with 
p,' in place of py, f(t; «); defined as in eqn. (1), for p, and f(t,; ) for 
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p,-* In the particular case where f is large, we must assume that Pa is 
not near unity for typical 7’, (otherwise no estimation problem involving 
more than one unknown occurs), and this now implies that p, must be 
small. Equation (26) becomes replaced by 


p=T pa (f(t,; )—-YT,)+UT,}, 


eo DEE ae i 
weer Pa = 5 (PAa)—TyBy+Pb pPAa)+2, ee 


where p= p,f is assumed the new unknown parameter. Hence 


p=IT{[pP,(a)f(t, 3 «)+UlleP (oe) +7] 
or 

16g p= Slog (pf--1)—log (pP,-P Ti), eae 
where f,=(e-*r/%)/a, P,=1—e-*®, 

In view of the difficulty even with only one type of particle in rejecting 
@=co from records only of decaying particles it seems very doubtful 
whether for such data we should feel justified in going beyond (31) to 
the full equation involving three unknowns. Even with (31), unless, in 
spite of the further unknown p, we obtain an estimate of « clearly distinct 
from ©, it is evident that the simpler hypothesis of only one mean life-time 
is adequate. In any case the interpretation needs a little care, for the 
assumption of 8 large gives a uniform distribution of decays from this 
type of particle over the chamber, and could be simulated by some 
mis-identifications or general ‘ noise’ effect in the chamber. 

From (31) we obtain the estimation equations 


0 log p a. dL Pres! a 


—_—-= — Fy —___ =( 

Op (ef, +1)(pP,+T,) i (32) 
0 log p =b 3 [Oe U,Q, =O 

Oc. on iG oat pie ee ta ee hides 


Instead of using maximum likelihood estimates obtained from (32) with 
approximate standard errors (in which case it would seem advisable to 
use the 1/x scale), we could extend the method of confidence limits to 
give a confidence region for p and « simultaneously. This would consist 
of constructing a region of values of p and « such that 


Janae +2 gonad ope SN: 1 ini ae 


a aa oe (; an a =) xe 
J x9 Ue I “o ts Lon, 
* Other modifications may be required in particular circumstances, depending 


on further relevant information from the track records. For example, P,(q) 
may be interpreted somewhat more generally as the chance of the observed 


where 


decay in the chamber (apart from the value of t,), given that the particle is of 


type a, and similarly for P,(f), (in this case the later eqns. (30), (31), (32) will 
need reconsideration). If, further, the conversion factor depends on the type 
of particle, then t,, Z', will no longer have the same values in f(f,; «) and 


Stes B) 


cae 


Statistical Estimation of Mean Life-Times 261 


the suffix p being used in connection with the unknown p, D, = 0 log p/d«, 
D,= 0 log p/dp, and x,” is the critical value of a ‘ x? quantity’ with 
two degrees of freedom, e.g. 4-6 for the 0-10 significance level (see, for 
example, Fisher and Yates 1938). As the quantities D,, D, are squared 
in (33), it is less necessary here to consider any skewness corrections. 
Of course, if this region ‘ explodes’ with respect to large « for all p, the 
hypothesis of two distinct life-times is again hardly justified. 


§6. DISCRIMINATION TESTS 


The difficulty of separating a mixture of two types of particle implies 
that a complete statistical analysis is hardly worth while unless the data 
are rather extensive. In this and other cases a comparison of the fit 
obtained by two different hypotheses may often be useful, and this is 
usually more conveniently performed on the L functions (log. probabilities) 
directly rather than on their derivatives or on the estimated constants. 
Thus let L=log p be the probability of (26) (with any infinitesimal 
elements omitted), and L’=log p’ the probability on the hypothesis that 
only one type of particle (in regard to life-times) is involved. Then we 
may take very approximately 


—2( Dinas qaae)> 


as a x? quantity with two degrees of freedom (cf. for example, Kendall 
1946, § 21.20), where the maxima are with respect to variation in the 
unknown constant(s). Unless we can find values of p,, « and f to make 
the quantity in (34) greater than, say, 4-6 (the 0-10 significance level 
of x? for two degrees of freedom), the more complicated hypothesis would 
hardly be justified. (In such a test it is essential that the most efficient 
estimate of the decay constant in L’ be used.) 

For an equation like (31), the corresponding approximate ° would 
have only one degree of freedom. 

Another possible application of the x? procedure is for testing the 
compatibility of two or more sets of data. Thus for two sets supposedly 
referring to one and the same type of particle we should have 


—2({L, Imax + [Lelmax—[L1+L2hmax), = slp 10 (34) 


as a x? with one degree of freedom, where L, and L, refer to the two 
sets of data. Subdivision of the same data, not involving selection of t,, 
e.g. with regard to values of 7',, might also sometimes be useful for testing 
internal consistency. 
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SUMMARY 


Measurements have been made of the damage caused by «-particle- 
irradiation of several materials; some results are given for electron and 
ultra-violet irradiation. A comparison of the fluorescence efficiency for 
%-particle and electron excitation for these materials is made. The 
implication of the results is discussed in some detail, and analysis leads to 
a new understanding of the abnormally low «-particle fluorescence yield 
in organic materials. 


It is known (Birks and Black 1951) that the fluorescence efficiency of 
an anthracene crystal decays during the exposure of the crystal to a 
strong «-particle beam. Investigations extending to other materials and 
other forms of bombardment are described briefly below and shown to. 
lead to a new view of the mechanism of fluorescence. 

The six organic materials used were in crystalline form and _ had, 
whenever possible, freshly cleaved surfaces. Sources of 7!°P and *5S were 
employed to provide beams of «-particles and electrons. The decay in 
the fluorescence intensity during the bombardment was found directly 
from measurements of the output current of a photomultiplier, changes. 
of 3% being detectable. 

The variation of the fluorescent intensity L under constant excitation, 
with the number N of «-particles absorbed per cm? was found to be 
similar for different organic materials, being expressed over a range 
depending on the extent of variation in energy and direction of the 
a-particles striking the crystal, by the equation 

AN ie acer oie men te CE) 

A is a constant, characteristic of the material. Small variations in the 
value of A for anthracene were found when different faces of a cubic 
crystal were presented to the «-particle beam. This effect might have 
been more marked had the beam been collimated. 

Some of the experimental results for «-particle bombardment are 
shown in fig. 1 and in the table. £4), is the energy dissipation per cm? 
required to reduce the fluorescence efficiency to half its initial value. 
No figures are given for CaWO, or ZnS-Ag. Experiments with these 
materials failed to reveal any perceptible change in fluorescence for 
Wa=45 101. 


* Communicated by Professor P. I. Dee, F.R.S. 
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Experiments with electrons as the bombarding particles were confined 
to anthracene. No change in fluorescence was detected, whether measured 
by «-particles or electrons, during a period in which the crystal absorbed 
1018 electrons/em? of average energy 40 kv. This was equivalent to an 
energy dissipation near the surface of 1-5 108 erg/em’. Irradiation with 
electrons thus appeared to have no effect upon the fluorescence comparable 
with that produced by «-particles. 


Fig. 1 
x 
x 
4 x 
A 

3 

a 

as 10 (Saeco € 
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The variation of Z-1 with N for anthracene bombarded by 3 Mev a-particles. 


Relative figures for various materials were obtained for the ratio, R, of 
the fluorescence efficiency for excitation by 3 Mev «-particles to that for 
excitation by electrons of average energy 40 kv. Values of R, relative 
to that for anthracene, are shown in the table. 


Table 


Anthra-| Naphtha-| Ter- | Stil- | Py- |Diphenyl}/ZnS—Ag|CaWO, 
cene lene |phenyl} bene | rene [acetylene] powder | crystal 


Fy), 10-8 
erg/cm? : . Bee 7:0 13 32 — — 


Ax1012 i : 37 2? 0-8 Os) ee Od 


R . . : : 5 2 0-7 24 3:2 


Irradiation by light of wavelength 3650 A was found to reduce the 
fluorescence efficiencies of anthracene, naphthalene and stilbene. The 
values of H,). in erg/em*® were very approximately 8 x 1011 for anthracene 
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and less than 4X 10!° for naphthalene. For naphthalene and sfilbene 
the decay was only temporary and the efficiency increased again when 
the crystal was left unexposed. Anthracene did not recover even when 
heated to 100°c for some hours. ‘The fluorescence efficiencies were 
measured by means of a standard weak 22°P gource. A preliminary 
experiment had shown that when an extremely thin anthracene crystal 
is damaged by «-particles, the efficiency of the crystal for various forms 
of excitation (electron, light and a-particle) all decrease by approximately 
the same factor if this factor is small. 

_ The relationship between L and N holding for a-particle bombardment 
and expressed in eqn. (1) was found for anthracene by another method 
by Birks and Black (1951). It also seems to hold for several inorganic 
phosphors according to calculations made from data given by Marsden 
(1910) for «-particle bombardment and by Grotheer (1939) for electron 
bombardment. This fact was not noted in the papers concerned. 

The results for the inorganic phosphors, however, were complicated by 
loss in transparency of the crystal during bombardment and by a recovery 
of the efficiency. These effects for the organic materials were too small 

_ in most cases appreciably to affect the values of L over the ranges of N 
for which the experiments were conducted. Hence for the organic 
materials the efficiency W could be expressed, relative to its initial 
value by 

Ve (1a) ere ecke Cee? eirce Nets 22(2) 


where [D], assumed proportional to N, is the concentration of damaged 
molecules in the bombardment layer and C is a constant which for 
anthracene cannot possibly be less than 400. 

As was noted by Birks and Black (1951) this result is to be expected 
if some form of intramolecular energy transport occurs, the excitation 
energy moving through the crystal until it is converted into fluorescent 
light or, on meeting a damaged molecule, into heat. Ifthe energy spends 
on the average a time @ on each molecule of a material of decay time 7, 
containing a concentration [D] of damaged molecules, where one damaged. 
molecule captures the energy with a probability & relative to that for an 
undamaged molecule, then the efficiency W is given by 


We IEUDY) | iGaty dh ees get (3) 
J0 


where y, the fraction of the energy initially created which survives tor 
a time ¢, is determined by the equation 

o =—y(7-1+kO-1[D)). rer ws) (4) 
When @ is independent of t, eqn. (3) reduces immediately to eqn. (2) with 
Taeent reveals that the power of prolonged «-particle henbaeiaent 


to produce a reduction in the fluorescence efficiency varies widely between 
different materials. Yet it shows also that the fluorescence efficiency for 
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excitation by «-particles is smaller than that for excitation by B-particles 
by a factor which is approximately constant for all the organic materials. 
This suggests that the saturation of the specific fluorescence for «-particle 
excitation found by Birks (1950, 1951) may not be caused by those 
damaging or fluorescence inhibiting effects which a single «-particle 
produces locally and permanently and which, when caused by large 
numbers of «-particles, become measurable as macroscopic reductions in 
the fluorescence efficiency. 

We are therefore forced to look for some other explanation of the 
abnormally low fluorescence yields of organic materials under «-particle 
excitation. Furst and Kallmann (1952) have suggested that in the 
a-particle track interaction between two neighbouring excited molecules 
may lead to non-radiative dissipation of the energy. It is known (Franck 
and Livingstone 1949) that energy can be transferred between molecules 
for distances of up to 60 A. This part of the process generally described 
as sensitized fluorescence will be most probable when the fluorescence 
spectrum of the donor molecule coincides with the absorption spectrum 
of the acceptor. This increased probability is more readily achieved for 
two excited molecules than for molecules only one of which is excited. 

The result of energy migration between two excited molecules is that 
one is reduced to its ground state while the other is raised to a higher 
excited state. A direct radiative transition from this state to ground 
does not seem to occur for anthracene (Pringsheim 1949). Instead the 
high excitation energy will probably lead to various photo-chemical 
processes such as dissociation, the net result being failure of the molecule 
to emit except in the infra-red. This process of double de-excitation 
furnishes a basis for the understanding of the behaviour of the specific 
fluorescence of organic materials when excited by ionizing particles. 

Let an ionizing particle raise m) molecules per unit length into their 
first excited state. The number n of these excited molecules remaining 
in the track at a time ¢ after the particle has passed will be given by 


dn 2 a 
Tp Ve — pn pn(1—e- *) io te ees 


The first three terms on the right-hand side represent respectively the 
probabilities of fluorescence, double de-excitation and energy migration 
outwards from the track. Not all the energy which proceeds from the 
track will lead to fluorescence as a fraction of it will, before it can escape 
be lost by double de-excitation with the main group of excited molecules! 
It can be shown that this fraction is (l1—e-*) where X is a constant 
less than «. The fourth term accordingly must be included. The specific 
fluorescence F is then given by 


P= |. {yn-+un—pn(l—e-™)} dt... 6) 


Briefly, i = dn, an =— (yN +pne 4”), 
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whence we obtain from (5) 
that y <p, dF 
dn, 2 ree—1}4, eee ee ves (' 7) 


on the assumption, almost certainly true, 


where z=An,y and a=«/du. 

Equation (7) was integrated graphically and XF as a function of z is 
shown in fig. 2 for various values of a. The scales were adjusted so that 
in each case the maximum value of AF was the same. The shape of the 
curve is independent of the value of a when this is less than 1. The points 
representing the experimental results of Taylor et al. (1951) for two 
kinds of heavily ionizing particles all lie within 20% of the curve for 
a—10. From a figure of 5°% for the efficiency of anthracene for electron 
excitation, m) can be calculated to be about 2x 107cem for a 5 Mev 
a-particle. This with a=10 gives A=7x 10-8 and “a=7X10-"p. This 
is for anthracene. Other materials will have different values of A and pw 
and consequently the shape of the F curve will vary slightly. Taylor 
et al. have shown that this is so for anthracene and stilbene. 


Fig. 2 


0: 


oie 0-1 ers as fo) 
The variation of F with Z. Curvel : a<1; curve2: a=10; curve 3; a=—100. 
" TThescale exhibited isfora=0. The crosses represent data given by Taylor 
for two kinds of heavy ionizing particles in anthracene, on the assumption 
that, is proportional to the specific energy loss. Thecircleisfor electrons 
of 120 kv. 


I am indebted to Dr. 8. C. Curran for many helpful discussions on the 
work and to Mr. G. M. Lewis for his advice and encouragement. 
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SUMMARY 

The track-density of the Bristol «,-meson has been determined photo- 
electrically. The measurements have been accurately corrected for 
uneven development and transparency of the emulsion. The mass has 
been calculated by comparison with tracks which have roughly the same 
slope as the meson-track and also by comparison with a proton track 
which crosses that of the meson. The result of the mass-determination 
is mk,=1518-+70 electron masses. 


§ 1. INTRODUCTION 


A HITHERTO unobserved phenomenon was recorded in a nuclear research 
plate which the Bristol group exposed on the Jungfraujoch under 30 cm 
lead (O’Ceallaigh 1951). A particle whose track is reminiscent of that 
of a proton, as far as grain-density and scattering are concerned stops 
in the emulsion and decays with the production of one charged and one 
or more neutral particles. Small-angle scattering measurements on 
the track of the charged secondary give pBc=250 Mev, while the grain 
density is the same as that for particles at minimum ionization. These 
facts set an upper limit of ca. 400 m, to the mass of the secondary particle, 
which, if it is of known type, must thus be an electron, a p- or a 7-meson. 

The range of the primary in the 400 » thick Ilford G5 emulsion 
employed is more than 4000 y, long enough to allow of a satisfactory 
mass-determination. Two estimates made by O/’Ceallaigh, the first 
depending on combined range and scattering measurements and the 
other on a combination of gap-counting and range observations, gave 
the results (1260+290) and (1385+200)m,, respectively. Though 
the agreement is satisfactory, it nevertheless seemed desirable to make 
a third estimate by photo-electric measurement of the track density. 


§ 2. THE MEASUREMENT OF THE DENSITY 
The photo-electric apparatus has been fully described by von Friesen 
and Kristiansson, and it has not been altered in the present work (von 
Friesen et al. 1950, 1952). All measurements were made with the same 
slit, its dimensions in the focal plane of the objective being 1-98x 30 
microns. Density observations were begun 20 « from the end of the 


* Communicated by Professor Sten von Friesen. 
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track and continued to the point where it left the emulsion. The 
density at a given point on the track was taken to be the mean of the 
densities in two equal intervals one slit-length long, one at each side of 
the point, that is, the density was determined at points 30 pe apart. 

In fig. 1, where the *;-meson is compared with one p-meson and 
three protons, the integrated density is shown as a function of the range. 
Particle No. 1 is especially interesting since it is a recoil particle, and 
hence it can be assumed with very great probability to be a proton. 


Fig. 1 
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Integrated density versus residual range for the x,-meson, three protons 
and a p-meson. 


It is situated at approximately the same place as the «,-meson track, 
with which it is roughly parallel. Thus it is conveniently used as a 
standard of comparison in determining the x,-meson mass. Particles 
Nos. 2 and 3 are also protons, as will be shown later, and their directions 
are not very different from that of the x,-meson. Figure 2 is a vertical 
section through the emulsion showing the position of the particles. 
When diagrams Nos. 1 and 2 are compared one sees that the values 
obtained for the integrated density is bigger the higher up the particle 
lies in the emulsion. Hence one cannot exclude the possibility that the 
different densities exhibited by the particles (fig. 1) depend more on the 
position of the particles in the emulsion than on possible differences in 
mass. To find out if this is the case the density was determined for 
a number of tracks which lie at different depths in the emulsion. In 
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this way a correction factor has been computed for each depth,fand by 
multiplication with this factor each density measured has{ been reduced 
to the middle depth in the emulsion. 


Fig. 2 
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Vertical section through the emulsion showing the position of the «,-meson 
and of the three comparison protons. 


§ 3. DETERMINATION OF THE DEpru-CORRECTION FOR AN EMULSION 

If particles with the same velocity and charge traverse a thick emulsion 
at different depths their tracks will not have the same appearance 
although the ionization is the same. The projected track area registered 
by the photo-electric apparatus is found to be a function also of the 
position of the track in the emulsion. This variation depends principally 
on three effects. 

Development 

It is practically impossible to produce absolutely even development 
throughout the whole of a thick emulsion. Irregularities are chiefly 
of two kinds: firstly the number of grains is variable, and secondly 
the average size of the developed grains can be different. Only the 
first factor comes into play in grain-counting, while both are of importance 
in photo-electric recording and gap-measurements. In the emulsion 
containing the «,-meson an additional fixing effect occurs. During 
the fixing process precipitated silver was dissolved, especially in the tracks 
nearest to the surface of the emulsion, and the area of the tracks was 
diminished. 

Fog 

Tracks which are not completely black are observed in thick emulsions. 
This is especially so for tracks which are situated deep in the emulsion. 
The effect occurs because the light is diffracted both by very small 
background grains formed during development and by grains in tracks 
which lie higher up in the emulsion. The result is that the track also 
appears to have a certain luminosity and light which apparently originates 
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in the grains measured. In these circumstances the calculated track-area 
becomes less than the true area. 


The Microscopic Effect 


The microscope probably introduces an effect also. The appearance 
of the imaged light-beam varies somewhat with the position of the track 
in the emulsion because the distance between the objective and the 
condenser is different when tracks at different depths are studied. The 
effect is undoubtedly small and possibly totally negligible. 

The degree of development, the fog, and the microscope effect are all 
functions of how deep the track lies in the emulsion and of that only. 
Hence the correction which will eliminate them and reduce all track- 
densities to a common depth must also be a function of the depth only. 

Fog and the microscope-effect depend on the emulsion and the 
microscope, respectively, while the variation in the degree of development 
is directly connected with the track under investigation. It is therefore 
not inconceivable that the degree of development might be a function 
also of the track density. The total correction would then be different 
for tracks of different densities. I have therefore determined the 
correction factor for three different densities. 

Protons, deuterons and tritons have been employed for the best 
determination of the correction factor. Their tracks have been studied 
in the range R=0 to R=1000 py, in which interval the density can to a 
good approximation be considered constant and identical for all three 
types of particles with the degree of development used. In fig. 3 the 
observed densities are set out as a function of depth in the emulsion. 
By dividing up the emulsion into equally thick layers and calculating 
the average density in each, a best-fit curve for the given degree of 
development has been drawn. The errors indicated in the diagram are 
standard errors. 

Eight »—mesons have been employed in determining the development 
correction for less dense tracks. In these p—meson tracks the region 
R=0 to R=600, has been used in part, as has also the interval 
R=600-1200 , to calculate two correction curves. In these cases it 
is impossible to plot the points directly in the same way as before because 
the density varies in the interval. Instead, the correction has been found 
by successive approximations. Figure 4 shows the relation between 
density and depth for the ionization densities in question. In certain 
regions the corrections are less exact than the proton correction because 
of the smaller number of points used in determining them. (The curves 
are dashed in these regions.) 

The calculated correction factor for the different ionization. densities 

is given in fig. 5 as a function of depth in the emulsion. This 
correction factor is defined as that factor by which the measured density 
must be multiplied in order to reduce it to half-depth in the emulsion. 
One sees that there is a little difference in the correction for the different 
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Density of proton-, deutron- and triton-tracks versus distance from the glass. 


Fig. 4 
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Development curves for low densities computed from p-meson tracks. 
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ionization-densities. It is not possible to tell from the measurements 
made whether or not this difference depends on the ionization or whether 
its origin is purely statistical. 


Fig. 5 


Correction factor 


Calculated development-correction factor versus distance between track and 
glass. 


Fig. 6 
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Corrected integrated density versus residual range. 


All the tracks employed in determining the mass of the «,-meson 
have been reduced to the centre of the emulsion. The integrated densities, 
after correction, are given as a function of the range in fig. 6. This 
is to be compared with fig. 1, in which the density before correction 
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is plotted. The proton-group in fig. 6 is an average derived from four 
protons. The separate particles lie so close together that it is not se ye 
to resolve the individual particles in the diagram. 


Fig. 7 
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«,-meson density, before correction, versus residual range. 


Fig. 8 
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Corrected x,-meson density versus residual range. 


Figures 7 and 8 demonstrate how important the correction is for a 
single particle. The «,-meson is chosen as an example. In fig. 7 its 
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density before correction is shown as a function of range, while the 
corrected densities are given in fig. 8. It is seen that the measured 
values have a large spread. The following method has been used to obtain 
best results in drawing the density-range curve, which can be approximated 
by a straight line. The track is divided up into intervals for each of which 
an average density value is calculated and inserted in the diagram. 
These points are then used in calculating the constants of the straight 
line, according to the method of least squares. The portion of the track 
closest to the point at which the particle enters the emulsion is not 
employed in the calculation since the density there may have been 
influenced by unconsidered effects in the surface of the emulsion. The 
density-range curves were drawn for the other particles in the same way. 


Fig. 9 
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Range-density diagram for the x,-meson and the comparison particles. 


§ 4, DETERMINATION OF THE Mass 


In the density-range diagram of fig. 9 curves are shown both for 
the «,-meson and for the particles with which it is compared. in finding 
the mass. The proton curve is a median curve of four protons of which 
three travel approximately parallel with the «,-meson, the fourth 
being a particle which goes from the glass towards the surface, where 
it stops. The »—meson curve is a median of the eight particles used 
in calibrating the plate. 

In table 1 the data will be found for those straight lines which 
approximate the density as a function of range. The slopes and their 
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intercepts on the axis R=0 are given. The figures give an indication 
of the accuracy with which a density measurement can be carried out. 


Table 1 


Density 
at R=0 


Particle Slope 


k,-meson 59-9 
Proton 1 59-5 
58-0 
59:7 
A 59-1 
Deuteron 59-9 


A relation R=(m/Z?) f(v) obtains between the residual range, mass, 
charge and velocity of a particle. Consequently if two particles of equal 
charge are compared with one another the ratio of their masses is that 
of their residual ranges provided that the comparison is made at points 
where their velocities are equal. Two particles with the same velocity 
and charge produce the same ionization, and therefore their track-densities 
should be the same also. . In order to obtain the relation between the 
masses of two particles it is thus sufficient to determine the relation 
between their residual ranges beyond points of equal density. 

The mass of the «,-meson can thus be found by comparing the 
proton- and «x,-curves in fig. 9. When the comparison is made for 
a density of 52 units the result is mx,;—1550 electron masses. This 
value is not, however, independent of which density is chosen in making 
the comparison, since the proton- and meson- curves intersect. The 
- statistical fluctuations in the densities of the individual tracks are so 
great that such an intersection is not improbable. The comparison 
is made at as great a range as possible because the relative error in the 
range-determination is there least. 

However, it is possible to eliminate the crossing of the proton- and 
meson-curves. The mass determination is then independent of the 
density chosen for the comparison. The average value of the intercepts 
of the curves on R=0 has been found and all the curves are made to pass 
through this point. The slopes are found by a new least-squares 
calculation, this time under the proviso that all curves should cut the 
line k=O at this point. It is theoretically correct that all the curves 
should cut the axis of ordinates at the same point since all density-range 
curves for singly-charged particles are similar in shape, with their 
f-axis elongated in a degree proportional to the mass. The intercept 
on the line R=0 is 59-4 units. The slopes are shown in table 2, which 
also gives the value of the «,-meson mass calculated from each individual 
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comparison particle. The mass-relation is aly found since the mass. 
is inversely proportional to the slope. 


Table 2 


Particle Mass of the 


k,-meson . 
Proton 1 . 1549 m, 


1533 

. 1503 

» 4 7 1509 
Deuteron ti 1512 


Mean value a 1521 


When all four protons and the deuteron are employed, the average 
mass of the x,-meson is found to be (1521+75)m,. The error is a 
probable error determined from the spread of the points which have 
been used in the calculation of the density-range curves. The error in. 
the slope 4k has been calculated from the formula : 


een, | 


where w; is the statistical weight of the different points, 4D, is the 
difference between the density measured and the density determined 
from the density-range curve for the range R;. The statistical weight 
of the points has been put proportional to the range. 


§ 5. SysteMATIC ERRORS IN THE Mass-DETERMINATION 


When the p-meson mass is calculated from that of the proton by 
comparing ranges a result ca. 20% too small is obtained. The same 
observation was made in earlier measurements with the photo-electric 
apparatus (von Friesen et al. 1951), and if the effect is the result of some 
systematic error peculiar to the method of measurement then one must 
assume that the same systematic error arises in determining the mass. 
of the x,-meson. The correct masses of the particles are shown as. 
abscissae in fig. 10, the ordinates being the measured masses, The 
departure from the correct value (represented by the dashed line) is large 
in the case of the p-meson but should approach zero asymptotically, 
and is, therefore, vanishingly small for the calculated «,-meson-mass.. 

In determining the meson-mass given in table 2 the density curves 
have been approximated by straight lines. This introduces a systematic 
error in so far as the measured ranges are not proportional to the masses. 
If such a proportionality did occur, the curves with their approximating 
straight lines would have been of similar shape, and the approximation 
would have been of no importance in the mass-determination. An 
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attempt to estimate the order of magnitude of the error was made by 
computing second degree range-density curves using measurements on 
protons and mesons. The curvature obtained, however, is so small 
that the standard error in it is of the same order as the curvature itself. 
Thus a reliable determination of the systematic error introduced by the 
approximation is impossible. An estimate has, however, been made 
which shows that it must be small in comparison with the probable 
error calculated from the spread of the density-measurements. 


Fig. 10 
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Comparison between the measured and actual masses of the comparison 
particles. 


It is not possible to decide what kind of systematic error is brought 
in by irregularities in the emulsion, though there is no reason to assume 
that any error arises because the «,-meson is quite far from the edges 
of the plate. Further, as already has been pointed out, proton No. 1 
is at the same place in the emulsion, and it shows no sign of anomalous 
density. 


§ 6. Mass-DETERMINATION BY COMPARISON WITH CROSSING PROTON 
TRACKS 

The «,-meson mass-determination just described was carried out 

following an appreciable correction of the density, and the accuracy 

therefore depends to a high degree on how good the correction is. It is 


oe ee aay 


Determination of the Mass of the Bristol K,-Meson 279 


possible, however, to find approximately the relation between the masses 
of two particles without considering exactly how uneven the development 
is, provided that it is a function only of the distance from the glass. 
Instead of comparing tracks which coincide as nearly as possible with 
one another in a vertical section of the emulsion. tracks which Cross 
each other are employed. The «,-meson travels from the surface 
through the emulsion and stops near the glass. It is compared with a 
proton which comes from the glass, goes upward in the emulsion and 
stops near the surface (see fig. 11). There is a point at which the 
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Mass-determination by means of crossing tracks. 
Above: vertical section showing the tracks. 
Below: density versus range for the same tracks. 


tracks are at the same distance from the glass, and in that point the 
degree of development is the same if it is a function only of the distance 
from the glass. Measurements on both tracks have shown that the 
crossing point is 2575 from the ends of both tracks. The density 
of the proton track at that point is 49-6 units and in the meson. track it 
is 47:7 units. Thus the density difference at the point of intersection 
is 1-9 units. The mass-difference to which this density difference is 
equivalent can be found from fig. 9, where the density curves for 
different masses can be constructed from the proton curve. A particle 
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with mass 1500m, has a density difference of 1-9 units at a range of 
2575 p. 

The error in the estimates made in this way are large, for several 
reasons. The greatest error is a consequence of statistical density- 
variations along the track. Furthermore, the curve cannot be drawn 
with as great certainty as it can after correction because the curve form 
is not simple. In a rough way the error in density has been found to be 
about 0-8 unit for each curve corresponding with a probable error in 
mass of 200 m,. 

Thus the result of the present «,-meson mass determination is 
(1500-200) m 

§7. THe Mass OF THE x,-MESON 

The mass of the x,-meson has thus been calculated in two ways from 
the measurements of the projected area of the track. In the first 
determination the mass has been calculated after corrections of the 
density. Taking account of the errors of the method the mass was 
found to be (1521--75)m,. The second method, wh2re crossing tracks 
were used, gave the result (1500+ 200) m 

The mean value of the two determinations has been calculated to. 
(1518-+-70) m 
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SUMMARY 


The paper describes calculations of the exchange energy of an assembly 
of nearly free electrons of low density (§ 2). Two differing one-electron 
wave functions are used, given by relations (2.7) and (2.22), both 
concerning a nearly empty Brillouin zone. It is found in both cases 
that the total exchange energy is much less sensitive to changes of the 
effective mass of the electrons than is the Fermi energy. Previous 
related work is referred to (§1) and the bearing of the results on some 
physical properties of metals is indicated (§ 3). 


§1. INTRODUCTION 

THE calculation of the exchange energy of an assembly of perfectly free 
electrons constituted one of the early advances of the modern electron 
theory of metals (see, for example, Sommerfeld and Bethe 1933). Use 
of one-electron wave functions of the form 

pete Vaal ex Ui tKpee ts), uemeroe ays, 3) +) (1.1) 
where k,,, r; represent position vectors in reciprocal and real space 
respectively and V is the volume of the metal, gives for the exchange 
energy for two electrons with parallel spins and wave vectors k,,, k, 
(Bloch 1929) 


Fn =2 | f ahr adn* Pad (Fan 2) dF Po) 


(1.2) 
as Acre? il 
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The total exchange energy of the metal in this approximation is given by 
3N\ 138 
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if the metal is non-ferromagnetic, N being the number of electrons in 
volume V. 

The relation (1.2) may be used in discussing such properties of the free 
electron gas as the temperature variation of the electronic specific heat 
(Wohlfarth 1950) and of the spontaneous magnetization for a ferro- 
magnetic metal (Lidiard 1951). In such discussions it appears that, 
because of the rapid variation of Jn» as |k,,—k, |>9, results are 
obtained which differ widely from observation. It is thought at present 
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that these deviations may be partly removed by taking some account 
of the effects of correlation forces between electrons with antiparallel spins. 
On the other hand, the relation for J,,,, based on one-electron functions 
of the tight binding type is of an entirely different form from (1.2) 
(Wohlfarth 1953) and more reasonable thermal and magnetic results 
may be derived. The present calculation, based on nearly free electron 
functions, makes an attempt to interpolate between the behaviour of 
the perfectly free electron gas and of an assembly of tightly bound 
electrons. 

Again, relation (1.3) may be used immediately in calculations of the 
cohesive energy of metals in the free electron approximation (Wigner 
and Seitz 1934) and in deriving the condition under which the metal is 
ferromagnetic (Bloch 1929). The present calculation of J, using nearly 
free electron functions, gives some information on such questions as the 
effect on the total exchange energy of an effective electronic mass different 
from the actual mass, of the presence of Brillouin zone boundaries and 
of the shape of the Fermi surface (cf. § 3). 

The problem of calculating the exchange energy for anything but the 
free electron gas is one of enormous difficulty and in this paper 
calculations are carried out for two limiting cases only (cf. § 2). Other 
attempts along somewhat similar lines have been made previously. In a 
calculation of the cohesive energy of lithium Seitz (1935, 1940) expresses 
the wave function in the form 

byl) =(a-yk,,.6,) exp (tk, tj). «5 eee eee 
where « and y, strictly functions of position inside each atomic cell, 
were taken to be constants. The calculated exchange energy does not 
differ appreciably from the free electron value (1.3), the correction being 
such as to increase J numerically. Herring (1951) has, however, recently 
found that the sign of this correction should be opposite to that obtained 
by Seitz, whose calculation is stated to involve a computational error. 
Herring and Hill (1940) carried out a rough calculation of this correction 
in their work on beryllium. This calculation applies to the limit of low 
electron density and thus bears some formal resemblance to the present 
calculation, although the wave functions used in the two cases are 
different. With q the number of electrons per atom (¢<1) Herring and 
Hill’s formula is of the form 


J=d-, {1+q??+0(¢4)}, 5 a+ ah bee oe Ge) 
where J;.. is given by (1.3) and ¢, a function of the radius of the atomic 
sphere, was obtained by numerical calculation for beryllium. Here ¢ is 
negative and small, in qualitative agreement with the present results. 


§ 2. CALCULATION OF THE EXCHANGE ENERGY 


A general relation for J,,,,, based on the nearly free electron approximation 
had been given in an earlier paper (Wohlfarth 1949, relation (4.8)) and it 
was stated that, owing to the complexity of this relation and of the shape 
of the Fermi surface, the summation process to give J as in (1.3) would 


Exchange Energy of Nearly Free Electrons in Metals 283 


be prohibitively difficult. It is, therefore, proposed to delimit the problem 
considerably. Two cases will be considered : 

(i) Interaction between electrons overlapping into the second Brillowin zone. 
It is assumed that the first zone is full and that overlap occurs in the 
neighbourhood of k=(-t7/a, 0,0), where a is the lattice parameter. 
With g, the number of electrons per atom in the second zone, small 
compared with unity, this situation is best visualized by reference to fig. 1, 


Fig. 1 


RIA 


(6) 


- (c) kx +a Chin Nae 


(a) Electron distribution in k space, k,=0; extended zone scheme. (b) Reduced 
zone scheme. (c) Variation of energy with k,, k,=—k,=0; AH, energy gap ; 
Ey=(h?/2m)(z/a)?. (d) Density of states curve. All figures are schematic. 


Figure 1 (a) shows the overlap, the two parts of the Fermi surface 
being such as to combine to form an oblate ellipsoid of revolution, since 
the states k=(-:7/a, 0, 0) are effectively identical. This ellipsoid may 
conveniently be regarded as lying at the centre of the zone shown in. 1 (5), 
which is that of 1 (a) displaced by +-x/a. 1 (c) shows the variation of 
the energy with k, and 1 (d) the density of states curve, all the figures 
being schematic. The calculation to be described is concerned with the 
few electrons distributed in the ellipsoid in 1 (b) or, alternatively, in 
the parabolic part of the M(H) curve in 1 (d) above the energy gap. 


UZ 


284 E. P. Wohlfarth on the 


As is well known (Mott and Jones 1936, p. 65) the energy is here given by 
h2 
B= 5 (abpthp+ke), . . . + . + (21) 


where «=1-++-4H,/4E, Ey=(h?/2m)(n/a)?. In (2.1) k, is measured relative 
to the centre of the zone in 1 (6) and relative to k,=-t7/a in 1 (a). 
The ellipsoid has semi-axes a, b, b (this a should not be confused with the 
lattice parameter), where 


G=a eke b=a Nek, ho (S07 Ni] V)EP se ee) 
Strictly the effective electronic mass m* has to be expressed in the tensor 
form mesh. ty 
m/m*= f 1 i gh tes es ea) 
OF Ose 


The Fermi energy, obtained by summing # in (2.1) over the occupied 
states, may, however, be written in the form 


hi2 
By=$N sah, 2 2. 2 + 6 (24) 


which is identical with the expression for quasi-free electrons with an 
effective mass equal to eT Re Re | ee ee CL 


Since « may be large, say ~ 10, the Fermi energy of the overlapping 
electrons may be 2-3 times larger than that of the same number of 
perfectly free electrons. The problem is to calculate the corresponding 
change in the exchange energy. 

The one-electron wave functions may be written as 


Pm i)= {Aom+A rm XP (—27iw,/a)} exp (tk, .6;). —«. (2.6) 
A simple transformation gives the alternative expression 
Pal ti) ={2/V UAB han?) 2 {sin (7727/4)-+4 Blam 08 (772,/a)} exp (iKyy « F4) 
ootaty mae 
correct to the second power of k,,,,, where k,, is measured relative to the 
bottom of the second zone as in 1 (a) and 
P=a(a—l)ome 90 a. Gay Se eee 
For k,,,=0, (2.7) represents a standing wave in the 2 direction (cf. Mott 
' and Jones 1936, p. 62). 
The calculation of J,,,, proceeds by standard methods. It is noted 
that J,,,, may be written in the form 


Jinn=eJ U (re )bn* (Pe )n (te) dr(ry), 
where V?U(r,)= — 4rreys,(Pa)bn*(ry), zat 
and it is found that 
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where p=27/a, K=(k,,—k,,)=(K,, Kz, K;). For K <m/a (2.10) may be 
expanded up to K? and J,,,,, becomes 


Ane? [ 1 pl at i Ke 2 
I m= Ee = at Sc Ona {ham + kan*)( B+ 1) 
—heaphin(2+208-+2°6)} |... + 1) 
The problem is now to sum J,,,, over the ellipsoid defined by (2.2). 
Only the first two terms in (2.11) present any difficulty and a brief outline 


is given of the method adopted of integrating the first term. It is required 
to evaluate 


dr (R,, Rg) XA YZ? 
fis) ae eee ie ear eee eee 281 
1 [RE ae te ow 
By use of Poisson’s equation it may be shown that 
1 dr(r) 5 
I=, { = JJ exe fa(R, —R,) . ¢} d7(R,, Ro). 


Consider the integral 
= iis exp (2R.r) d7(R). 
Put X=au, Y=bv, Z=bw, then d7(R)=ab? dr(e), where p=(u, », wv). 
The domain of integration is now the unit sphere M’=w-+v’?+w=1. 
Also R.r=p.r’, where r’=(ax, by, bz). Hence 
T=al? | exp (ie. r’) dr(p)=47ab? ( 
ee 


ae T=4ra°b? | eu (my 


sin r’—r’ cos 7” 
St ESS ae ie 


72 
or, transforming r=(a,y,2z) to r’=(ax, by, bz)= (airy ne), 
gl — y/212/2\—1/2 /gin r’—r’ cos r’\? 
= 2 , alae’ SS See 
I=4rab {ane (S Sle be ) ( 73 ) p 
The integrations are elementary and finally give 


1+ (19?) 2 
(yy? 


[=2n2a2b?(1—7?2)-1/? In (2.13) 


where n=a/b. Hence, using (2.2), 
T=47?k,Afy(«), — f(«)=aN8(a—1)-1? cosh-*a?, (2.14) 


Since f,(1)=1 the corresponding free electron expression is obtained 
correctly for «=1. The limiting expressions are 


a1, fi(a)=I—F(e—1)+-0((a—1)9, sits 
Garton aft) hone In do, 
Tt follows from (2.14), (2.15) that f,(«) is a very slowly varying function. 


Even for «!%=m/m* as large as 3, f;(%)= 0-79, so that a change of 
300% in the Fermi energy produces a change of only about 20% in the 


ial 
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leading term in the exchange energy, compared with the corresponding 
free electron values. "1. 

The integration of the second term in (2.11) proceeds by similar means 
and only the final expression for J need be given. It is convenient to 
express J as a power series in q?/*, where q is the number of electrons 
per atom. This is achieved by noting that for the simple cubic lattice 
the lattice parameter a is given by a=(qV/N)"%, where N is the total 
number of electrons, slightly modified relations applying to other lattice 
types. Then, from (2.8), B=(qV/N)"3(«a—1)/27. The final expression 
then becomes, with J;,, given by (1.3), 


Jmdead Het (3)"" geatthtort (2) ae atrher+ou}. 2.16 
Here f,(«) has already been defined in (2.14) and (2.15). Also 
fo(a)=1—2(a—1)(a"Ff;(a)—1), 
a>1, f,(«)=1—2(«—1)?+O0((«—1)4), = See) 
a—+>oo, f,(x)=1+2«—« In 4a, 
so that f,(«) varies much more rapidly than f,(«), soon becoming negative. 
eoay j,(a)=at(2a 1), >. 


again a fairly rapidly varying function. Since, however, the numerical 
coefficients in (2.16) become rapidly very small for increasing powers 
of qg?/> the total exchange energy differs only slightly from the free 
electron value, since the effective mass has little influence on the leading 
term and the convergence of the series is rapid, especially if g<1 as is 
assumed in the calculation. 

(ii) Interaction between electrons near the centre of the first Brillowin zone. 
Here, if g<1, the electrons are, to a first approximation, described 
correctly by free electron functions (1.1). To a second approximation, 
if only the perturbation arising from the first zone boundary is important, 
the energy in the one-dimensional case becomes (Mott and Jones 1936, 
p. 61) ys 2m __(AE/2? 

2m RE (k—2rr]a)¥—K2* 
Here, since ¢<1 implies k<z/a, expansion gives 
HE’ =E+8E y= (h?/2m)(k?(1—8")—782k/a+. . .), 


where 6= 4EH/4E,, cf. fig. 1. The energy has to be measured relative to 
the lowest state, k=0, i.e. relative to H’=0. The effective mass is given by 
mm d?K’ 
m* fh? dké 


(2.19) 


=1—8—8*0(q¥8), 2... (2.20) 


In the three-dimensional case, unless 5=0, the Fermi surface is no longer 
spherical. For the purpose of the calculation of exchange energy the 
deviations may, however, be neglected, since the above calculation 
shows how insensitive is the exchange energy to even quite large 
deviations of the Fermi surface from a sphere. With this approximation 
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the effective mass near the centre of the zone is isotropic, in contrast to 
the above case, and it is given by (2.20). The wave function is 


Pl ¥:)=(Aom+2A 1m (COS PE; COS py; Cos pz,)) exp (ik,,.¥;), (2.21) 


where p=2n/a. It is required to calculate the exchange energy correct 
tog. To this approximation (2.21) becomes 


Pn t= VL dyn?) {1 +-Yn (COs px; +-cos py; +cos pz;)} exp (ik, r,), 
where Var (9 lTT) (erat Maen ena es ne 5h ae id (2:22) 
The relation for J,,, is obtained as before and after expansion to terms 
In (46/7), i.e. in 82g2/8, it is found that 

47re” 1 
mn = WT Tk, —k, B80 m + Yn? —Ym¥n)} + + (2.23) 


The total exchange energy J is obtained by integrating (2.23) over the 
occupied states, the Fermi surface being here taken to be spherical as 
discussed above. The integrations are elementary and the final 


expression is 7 is 
7 
T=Jee41— os (5) argr?+-sro(qt)} berate Bch ~(2e24) 


where J,,, is given by (1.3). Comparison of (2.20) and (2.24) shows 
that the ratio of the relative changes of the exchange and Fermi energies 
on changing the effective mass is given by 

I —JSeo.)/F £0. 

cp ute == 1329724-O(g). 9. «Ce (2.28) 
Since ¢<1 it appears again that the exchange energy depends much less 
sensitively on changes of the effective mass than does the Fermi energy. 
As already noted relation (2.24) bears some formal resemblance to Herring 
and Hill’s relation (1940) given by (1.5). In the present calculation the 
correction term is closely connected with the magnitude of the energy 
gap at the first zone boundary. 


§ 3. Discussion 


The present calculation is based on an approximate treatment of the 
metallic state whose validity in describing any particular metal is doubtful 
since the potential energy of the electrons is not everywhere small 
compared with their kinetic energy. Further, the results obtained in § 2 
are strictly applicable only in the limit of low electron density. The two 
calculations described above are, however, based on two entirely different 
types of one-electron wave functions, (2.7) representing a slightly 
perturbed standing wave and (2.22) a slightly perturbed plane wave, 
both for electrons in a nearly empty zone, but with differing effective 
electronic masses. The fact that qualitatively similar results are obtained 
in both cases indicates that conclusions drawn from these calculations 
may have a wider validity than is justified by the premises underlying 
the calculations. 

(i) Relations (2.11) and (2.23) show that the leading term in J,,,, is in 
both cases equal to the value (1.2) for perfectly free electrons. Hence 
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the present calculation still predicts an anomalous temperature variation 
of electronic specific heat and spontaneous magnetization (cf. §1). The 
calculation has, in fact, failed to provide an intermediate case between 
free electron and tightly bound electron behaviour. 

(ii) Relations (2.16) and (2.24) show that the calculated exchange energy 
varies much less rapidly on changing the effective mass of the electrons 
than does the Fermi energy. This conclusion may provide some 
explanation of the fact that in calculations of the cohesive energy of 
many metals good agreement with experiment is usually attained by 
assuming the exchange energy to have the free electron value Jee, 
although great care is often necessary to correct the Fermi energy if 
m/m* differs appreciably from 1 (see, for example, the discussion of 
beryllium by Raimes 1952). In particular, relations (2.14)-(2.16) show 
how insensitive is the exchange energy to deviations of the Fermi surface 
from spherical shape. This result may have some bearing on the problem 
of the contribution of exchange forces to the elastic constants of metals 
(Leigh 1951). 

(iii) The condition for ferromagnetism of a perfectly free electron gas 
is that the ratio of exchange energy J to Fermi energy Ey given by 
(1.3) and (2.4) («=1) be larger than 2-5 (Bloch 1929). The second 
calculation in §2 shows that for nearly free electrons, to the second 
approximation, the Fermi energy is reduced by a factor (1—82) (ef. 
eqn. (2.20)), while to the same approximation (2.24) shows that the 
exchange energy is practically unchanged. Decrease of the Fermi 
energy, i.e. increase of the tightness of binding and narrowing of the 
energy band, thus increases the ratio J/E, and therefore the tendency 
towards ferromagnetism, a well known result. 


It gives me great pleasure to acknowledge the help I have had from 
Dr. A. N. Gordon in carrying out the evaluation of some of the integrals 
used in the calculations. 
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ABSTRACT 


The experiments were made (i) with a new apparatus in which a bath 
of liquid helium can be cooled to 0-2° K and (ii) with a two-stage demag- 
netization apparatus in which temperatures of 0-006° K can be reached. 

The results showed that (i) for specimens in different states of physical 
strain down to 0-2° x and (ii) for a specimen down to 0:006° x the resis- 
tance is linear with log 1/7. This is in contrast to the relationship to be 
expected from Slater’s theory, i.e. R proportional to exp 6/7’. 


THE anomalous increase in the resistivity of certain ‘pure’ metals 
with falling temperature which occurs at temperatures below about 
6° is now well established experimentally, especially in the case of 
gold (for references see Mendoza and Thomas 1951). Slater (1951) 
has suggested that this behaviour may be due to a rearrangement of the 
gold lattice such that an energy gap occurs in the occupied Brillouin zone 
just at the top of the distribution of the conduction electrons. The metal 
would then behave like an intrinsic semi-conductor and a rapid rise in 
the resistivity with falling temperature would be expected. Some 
limitations of Slater’s theory have been discussed by Katz (1952). 

This paper reports the results of measurements of the electrical 
resistivity of gold made 

(I) in a new apparatus in which a bath of liquid helium can be cooled 
to 0-2° x, and 

(II) in the two-stage demagnetization apparatus described by Darby 
et al. (1951), in which measurements have been made down to 0-006° kK. 


Part I. Experiments with the Liquid Helium Demagnetization Cryostat 
(A. J. Crorr and E. A. FavtKNER) 


The main difficulty in measuring electrical resistivities at temperatures 
below 1° K is to get satisfactory thermal contact between the specimen 
and the paramagnetic salt which is used for achieving the cooling and 
measuring the temperature. The usual methods, such as compressing 
the salt round the wire, are liable to strain it. A better course is to sub- 
merge both in a bath of liquid helium, for this still has a useful thermal 
conductivity at 0-1°K. This is made possible in this apparatus by the 
use of the method of Hudson, Hunt and Kurti (1949) for reducing the 
heat influx due to recondensation of helium transported by the film 


* Communicated by Professor F. E. Simon, F.R.S8. 
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to the warmer regions. A container with the metal specimen, the 
paramagnetic salt and liquid helium inside it is suspended on a capillary 
tube (0-2 mmi.d. and 100 mm long) which limits the amount of 
vapour able to recondense and yet provides a way for getting the helium 
in and out. 

The complete apparatus is shown in fig. 1. The inner bath is suspended 
inside an evacuable space which is at the bottom of a liquid helium dewar 
attached to the pumping line by a copper/glass seal. This outer bath 
holds sufficient liquid helium when it is full at 4° K to leave enough for 
an experiment several hours long after it has been pumped down to 1-1° K 
and has cooled down the inner bath as well. (The use of glass has many 
advantages, but even though Monax glass is used, the dewar has to be 
repumped after a few experiments.) 

Of the paramagnetic salts for which the absolute temperature scale 
has been determined, the most suitable is gadolinium sulphate ; it gives 
a relatively large entropy change/unit volume on magnetization and has 
a very large specific heat between 0-2 and 0-7°K. Using a 100kw 
water-cooled magnet which gives 11 ke, and demagnetizing from 1-1°k, 
we reach a temperature just below 0-2°K. The ballistic method of 
measuring the magnetic moment of the salt is used. To deduce the 
absolute temperatures from these measurements we have applied the 
corrections given by*van Dijk (1946). As cylindrical samples were 
used in both his work and ours, it is difficult to allow for the differences 
between the sample and secondary coil shapes in the two cases. For 
this reason the values of the absolute temperature are probably not very 
accurate at the lower temperatures. 

The experiments so far done in the cryostat have been measurements 
of electrical resistance for which we have used the galvanometer amplifier 
circuit developed by other workers in this laboratory (MacDonald 1947). 
The leads are brought intothe inner bath through holes sealed with Araldite, 
a thermo-setting plastic ; a tungsten pinch sealed to a small copper/glass 
seal is used to bring them through the top of the vacuum space. The 
cryostat has proved easy to use ; the normal rate of warming up is about 
0-15°K/hr at 0-3°x. If one wanted to go to lower temperatures than 
0-2°k by using, for example, iron ammonium alum, the heat input would 
have to be cut down; this could be done by using an extra liquid helium 
bath at the top of the capillary tube. 

The gold wire was 0-1 mm in diameter and had been drawn by Messrs. 
Johnson, Matthey (their no. 3459), using first a tungsten carbide die 
and then a diamond one. 

Their analysis is 
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A length of 1-25 m was loosely wound on a former; the ends were 
soft-soldered to connecting lugs. Experiment I was made with the 
wire taken straight from the reel. Experiment II was done on the same 
length after it had been annealed in a vacuum at red heat for three 
hours. Before Experiment III it was annealed ‘again, but at a higher 
temperature. In each case the measuring current was 5 mA. 


Results 


The results of these experiments are shown in fig. 2; above 1°K 
each point represents the mean of several resistance measurements made 
at a constant temperature ; below 1° Keach point represents one reading 
of temperature and resistance. 


Fig. 2 
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The straight lines conform to 
R/R,=const.—a log T, 


where R, is the value of the resistance at Or ker 
The magnitude of the increase can be expressed by a coefficient « 
which is defined by 
103 dR 
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The following table shows the residual resistance ratio and the value 
of « found in the three experiments :-— 


Experiment I: R/R,= 97xX10-4 at 4°K a=11 
Experiment IT: = 48 a=58 
Experiment IIT: =103 a=23 


Comment 


(1) The experimental points lie quite close to the —a log T line between 
0-25° K and 2°x. The deviations from the straight line below this tem- 
perature which appear in II and III are probably due to the uncertainty 
in the corrections to the absolute temperature scale which has already 
been mentioned. It does not show up in I because of the lesser slope 
and greater scatter. (Similar unexpected deviations have been found in 
other work not yet published.) 

(2) The first annealing reduced R/R), as one would expect, and 
increased «. The wire seems to have become contaminated during 
the second annealing because R/R, became larger than it was originally ; 
it is interesting to note that the value of « is between the values found 
in the first two experiments. 

(3) The purity of the wire in Experiment II as indicated by the value 
of R/R, is between that of Mendoza and Thomas (1951) (R/Ry=83 x 10~*) 
and that of de Haas, Casimir and van den Berg (1938) (R/Ry=22 x 10~*). 


Part II. Experiments with the Two-Stage Demagnetization Apparatus 
(J. Harron and E. F. W. S—ymour) 


In the measurement of electrical resistance at temperatures too low 
for the specimen to be immersed in liquid helium, special attention 
must be paid to the thermal contact between the metal and the para- 
magnetic salt used to produce the low temperature. In the experiment 
described here, the gold wire was embedded in the paramagnetic salt. 
It was decided that a potential difference of about six microvolts across. 
the specimen could be measured with sufficient accuracy (about 0-3%). 
With the expected resistance of the gold at the low temperatures, this 
potential difference was calculated to produce a power dissipation of 
only a few ergs per minute. Further, since the measuring current would 
be on for only short periods, it was considered that the heating due to it 
would not be too high. Unfortunately it is not possible to make any 
quantitative estimate of the tolerable power dissipation in the wire 
since neither the thermal conductivity of the paramagnetic salt nor the 
efficacy of the contact between it and the wire is known at the temperature 
involved. 

The gold wire used was from the same reel as that used in the experi- 
ments described in Part I of this paper. 
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The first stage of the demagnetization apparatus shown in fig. 3 
was a cylinder pressed from powdered iron alum. The second stage 
consisted of crystals of a mixed salt of potassium aluminium alum and 
5% potassium chrome alum. The crystals were powdered and compressed 
around the gold wire in such a way that there was no contact between 
different parts of the gold wire or between the gold and the various 
copper wires which were also embedded in the salt. 

Fig. 4 
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Current and potential leads passed through the first stage and then to- 
the second stage via the lead wires which provided the heat switch between 
the stages. Each of the four leads was soldered to copper wires embedded 
in the second stage to increase the thermal contact. The ends of the. 
gold wire were soldered directly to the current leads, but the potential 
connections were made by small gold strips clamped on to the gold wire 
a ner ae “ the potential leads. This precaution was adopted 

o avoid the possibility of disturbing eff i i 
epee nah ni ease g elfects which alloying of the gold 

The demagnetization procedure and the temperature measurement of 

the paramagnetic salts were as described in Darby et al. 1951. 
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The resistance measurements were made in a conventional manner. 
A current of 0-88 ma was passed through the specimen. The voltage 
across a standard 1 ohm resistance in the current circuit, measured 
with a potentiometer, enabled accurate measurements of the current 
to be made. The voltage across the gold specimen due to this current 
was measured by means of a feedback galvanometer amplifier of similar 
design to the one mentioned in Part I. 


Results 


The resistance of the gold wire at room temperature was approximately 
0-75 ohm and at liquid helium temperatures 0-0075 ohm. 

The lowest temperature at which resistance measurements were made 
was approximately 6-5x10-?degk. As the specimen warmed up, tem- 
perature and resistance measurements were made each minute to begin 
with, and later on at longer intervals. The results are shown in fig. 4. 
The value of « as defined above is approximately 10, which is in good 
agreement with the results of Experiment I of Part I of this paper. 

The total heat inflow to the second stage from all causes was 25 erg/min. 
The electrical power dissipation in the gold wire during measurements 
was approximately 3-5 erg/min, but the measuring current was on for 
only about 10 seconds per measurement. It seems reasonable to suppose, 
therefore, that the temperature of the gold wire would not be appreciably 
different from that of the surrounding paramagnetic salt. This conclusion 
is supported by the fact that the measured value of the resistance at a 
given temperature was independent of whether the measuring current. 
was on for only a few seconds or for periods of about thirty seconds. 


Conclusions 


The results of these two series of measurements show that the electrical 
resistivity of gold at low temperatures increases with falling temperature 
approximately according to 


Rx log (const./7’) 


down to about 0-006° x. If the increase were due to behaviour similar 
to that of a semi-conductor, a much more rapid rise would be expected, 


that is Recexp6/T 


with 6 of the order of magnitude of 1°x. The resistance at the lowest 
temperature would then be orders of magnitude greater than its value 
at the minimum, whereas the observed increase is only about 3%. 
The results for wires subjected to different annealing procedures show 
no clear correlation between the magnitude of the resistance increase 
and the physical state of the metal. Unfortunately it is not possible to 
decide whether the effect is a property of perfectly pure gold or is deter- 
mined by the presence of impurities (see, however, MacDonald 1952). 
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SUMMARY 


A way of simplifying the microscope measurements and the resulting 
calculations, in the photographic emulsion method of studying neutron 
spectra, is described. The energy distribution of neutrons from the 
reaction /B(d, n)'!C has been measured, to show that the saving in time 
is accompanied by no appreciable loss of accuracy. 

The method has been used in a study of the angular distribution of 
neutrons from the reaction !4B(d,n)"C at 8-1mev. The results are 
consistent with the theory of the stripping process and confirm that the 
4-47 mev level in C, like the ground state, has even parity. 


§1. INTRODUCTION 


THE use of photographic plates for measuring neutron energies has been 
described by many authors (e.g. Powell 1940, 1942; Richards 1941). 
Measurements are made on the tracks of protons which have moved 
forwards, after neutron collision, in directions making small angles 
with that of the neutron flux. The energy of each neutron, H#,, is then 
deduced from the proton energy, H,; the angle between the directions 
of motion of the neutron and the proton, ¢, is obtained from the relation 
H,=E, sec* ¢. 

The usual procedure has been to allow the neutrons to enter the plate 
parallel to the surface of the emulsion, and on each track to measure the 
following : the length, a; the angle between the direction of the neutrons 
and the horizontal projection of the initial direction of the proton, 0; 
and the vertical component of the first portion of the track, h. h gives 
the angle of dip 8, a correction being made for the shrinkage of the 
emulsion in the vertical direction during processing ; the range of the 
proton, R, is equaltoasec5. The angle 4 is determined from the equation 

cos ¢ = cos @ cos 6. 
The range-energy relation gives H, when # is known, and £,, is calculated 
from H, and ¢. 

This method has proved satisfactory in many studies of neutron energy 
distributions. However, the microscope work is slow, and each track 
requires several calculations ; these facts have limited the number of 
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tracks which it has been practicable to measure in any particular plate ; 
consequently, the statistical accuracy of the value of the neutron flux, 
obtained from the number of tracks found in a known volume of emulsion, 
has not usually been great. 

A new technique has therefore been developed to reduce the time used 
in measuring and calculating. It is hoped that the measurement of 
neutron fluxes, and the comparison of fluxes to obtain angular distri- 
butions of neutrons from nuclear reactions, will be made less arduous by 
this technique. 


§2. Tur IMPROVED METHOD 


Suppose that the beginning of a track is in focus at a fixed point in the 
field of view, and that the objective lens or the plate is moved through 
a small vertical distance F. The part of the track then in focus is defined. 
by its intersection, as it was formed in the unprocessed emulsion, with a 
plane SF above or below the beginning of the track ; S is the factor by 
which the thickness of the emulsion is reduced during processing, and 
is about 2-1 for Ilford Nuclear Research emulsions. A track with a 
given ¢ will appear somewhere on the hyperbola formed by the inter- 
section with this plane of a cone with semi-angle ¢, apex at the fixed 
point where the beginning of the track was set, and axis parallel to the 
neutron direction. A special graticule carrying a series of these hyper- 
bolae could therefore be made and used in a microscope eyepiece for the 
determination of ¢ in a single operation. 

It is possible, however, to reduce the labour of calculation still further : 
for protons of energy about 10 Mev, the range in Ilford C2 emulsion 
varies as H'°, Since H,=E, sec*¢ the range which would have been 
observed for a proton with the whole of the energy £,, is 


Ro =R(sec Pees 


f being the range observed for a track at an angle ¢. FR is not usually 
measured directly, since it is easier to measure L, the component of the 
length in the direction of the neutrons ; R is then equal to L sec ¢, and 


Ry =L(sec ¢)**°, 


Hence it is best to use a graticule with hyperbolae plotted at equal inter- 
vals of (sec ¢)*4°, since with such a graticule we need measure only LZ 
and (sec ¢)*4°, The product Ry may then be plotted, and the distribution 
of ranges re-plotted as a distribution in energy if this is required. 

It has been assumed throughout that the tracks are straight: in 
fact there is some scattering, and it is difficult to estimate the inaccuracy 
this introduces. The results show, however, that the errors are no 
greater than in the older method. 

The appearance of the graticule used in Bristol is shown in Dgeedes 
it was constructed by two-stage photographic reduction of an ordinary 
drawing, onto a plate of very high resolution. The outer hyperbola 
has (sec ¢)*4=1-30; this corresponds to ¢=19-5°. The overall size of 
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the graticule, and its position in the field of view, are chosen so that no 
appreciable error can be caused by the curvature of the microscope 
field. It must be confirmed that there is no excessive curvature of the 
field of any microscope which is to be used for this work, but no difficulty 
of this sort has yet been encountered. 

The known vertical movement is provided by a special attachment 
on the fine-focussing adjustment. The ordinary control knob (I) is 
connected to a second knob (II) on the same axis, in such a way that 
there is free relative movement through a small angle. Because of the 
presence of backlash in the microscope mechanism, the setting of the 
start of the track on the origin is done by moving II in one direction 
and pushing I with it ; the known further movement is made by moving I 


Fig. 1 


Special eraticule for measuring the factor (sec ¢)** for tracks of recoil protons. 
z =: The numbers are % corrections, i.e. L00{(sec ¢)*s-1}. 


in the same direction through the whole of the interval of free relative 
motion. This interval defines F’, and its exact magnitude is adjusted to 
suit the particular objective and graticule in use. Since the scale of 
each hyperbola is proportional to the magnification, to the vertical 
distance F, and to the shrinkage factor S, a given graticule can be used 
with any objective provided that F is suitable adjusted. For the 
graticule shown in fig. 1, used with a X45 objective, F is 2-9y for plates 
: ne 
ne aa L, is measured on the micrometer head of the stage 
adjustment, modified by the addition of a ratchet device ; the zero is 
set when the beginning of the track is moved up to the vertical line in 
the graticule, and the length is read directly when the end of the track 
is on the same line. 
x2 
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The accurate measurement of dip is always difficult, and in this 
procedure it is the measurement of dip implicit in the estimation of 
the factor (see ¢)!4 which leads to the chief inaccuracy. However, 
it is found that the mean error in this factor is about 1%, and this is 
not excessive. The assumption of a simple power law for the range— 
energy relation can lead to error if widely differing neutron energies are 
being measured: calculation shows that the angle 19-5° for which the 
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Neutron energy distributions from the reaction 1°B(d,n)14C. Curve (a) shows 
the new results, curve (b) those obtained in 1949. Each ordinate is 
proportional to the number of neutrons in a 0-2 Mev interval. The factor 
by which the number of proton tracks observed was multiplied to give 
the neutron intensity, in the arbitrary units used, varied from 1-0 at 
7 Mev to 0-6 at 4:5 Mev. 


outside hyperbola is drawn and which corresponds to a factor 1:30 at 
10 Mev, gives factors 1:305 at 16 Mev, 1-290 at 6Mev and 1-282 at 
3 Mev; the difference is correspondingly less for the other hyperbolae 
This graticule can therefore be used without appreciable systematic ays 
between about 6 Mev and 20 Mev; at lower, or much higher, energies 
it would be necessary to use slightly different graticules if the greatest 
possible accuracy were to be obtained. 


ee 
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The method was tested with plates exposed to the neutrons emitted 
at 90° from a 1°B target bombarded with 930 kev deuterons. Similar 
plates had been made the basis of an earlier investigation (Gibson 1949). 
In fig. 2 the new results are compared with those obtained in 1949; the 
resolution and the accuracy of the new method are at least equal to those 
of the older and more laborious technique. 
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P d . : 7 11Bi(d n)??C. (a) 120, 
C-m-s angular distributions for neutrons from the reaction ,n)™ 
rasa g in its ground state, (b) ?C formed in 4-47 Mev excited pies: 
Full lines are the curves given by Butler’s theory, with ly= ie an. 
r—4-4 10-13 em; broken lines are given by the Born-approximation 
theory with 1,=1 and r=—6-4 x 10-" cm. 


The number of tracks in the group attributed to the formation of 4C 
in an excited state at 2-02 Mev was not great in the original results. 
The new figures confirm the existence of this level; since the groups 
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appear at the same energies as before, the energy of the excited state 
may still be taken as 2-02 Mev. 


§3. 11B(d, n)#C 

The method, having been tested and found satisfactory, was used in a 
study of plates exposed to neutrons from the reaction 1B(d,n)¥C ; 
in this case the deuterons were accelerated to an energy of 8-1 Mev in the 
University of Liverpool cyclotron. The experimental arrangements were 
the same as those described by Middleton, El-Bedewi and Tai (19583), 
who have discussed fully the importance of the angular distributions of 
neutrons from this and similar reactions. 

Two groups of neutrons were clearly resolved, and the energies indicated 
that they were emitted when ”C was formed in its ground state and in 
the excited state at 4:47 mev (Gibson 1949). Measurements were made 
on these two groups at a number of angles, and the differential cross- 
sections were calculated. The calculations included allowances for the 
loss of proton tracks through the surfaces of the emulsion and for the 
variation of n, p scattering cross-section with energy, and the conversion 
to the centre-of-mass system. The angular distributions obtained are 
shown in fig. 3. 

The results for neutrons from formation of ?C in its excited state fit 
well with theory, if it is assumed that an angular momentum /,=1 is 
transferred to the nucleus by the proton which it absorbs. Two theore- 
tical curves are shown in fig. 3: the one with which agreement is closer 
is based on Butler’s theory of the stripping process (Butler 1951), with a 
radius of 4:4x10-%cm (the Gamow radius); the second curve is 
obtained from the Born-approximation theory (Bhatia et al. 1952) with 
a radius of 6:4x10-%cm?. Both theories, however, show quite 
definitely that the angular momentum transferred is 1,=1. If the 
ground state of 14B is assumed to have spin 3/2 and odd parity, the 
4-47 Mev level in #C must then have even parity; the results give no 
information about the spin, except that it must be 0, 1, 2 or 3. Ajzenberg 
and Lauritsen (1952) suggest spin 2 and even parity. 

Theory suggests that formation of ”C in its ground state, by a 
stripping process, is less likely than formation in the first excited state 
(W. M. Fairbairn, private communication). This is confirmed by the 
results shown in fig. 3: the differential cross-sections are lower, and 
it is probable that competing reactions partly mask the characteristic 
form which the stripping process gives to an angular distribution. The 
behaviour below 30°, however, is consistent with the value L,=1; this 
value is expected, since the ground state of 12C has even parity and spin 0. 
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SuMMARY 


Methods are described whereby stars produced in nuclear emulsions 
by protons, in the energy range 110 to 146 Mev may be separated into 
disintegrations produced in the light (carbon, nitrogen and oxygen), and 
heavy (silver and bromine) components of the emulsion. The cross 
sections obtained for nuclear disintegrations are compared with those 
predicted by the ‘ transparent nucleus ’ theory of Fernbach, Serber and 
Taylor (1949). 


§1. INTRODUCTION 
Tur three components which determine the total cross section for the 
nuclear scattering of protons of high kinetic energies are : 
(a) large angle coulomb scattering, 
(b) diffraction or shadow scattering, 


(c) inelastic nuclear scattering (Fernbach et al., 1949; Gatha and 
Riddell 1952). 


In a nuclear emulsion it is generally possible to separate inelastic 
events of the type (c) from events of the types (a) and (b). For events 
in category (c) a nuclear disintegration (‘ star’) is produced, and/or the 
incident proton loses energy in the scattering ; for coulomb and diffraction 
scattering there is no large change of energy, and therefore no appreciable 
change in ionization of the track produced by the proton following the 
scattering process. Thus in a photographic plate it is possible to identify 
events of the type (c), and to submit them to study without the undesirable 
inclusion of events of the types (a) and (b). However, the photographic 
plate has the disadvantage that it has more than one constituent element. 
This disadvantage can be partially removed by comparison of the effects 
produced in normal and diluted emulsions ; the stars observed can then 
be separated into two main groups, namely those produced in the light 
elements carbon, nitrogen and oxygen, and those produced in the heavy 
elements silver and bromine. 

In this laboratory an examination of the stars produced by protons 
with energies in the region of 130 Mev is being carried out. Separation 
of the stars into two nuclear groups, as a preliminary stage in their 


examination, will be described in this paper. 


* Communicated by Professor P. I. Dee, F.R.S. 


On the Cross Section for Star Production in Nuclear Emulsions 3065: 


§ 2. EXPERIMENTAL PROCEDURE 


(a) Exposure of the Plates 


Ilford normal and x 4 diluted G5 plates were exposed to the external 
proton beam of the Harwell cyclotron. The exposure was carried out 
using the early uncollimated and unducted beam. The protons entering 
the plate assembly possessed an energy of 146-+3 ev. After allowing 
for energy loss by ionization in the plates, the protons producing the 
stars may be considered to have kinetic energies between 110 and 146 Mev. 


(6) Hxamination of the Plates 


The plates were examined using x50 and x 90 oil immersion objectives 
and x6 eyepieces. ‘The search for stars was carried out both by the 
‘along the track’ and ‘ area’ scanning methods (Bernardini et al. 1952). 

A proton-induced star with n secondary charged prongs will be referred 
to as an np star. All scatterings of the type lp in which the ionization 
of the secondary particle was not statistically greater than the ionization 
of the primary proton were rejected, as diffraction and large angle 
coulomb scatterings appear as events of this type. With this criterion, 
251 and 122 stars were found by the ‘ along the track’ method in the 
normal and x4 diluted plates, for total proton path lengths of 8 460 
and. 7290 cm respectively. Thus the mean free paths for star formation 
were 35-1128 and 59-51§. cm respectively (the errors quoted throughout 
this paper will be the standard deviations). Before the cross sections 
for nuclear interaction in heavy and light nuclei can be calculated from 
the data given above, the mean free paths must be corrected for two 
factors. 

(a) Sometimes fast protons undergo collisions with free protons, present 
in the emulsions as the nuclei of hydrogen atoms. These collisions yield 
events of the type 2p. The contribution to be expected from events of 
this type was estimated using the published values for the cross section 
for p-p collisions (Christian 1952, Pickavance and Cassels 1952). This 
contribution was subtracted from the number of 2p events observed. 

(b) By rejecting all events in which a proton changes direction without 
observable change of ionization, some genuine inelastic collisions are 
rejected. Since the wavelength of the incident proton (~10~™ cm) is 
small compared with nuclear dimensions, the proton may be assumed to 
make its first collision with an individual nucleon inside the nucleus 
(Serber 1947). This collision can be of the form p—p (proton—proton) 
or p—n (proton—neutron), and following this encounter the energy of the 
two nucleons may vary from equality to the point where one has most 
of the available kinetic energy. If in a collision one of the nucleons 
receives most of the available kinetic energy, and this nucleon escapes 
from the nucleus without further collisions, the residual nucleus can 
receive only a small excitation energy and it is unlikely that slow charged 
particles will be emitted. In this case, if the fast secondary particle is 
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a proton, a Ip type event can then occur, but is rejected due to our 
criterion for the selection of stars. If the fast secondary particle is a 
neutron, the incident proton can appear to stop in the emulsion with 
no resultant prongs (‘Op star’). The recognition of Op events is 
unambiguous. A correction for the genuine Ip stars rejected was therefore 
made by estimating the ratio of the number of stars with a single secondary 
proton of high energy to the number of Op stars. This was done by using 
the published cross sections for free pp and n—p scattering (Christian 
1952, Pickavance and Cassels 1952), and assuming as a rough approxi- 
mation that the target nucleons could be considered as stationary 
in a potential well. A value of 2-5* was obtained for the ratio. Therefore 
2-5 times the number of Op stars observed was added to the number of 
lp stars accepted in the original measurements. 

Following these corrections, the mean free paths for star production 
became 32-6733 and 54*§ cm for normal and x4 diluted emulsions 
respectively. 


§ 3. SEPARATION OF STARS 
(a) Method of Prong Distribution in Normal and Diluted Emulsions 


In the conditions of our experiment the mean free path, d,,, for the 
production of a star with n prongs is related to the cross section for 
production by the equation 

1 

YY ET aye eee 6) oe ae ameme (1) 
where N,=number of silver and bromine atoms per cm? in the emulsion 
(Ny=2-06 x 10° and 0-92 102 atoms per cm? for normal and x4 
diluted emulsions respectively), N;—number of carbon, nitrogen and 
oxygen atoms per cm* (N;=2-69x 10” and 4:33 10% atoms per cm? 
for normal and x4 diluted emulsions respectively), oj —mean cross 
section for the production of a star with n prongs in silver and bromine 
nuclei, c; =the corresponding cross section for carbon, nitrogen and 
oxygen nuclei. 

Thus if 4, is known, the cross section oy, and o7 may be derived by 
applying eqn. (1) to normal and x4 emulsions. The total cross section 


for star production in light and heavy nuclei may be obtained in a 
similar fashion. 


eee 
* There is little data on the maximum energy that can normally be trans- 
ferred to a nucleon inside a nucleus without giving rise to the emission of a 
slow charged particle. An upper limit of 25 Mev was used in this calculation. 
This figure is suggested by the absorption of negative u-mesons by protons in 
nuclei. The maximum energy transferred to the nucleus in this case is about 
25 Mev (Tiomno and Wheeler 1949), and about 5% of the interactions are 
accompanied by the emission of slow charged particles (Fry 1952). The figure 
of 25 Mev also represented a rough lower limit on the energy loss that could 
be detected by grain counting on the incident and scattered tracks in lp events. 
The factor 2-5 was found to be not very sensitive to the value assumed for the 
maximum energy transferred in the collision. 
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In order to improve the statistics for the prong distribution, data on 
stars found by area scanning was added to that obtained by Blone the 
track scanning. It was found that although the mean free path obtained 
from area scanning was considerably longer (1-5) than that obtained 
from along the track scanning, the relative prong distributions were 
identical for stars other than of the Op type. We attribute this result 


Fig. 1 


NUMBER OF PRONGS 7 


Table 1 


Osi 


Type of star | (x 10? cm?) | ( 


to the difficulty of properly overlapping scans, when using high power 
objectives with very small fields of view. The data obtained from area 
scanning was therefore normalized to the overall mean free path obtained 
by along the track scanning, for all stars other than of the type 0p. The 
frequency distribution of prong size for stars produced by 130 Mev protons 
jn normal and x 4 emulsions is given in fig. |. The number of stars with 
n prongs is given in terms of the parameter 1/A,,. 
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Using the values for 1/A, given in fig. 1, oz, and o;, were calculated 
using eqn. (1). The results are shown in table 1; it can be seen that 
oz, only becomes larger than og, for 5 and 6 prong stars ats that for 
all practical purposes the 1, 2 and 3p stars are all produced in heavy 
nuclei. 

The overall cross sections obtained for star production were found to 
be (128-416) x 10-*6 cm? in heavy nuclei, and (16-47) x 10-26 cm? in light 
nuclei. 


(b) Alpha-particle Method 
Considerable work has been performed on the numbers of «-particles 
which emerge from stars formed by the cosmic radiation. It has been 
found that the ratio of «-particles to all tracks of energy less than ~30 Mev 
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is 0-27 for silver and bromine nuclei and 0-5 for gelatine nuclei (Rochester 
and Rosser 1951). We have examined the numbers of «-particle tracks 
appearing from stars initiated by 130 Mev protons in normal emulsions. 
In an analysis of this type an «-particle is given the overall definition 
of a particle whose track has a range greater than 5 microns, and whose 
charge is >2e; we have followed this convention. 

The ratio of the number of «-particle tracks to total number of prongs 
as a function of the latter figure is given in fig. 2. It can be seen that the 
ratio increases steadily with the number of prongs associated with a star 
and appears to flatten at the value of ~0-4 for 5 and 6 prong stars. The 
energy distributions of «-particles also show a difference for 1, 2 and 3p 
stars compared with 5 and 6p stars. This is illustrated in fig. 3. 
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The data presented in figs. 2 and 3 support the conclusion reached in 
the previous section regarding 5 and 6p stars—namely that they are 
mainly formed by the disintegration of the light gelatine nuclei. The 
energy distribution of «-particles from these stars shows a peak in the 
region of 4mMev. This result is consistent with the low potential barrier 
confronting the «-particles in light nuclei. In addition the value of 0-4, 
for the ratio of «-particles to total prong number for 5 and 6p stars, is 
in good agreement with that found for stars formed by the cosmic 
radiation in pure gelatine (Harding 1951). 


Fig. 3 
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In the 1, 2 and 3p stars the ratio of «-particles to total prong number 
is considerably lower than that found in stars formed in silver and 
bromine nuclei by the cosmic radiation. This result is to be expected 
since the cosmic ray stars involve a considerably greater excitation of 
the nucleus than that produced by 130 Mev protons. Most of the particles 
of energy <30 Mev ejected in the cosmic-ray stars are evaporation particles 
and have a nearly isotropic prong distribution, whereas in our case the 
ratio of the number of slow protons travelling forward in the direction 
of the incident particle to the number in the backward direction was 
found to be 2:5. Seed (1952) has observed a similar effect for stars 
produced by protons with a kinetic energy of 100 Mev; the result may 
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be interpreted in terms of protons directly knocked on in collisions. 
produced by the primary particle. The energy distribution of «-particles 
from 1, 2 and 3p stars has a peak in the region of 12 Mev. This result is. 
consistent with the high potential barrier presented to a-particles by 
silver and bromine nuclei. Thus the data obtained on the 1, 2 and 3p 
stars again supports the conclusion reached in the previous section. 

A rough cross section for star production in light nuclei was estimated. 
by assuming that all the 5 and 6p stars, and all the other stars with at. 
least one «-particle of kinetic energy less than 9 Mev were produced in 
carbon, nitrogen or oxygen (Menon, Muirhead and Rochat 1950). This 
yielded a value of ~15x10-°6 cm? for the total cross section for star 
production in light nuclei. This figure is probably an underestimate, 
since reactions of the type “C(p, pn)4C (Dickson and Randle 1950) 
would probably not be recognized as the disintegration of a light nucleus. 


§ 4. CONCLUSIONS 


Protons with a kinetic energy of 130 Mev produce larger stars in light 
than in heavier nuclei. Although the comparison of stars in normal and. 
<4 diluted emulsions allows a statistical separation of the stars to be 
made, it does not show the characteristics of the individual stars. The 
individual stars formed in the light nuclei may be recognized in many 
cases, however, by the high «-particle to total prong number ratio, and 
by the low energy of the «-particles emitted. 

It is of interest to compare the measured cross sections for star 
production in light and heavy nuclei with the ‘transparent nucleus ’ 
theory of Fernbach et al. (1949). These authors derive cross sections for 
the absorption and diffraction scattering of energetic nucleons by nuclei 
with the aid of an optical model, in which the effect of the passage of a 
nucleon wave through a potential well is examined. They derive a 
formula for the nuclear interaction of the following form : 


2K?R? 

where o=cross section for nuclear interaction, R—nuclear radius 
K=absorption coefficient (this term is the inverse of the mean foe 
path, Ao, of the incident particle in nuclear matter). A similar expression 
has been obtained by Bethe (1940). | 

; Perry (1952) has shown that K and hence o can be obtained from the 
single measurement of A (the total mean free path for star production in 
any emulsion), if a value of R and the correctness of eqn. (2) are assumed 
His method is equivalent to inserting the expression for o in eqn. (2) ate 
eqn. (1), giving 

while —2KR; l ! 
aus E {1 (1+2KR,)e ee see ey 5 (3) 

(the subscripts 7 refer to the atomic constitution of the emulsion). Since 
both NV; and A are known in this equation, K may be obtained if a abe 
figure for the nuclear radius is assumed. 


o=nk*| 1—{1—(14 2K R) e“ 2") : |: ee ees: 
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The term K is related to the cross sections for free proton—proton and 
proton—-neutron scattering by the expression 


ie 3A« [Zon t(A—Z) oy al (4) 
ML emrimincd we ge a AS koa 
where A=atomic weight of the nucleus, Z=atomic number of the nucleus, 
o,,—total cross section for free proton—neutron scattering, o,,—total 
cross section for free proton—proton scattering, «=exclusion principle 
factor which reduces the free particle scattering inside the nucleus. 

The quantities o,, and Fp», have been measured (Christian 1952, 
Pickavance and Cassels 1952), and « may be calculated with the aid of 
a Fermi gas model for the nucleus (Goldberger 1948). Hence K may be 
obtained if a value for R is assumed. In table 2 the value of K expected 


Table 2 


TOA)? oni! 


Nuclear radius | Energy of protons 


LOAS cm in Mev Transparent | Expected from 
nucleus eqn. (4) 

130 20000 2-6 

1:37 A18 220 ao 2:5 
375 Tyas 2-4 
130 2-0793 2-1 
220 1eGe 2-0 
375 1-1533 1-9 

Table 3 
Cross sections x 107° cm? oH OL 
Measured 128+16 16+7 
aeaneperct R1-37 A13x10-% em | 11043 ooo. 
Nucleus theory R 1-47 A1/3x10-%%em | 11178 2345 


from this formula is compared with that obtained from eqns. (2) and (3). 
Values for R equal to 1-37 A1/?x 10-43 cm and 1-47 Al? x 10-43 cm were 
assumed. | ; 

In addition to our data, values of K are also compared with those 
obtained from the published values for A at 220 Mev (Perry 1952), and 
at 375 Mev (Bernardini et al. 1952). In the latter case we have assumed 
that the values for the cross sections of the free proton—proton and 
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proton-neutron scattering at 375 Mev may be estimated. by extrapolation 
of the results obtained from measurements at somewhat lower energies. 

In table 3 the experimental values given in §3 for the cross sections 
for star production in light and heavy nuclei are compared with those 
obtained by substitution of K from table 2, column 3, into eqn. (2). 
The cross sections og and o, obtained by both methods are in satisfactory 
agreement, thus indicating that eqn. (2) is reasonably accurate for both 


heavy and light constituents of the emulsion. 
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ABSTRACT 


The character of the radio echoes obtained from the ionized trails 
produced by meteors has been investigated. The dependence on the 
polarization of the incident wave has been studied and the existence of 
resonant scattering with transverse polarization verified. The magnitude 
of the resonance and its dependence on the electron line density, «, is 
found to be in good agreement with theory, as is also the variation of echo 
amplitude with «. 


§1. INTRODUCTION 


Various authors (Herlofson 1951, Feinstein 1951, Kaiser and Closs 1952) 
have made theoretical studies of the scattering of electromagnetic waves 
by the ionized column of air created by the passage of a meteor through the 
upper atmosphere. All predict that, under certain circumstances, resonant 
scattering may occur when the incident electric vector is normal to the 
axis of the meteor trail (transverse scattering). According to the theory 
of Kaiser and Closs, the echoes fall into two categories, depending on 
whether the electron line density produced by the meteor is greater or less 
than a few times 10" cm~!. The present investigation is concerned with 
an experimental study of the echo characteristics for comparison with this 
theory, the main predictions of which may be summarized as follows :— 
(a) Electron line density, «, <10!*. If the incident electric vector is 
parallel to the trail axis (parallel scattering) it is sufficient to regard the 
electrons as an assembly of free scatterers, hence when the radius, 75, of the 
column of ionization is small compared with A/27 (A=free space wave- 
length) the scattering formula of Lovell and Clegg (1948) is valid. As ry 
increases (through. diffusion), the echo amplitude will decay in a manner 
depending on the variation of electron volume density, n cm~, with radius 
ry cm from the axis and time ¢ sec. If, initially, r7 <A/27 and if heating 
of the column may be neglected, solution of the diffusion equation gives 
D=—MrexPise Tilo) la eee a st we ex (1) 
where n)=«/m1'9, 79° =4Dt, D cm? sec-1=diffusion coefficient. The echo 
amplitude varies with time according to the law exp (—4k?Dt), k=27/. 
The duration to amplitude exp (—1) times the maximum is thus 
gee SL) ee ee eh 0 cece, cues (A) 


* Communicated by Professor A. C. B, Lovell, 
SER. 7, VOL. 44, NO. 350.—MARCH 1953 a 
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The assumptions made concerning the diffusion process lead to results 
consistent with observation. The effect of neglecting factors such as 
heating of the column and of electron loss through recombination and 
attachment will be considered elsewhere. 

With transverse scattering the reflection coefficient passes through a 
maximum at resonance, subsequently falling to the same value as for 
parallel scattering, followed by the exponential decay. The resonance 
occurs when the axial dielectric constant equals —1-4, and its theoretical 
magnitude, 7 (measured as the maximum ratio between echoes obtained 
from a single meteor with transverse and parallel scattering) is 2-0 when 
«<10!, decreasing to unity when « becomes a few times 10122 cm-!. We 
have reproduced, in fig. 1, the predicted ratio between the amplitude 
reflection coefficients, g,/g,, as a function of 4k?Dt/« for the case « <10”. 


Fig. 1 
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Theoretical polarization ratio, g,/gi, as a function of 4k?Dt/x, for a<10®. 


The additional assumption implicit in the above is that the meteor 
velocity is such that several Fresnel zones are traversed in a time small 
compared with 7; this is satisfied by the echoes obtained in the experiment. 

(6) «>10". The incident wave no longer penetrates throughout the 
column, which reflects similarly to an expanding metallic cylinder. When 
t>7, (79 >A/27), it is sufficient to assume total reflection at the surface of 
critical density, as in the case of ionospheric reflection, hence the echo 
persists until the axial electron density falls below the critical value. 
Except initially (t< 7), when g,<g, the theoretical reflection coefficients are 
equal and substantially constant in magnitude. In fact, echoes of duration 
greater than a few tenths of a second exhibit rapid amplitude fluctuations 
which Greenhow (1952) has ascribed to distortion of the trail b 
atmospheric turbulence and wind gradients. _ 

For an echo of class (a) the duration is equal to 7, given by (2), ie 
independent of line density, while the amplitude is proportional to y 


poe iieinaie 
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For class (b) the amplitude varies only as «1/4 while the duration, 7’ is 
greater than r, being proportional to x. Both rt and 7’ vary as A?D-1, 


§ 2. THE APPARATUS 


In order to test the above theoretical predictions it was decided to make 
independent measurements with transverse and parallel polarization of 
the echo amplitude obtained from a single meteor trail. The equipment 
was designed to operate during the activity of a major meteor stream with 
a fairly well defined radiant point, so that at any given time the meteors 
were travelling in parallel paths. 

The complete aerial system consists of an array of six half wave dipoles 
mounted A/8 in front of a conducting sheet which itself contains an array 
of six A/2 slots. A second conducting sheet is mounted X/8 behind the 
slot plane. The aerial is steerable in azimuth and elevation and during 
the experiment the fan-beam was continuously adjusted to maintain the 
polarization of the radiation from the slot and dipole arrays, parallel and 


Fig. 2 
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Short duration echo. 
0 Echo pulses received on slot array. 
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transverse respectively to the meteor paths. The minimum elevation of the 
aerial beam, 30°, is such that effects of ground reflection may be neglected. 
Preliminary experiments had indicated that an aerial of this nature would 
ensure identical polar diagrams for the two polarizations ; unfortunately the 
experimental results indicate that this ideal was not realized (see below). 
The transmitter output of 40 kw, peak power, in 8 psec pulses, is 
delivered equally to the two aerial arrays through twin T—R switches. By 
means of switched preamplifiers the slot and dipole array are connected to a 
600 ke/s bandwidth receiver in alternate intervals between transmitter 
pulses, and thence to an automatic recorder of the type described by 
Ellyett and Davies (1948). The twin cathode ray tube recorder display has 
a fast time base for range measurements and a slow time base (duration 
0-1 sec) which displays the individual echo pulses in sequence. A record of 
a short duration echo showing the two sets of echo pulses, is reproduced in 


fig. 2. 
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The mechanical design of the aerial system made it impossible to measure 
the polar diagrams by conventional methods. An overall check on the 
equipment was, however, attempted by attaching to a hydrogen filled 
balloon, a wire cage reflector so constructed as to be insensitive to polariz- 
ation. Two successful flights were obtained but owing to the short time 
for which the reflector was in range of the equipment, sufficient measure- 
ments could not be obtained to predict the variation of the error over the 
complete polar diagram. Measurements made over a range of +10° (ina 
vertical plane) from the supposed centre of the aerial beams showed an 
error which varied between 1-5 and 3-0 (in favour of the dipole array). 
It will be seen that this is consistent with the results discussed below. 


§3. EXPERIMENTAL RESULTS AND ANALYSIS 


The complete equipment was first operated during the Geminid meteor 
shower of December, 1950, with the results described by Clegg and Closs 
(1951). Almost all echoes appeared with greater amplitude on the dipole 
aerial which was interpreted as indicating enhanced tratisverse scattering 
due to the plasma resonance effect. The transverse/parallel echo ratio 
ranged between 2 and 10. 

Observations were continued in 1951 during the summer daytime 
meteor showers (Clegg et al. 1947, Aspinall e¢ al. 1948, Aspinall and Hawkins 
1951) and the analysis of these results is discussed in the following. The 
period covered is between June 4th and 12th when the predominant 
activity was from the Arietid radiant. In ‘order to eliminate partially 
the sporadic meteors and meteors which may have been observed in the 
aerial side lobes, echoes with ranges outside the limits 100-250 km were 
rejected as were also those for which the larger signal was less than five 
times receiver noise level, leaving a total of 140. These were divided into 
two groups (i) and (ii), the former being echoes which did not exhibit a 
decay during the length of the slow time base (0-1 sec). A consideration 
of the expected range of values of the diffusion coefficient in the meteor 
ionization region, showed that this should result in the echoes of group (i) 
being exclusively of category (b) of § 1, while group (ii) would contain a 
large proportion of category (a). 


Group (i) 


These echoes exhibited in all cases a substantially constant ratio, Ps 
throughout their lifetime, however p was often considerably greater than 
the theoretical value of unity, assuming the echoes to fall into category 
(6) of §1. A histogram showing the distribution of polarization ratio 
among the 46 echoes is given in fig. 3. It is our belief that the deviations 
from unity of p are due to experimental errors, mainly arising from defects 
in the slot aerial array, which result in a reduced sensitivity and a distorted 
polar diagram. This is substantiated in the following discussion, hence 
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fig. 3 may be interpreted as a distribution of experimental errors, in which 
case the 46 echoes correspond to electron line densities large compared with 
O22 cm 
Group (ii) 
The majority of these had an echo ratio which varied during the echo 
lifetime. Each echo was plotted and the ratio, p, was evaluated as a function 
of time, t. Figure 4 gives p vst for the echo of fig. 2. Ofthe 94 echoes, 45 


Fig. 3 
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Distribution of measured polarization ratio, p, for the 46 echoes of 
group (i) (long duration echoes). 


Fig. 4 
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Experimental polarization ratio, p, vs ¢ for the echo of fig. 2. Each point 
represents an average over several echo pulses. 


were of this general type for which p first increased to a maximum value, 
Pm, and then decreased to a fairly constant value, p... _Except for the 
deviations of p., from unity these are of the form predicted for echoes 
exhibiting the plasma resonance with transverse scattering. Of the 
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remainder, 31-were similar to the echo shown in fig. 5. In these cases, p Was 
either substantially constant or else increased monotonically to a steady 
value, p, While the echoes had an initial rapid rise, a short period fa 
fairly constant amplitude (in most cases) followed by an approximate y 
exponential decay. Except again for the deviations of px from unity, 
these echoes have the general form expected from line densities greater 
than a few times 10! em— namely category (6) of $1. 

Thirteen echoes of doubtful character were eliminated from this part 
of the analysis, and for 5 which did not exhibit a signal on the slot aerial, 
only a lower limit could be set for p... 


Fig. 5 
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(a) Non-resonant echo. 
(b) p vst for the echo of (a). 


Figure 6, which gives the distribution of p,, amongst these echoes is 
remarkably similar to fig. 3. To facilitate comparison, fig. 3 has been 
normalized by the factor 81 : 46 (the ratio of the number of echoes in the 
two groups) and reproduced as the broken line of fig. 6. The agreement 
strongly supports the assumption which we make in the following analysis 
that these distributions represent experimental errors. The ratios of fig. 6 
are slightly greater, on the whole, than those of fig. 3 which may, in part, 
be due to decay of the slot echo into receiver noise before the asymptotic 
limit, p., is reached. 

On these assumptions, the ratio p,/p,, represents the true resonant 
polarization ratio, 7. The values of 7 obtained for the 76 echoes are given 
in fig. 7. The 45 echoes which appear to exhibit the resonance effect show 
a maximum in the distribution corresponding to y= 1-7, in good agreement 
with the theoretical value of 2-0. 


—— 
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§4. Exectron Ling DENSITIES 


The existence or otherwise of the resonance enables the echoes of group 
(ii) to be divided into the two electron line density categories (a) and (b) 
of §1. In addition, since the resonance should occur when the axial 

Fig. 6 
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—_ Distribution of values of p,, for the echoes of group (ii). 
= {> — As for fig. 2, with ordinates increased in the ratio 81 : 46. 


Fig. 7 
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Distribution of values of the true resonant polarization ratio, y, for the echoes 
of group (ii). The shaded. portion represents the echoes exhibiting the 


resonance effect. 
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dielectric constant is —1-4, the line density of an individual trail may ie 
estimated, if the instant, t,, when p=p,, and the time constant, 7, of the 
exponential decay are known. The appropriate relation 1s 


a=2-14X 10" t,/7. ie Sr 


For the non-resonant echoes (category (0)) a similar estimate may be 
obtained from a knowledge of 7 and of 7’, the total echo duration. — Figure 
8 shows the dependence of « on 7’/7 derived from the theory of Kaiser and 
Closs. Using (3) and fig. 8, line densities were estimated in 66 cases. ; 

The echo amplitude is also a function of «, as well as of range (oak) 
and of the position of the meteor trail in the aerial polar diagram. The 
amplitude of the slot and dipole echoes were recorded in arbitrary units, 
reduced to a common range of 100 km and for the former group, corrected 
by the factor p,. In order to eliminate partially the polar diagram 


Fig. 8 
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avs7’/r according to Kaiser and Closs. 
7 =decay constant defined by (2). 
7’ =total duration to exp (—1) times maximum amplitude, 


effect, the results were divided into groups according to the estimated line 
density, namely, 1 —<2, 2 —<4, 4 —<8, 8 —<16, 16 —<32, 32 
—<64x10!1cm~! containing 4, 9, 17, 12, 9, 15 echoes respectively, 
and the line densities and corrected amplitudes were averaged for each 
group. The result is plotted in fig. 9; the amplitude scale of the 
theoretical curves given for comparison has been adjusted so that they fit 
the experimental points for the smaller line densities. The agreement is 
satisfactory with the possible exception of the group 16-32 x 1011 em-1.* 
The saturation of echo amplitude for line densities above 10! cm! is 
particularly apparent. 

The variation of with «, obtained from the averaged results, is given in 
fig. 10, together with the theoretical curve. 


* Due to various factors (e.g. finite meteor velocity and finite transmitter 
pulse interval) it is possible that some category (a) echoes arising from low 
electron densities, would not exhibit the polarization effect. The above 
analysis, however, would place these in the line density group 16-32 x 1044cm-1! 
with the effect of depressing the average amplitude for this group. 
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The amplitude scale of fig. 9 is such that the value at line density 
1-65 x 101 corresponded to a signal received on the dipole aerial with a 
maximum amplitude about 5 times the receiver noise level. The mean 


Fig. 9 
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Echo amplitude (reduced to range 100 km and corrected for experimental 
errors) vs electron line density. 
x Transverse scattering. o Parallel scattering. 
Theoretical curves adjusted to fit the experimental 
jae for small line densities. 


Fig. 10 
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range for this group of echoes was 130 km giving an estimate of the 
equipment sensitivity. An independent estimate may be obtained from 
the equipment parameters previously given. Assuming a receiver noise 
factor of 10db, a mean integrated sky background temperature of 
2000° kK, an aerial gain of 20 (referred to an isotropic radiator) and taking 
account of the polarization ratio, a 5x noise echo on the dipole aerial 
would require a line density of 1-5 x 101! at the above range. 


§ 5. Dirrusion COEFFICIENTS 


For 76 echoes of group (ii) it was possible to estimate the decay constant, 
7, and hence, using (2) the diffusion coefficient D. The results are given 
in fig. 11, the most probable value being 2-2-5104. This value should 
be related to the conditions of the atmosphere at the corresponding height, 
which, for the Arietid shower, has been determined independently by 
Evans and Kaiser (unpublished) at 88 km. 


Fig. 11 
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(e) 
Distribution of diffusion coefficient for the echoes of group (ii). 


A simple theory of ambipolar diffusion suggests that, as a result of the 
space charge electric field within the ionized column, the effective diffusion 
coefficient is related to the classical kinetic theory value of the ionic 
diffusion coefficient, D,, by 

D=D,;(1+T,,|T;) ; oS) everest 
where 7’, and 7; are the electron and ion temperatures respectively. 
Putting 7’ .=7;—T and assuming the ionic mass, m g, to equal the mass of 
the neutral air molecules (this should be approximately true, even for 

ionized meteor atoms) we obtain 


D 4 8kT ; 
= 3 d; are, =3-6x 108 A; Pr1/2 em? sec 1 4 4 a A (5) 


where A; cm=ionic mean free path, k—Boltzmann’s constant, m= 29 a.m.u. 
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Using the above data for D, with T~200° K, we obtain A,;~0-5 cm. If 
the cross section for ion-neutral molecule collisions is assumed to be the 
same as for collisions between neutral molecules (~10-15 cm?), the atmos- 
pheric particle density, N, can be calculated from the kinetic theory 
formula for \;. Table 1 compares the result with the values at 88 km 
height, interpolated from the data of other authors. It is seen to be 
closest to that of Whipple (1943); however an increase of about three 
times in the assumed collision cross section would bring it into agreement 
with the value obtained by Havens et al. from V-2 rocket flights. 


Table 1 


Source A B C D 
Particle Density, Ncm-? 41014 3x 1014 221014 1-2 1014 
A From decay of Arietid meteor trails. 
B- Whipple (1943). 
C Gerson (1951). 
D_ Havens et al. (1952). 


§6. SUMMARY 


Although the majority of Arietid meteors observed gave larger signals 
with transverse than with parallel scattering, this has been shown to have 
been partly due to errors inherent in the equipment. 

It proved possible, in the majority of cases, to allow for these errors and 
hence to derive the true polarization ratio. For a substantial fraction of 
the echoes, this ratio varied approximately as predicted for trails with line 
densities less than a few times 10” cm}, exhibiting the plasma resonance 
effect.. The most probable value of the resonant polarization ratio was 
1-7 in fair agreement with the theoretical value of 2-0. Almost all of the 
remaining echoes were of the form expected from higher line densities. 

Individual line densities were estimated independently of the equipment 
parameters. Both echo amplitude and resonant polarization ratio were 
found to vary with line density in the manner predicted by theory. 

The sensitivity of the apparatus, determined from the estimated line 
densities and echo amplitudes, is in good agreement with the value 
calculated from the equipment parameters. te 

From the decay of the echoes, the coefficient of diffusion of the Arietid 
meteor trails was determined. The most probable value lies between 


2-0 and 2-5 cm? sec !. 
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ABSTRACT 


Thin targets containing separated 1°B have been bombarded by deuterons 
of energies 0-59, 0-78 and 1-07 Mev, and the emitted alpha particles 
observed using a differential ionization chamber. Transitions proceeding 
via the ground state and first excited state of 8Be have been identified. 
On the basis of the dispersion theory of Wigner and Eisenbud, the shape 
of the alpha particle spectrum is used to indicate properties of the first 
excited state of ®Be. 


A stupy has been made of the groups of alpha particles emitted under 
deuteron bombardment of the separated isotope 1°B. Using deuteron 
energies between 0-59 and 1-07 Mev, alpha particles corresponding to 
transitions to the ground state and first excited state of §Be have been 
observed. The energy release in the ground-state transition has been 
measured. From the shape of the distribution curve corresponding 
to the excited state of *Be, properties of this level have been assigned 
on the basis of a form of dispersion formula due to Wigner and Eisenbud. 
Thin targets of B,O, (approximately 0-3mm air equivalent) were 
evaporated on brass backings by slowly heating boric acid in a vacuum. 
This boric acid was prepared originally from metallic 1°B. The beam 
from the 1-25 Mev Philips High Tension Set was allowed to strike the 
target through a 5/16 in. stop, and alpha particles were detected at 
90° using a differential ionization chamber (French and Treacy 1950), 
with an air absorption cell and mica foils as stopping media. The 
chamber subtended an angle of 16° at the target. Chamber and foils 
were calibrated with «-particles from ThC and C’, the mean ranges 
being taken as 4-72 and 8-57 cm respectively. The chamber resolution, 
after allowing for straggling, was found to be 0-25+-0-02 cm of range. 
Allowance was made for variation of stopping power of the mica foils 
with alpha energy (Livingston and Bethe 1937). The bombarding 
voltage scale was calibrated by using gamma rays from well-known 
resonances in the proton bombardment of 1°F (Chao eé al. 1950). 
Results are shown in fig. 1, for a deuteron energy of 0-78 Mev. Peaks 
labelled B and C are due to alphas accompanying transitions to the first 
ok I a ic pl Leet cet Bl laa ee ee ee 2 
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excited state and ground state of *Be respectively, and peak A to alphas 
from a contamination of carbon formed in the target in the evaporating 
process. Under bombardment the targets were maintained at a 
temperature of 120° 0, and there was no evidence for an increase of 


Fig. 1 
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Alpha particles emitted in the bombardment of 1°B by deuterons of energy 
0-78 Mev. 


Fig. 2 
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Alpha particles of ranges near 6 cm emitted in the B+-D process at deuteron 
energies of 0-59, 0-78 and 1-07 mev. 


carbon contamination throughout each run. Peak © was located 
accurately at different bombarding energies and, allowing for target 
thickness, occurs at a mean range of 15-05 cm for 0-78 Mev deuterons. 
With the range-energy relation of Livingston and Bethe (1937) this 
gives an energy of 12-33 Mev, and a Q-value of 17-91+-0-06 mev, for the 
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ground-state transition. The energy of *Be corresponding to the top 
of peak B is 2-88+0-08 mev, above the ground state. No well-defined 
peak, such as has been observed previously in this process, appears in 
the vicinity of 6cm range. This was confirmed in separate runs, shown 
in fig. 2, in which alphas were observed at deuteron energies of 0-59, 
0-78 and 1:07 Mev. (The peak appearing below 5 cm range in fig. 2 can 
be attributed to a trace of UB present in the target.) Separate 
bombardments at 0-78 mev, on cold 1°B and natural boron targets 
showed no evidence for further peaks. This agrees with results of 
Cockcroft and Lewis (1936) but not with those of Smith and Murrell 
(1939) and Whitehead (1951). Counting rates used in the present 
experiments were not sufficiently great for the detection, by ‘ pile-up ’, 
of protons at 53 cm range from the 160 (d, p) process (French and Treacy 
1950), which might account for this apparent discrepancy. 

If no transitions are involved other than those to the two lowest 
levels of *Be, it should be possible to describe completely the shape of 
the observed alpha particle spectrum between ranges of 24 and 14cm, 
assuming only a form of curve for the shape of the broad peak (B of fig. 1) 
and the angular distributions of the disintegration products in each stage 
of the process. Effects of straggling and chamber resolution may be 
neglected in so broad a peak. We have assumed that the peak follows 
the form 

Ta- Ib 
© (HgF AEP TT Ye? 
where 2ky* 
Lo= PG 
v2 ka OCHS) t2 
a meet ~~ d (ka) ) 
suggested by the work of Wigner and Hisenbud (1947) and Adair (1952). 
The width [a was assumed constant, as it varies slowly with the Q-value 
for formation of *Be, which is itself not a sensitive function of H, the 
energy of *Be with respect to emission of two alpha particles. In these 
formulae H, and y* are constants, k the wave-number of the reduced mass 
particle and / the spin of the excited state of SBe, assumed 0 or 2. Coulomb 
wave functions / and G, evaluated at an effective *Be radius a, for 
orbital angular momentum /, were obtained from tables (Bloch et al. 
1951). The peak was fitted at its low H side by a method of successive 
approximations, using values 4-48 10-"%cm for the radius a (Christy 
and Latter 1948) and 0-078 mev for the energy of the ground state of 
8Be (Carlson 1951). The following values of the parameters H, and y* 
were obtained : ~ 
For J=0: H,=2-90 mev, y?=2-29 x 10-8 Mev cm, 
_For 1=2: H,=5-29 Mev, y?=13-4 x 10-1 Mev cm. 


The shape of the continuum of alphas resulting from the disintegration 
of the excited ®Be was calculated in each case, assuming that all alphas 


Z2 


n (E) 


A= 
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involved are emitted isotropically with respect to their parent nuclei. 
The resulting total distributions are shown in fig. 3, together with the 
rimental energy spectrum. 
The Caner Sees agree better with the curve for j=2. ‘The 
discrepancy between either calculated curve and the experimental one 
may be due to a contribution from a direct disintegration into three 
alpha particles, or to the presence of a very broad peak in the continuum, 
not sufficiently intense to be observed in an experimental range spectrum. 


Fig. 3 
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Experimental energy spectrum (full curves) and theoretical distributions 
(dotted curves) for 8Be spins / of 0 and 2. 
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SUMMARY 


A new variational method is developed to treat the properties of slow 
electrons in ionic crystals, and similar problems. It is shown that the 
true energy of the ground state is always lower than the energy obtained 
from second order perturbation theory, however strong the interaction. 
The corresponding wave functions differ, however, from those of pertur- 
bation theory unless the interaction is sufficiently weak. For NaCl the 
energy of the ground state is found lower than by any previous method, 
and hence closer to the true energy. The effective mass is of the order of 
twice the electronic mass ; no trapping takes place, therefore. 


§1. INTRODUCTION AND DiscussION 


It is well known that slow electrons in strictly periodic crystal lattices 
behave in many respects like free electrons. A slow electron, however, 
polarizes an ionic lattice so that the field in which the electron moves is no 
longer periodic. Fréhlich, Pelzer and Zienau (1950, quoted as I) have 
pointed out that this energy corresponds to the self energy of the electron 
in a field theory if the polarization of the crystal is considered as a field. 
They found that although the energy of the ground state is depressed by 
this interaction, the electron together with the polarization attached to it 
can still move nearly freely through the lattice. This result differs from an 
earlier suggestion by Landau (1933) who neglected the dynamic properties 
of the lattice. He thus treated the electron as moving in a static field due 
to the polarization produced by the electron. This would lead to the 
trapping of an electron near a given point in space. Although such a 
trapped electron could still move through the lattice, its effective mass 
might be very large, in contrast to the result of I. 

The calculations in I hold only for relatively weak interaction. The 
purpose of the present paper is to develop a new variational method which 
applies also to stronger interaction. The result shows that for a very 
general type of interaction the energy of the ground state is the same as 
found from second order perturbation theory. Since a variational 
method has been used now, however, it follows that this energy is always 


* Based on E. R. A. Report L/T 284. 
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notes left by him.—H. Frohlich. 

{ Communicated by Professor H. Frohlich. 


330 M. Gurari on the Self Energy of 


higher than the true energy, however strong the interaction. The wave 
functions used in the variational method tend towards those of pertur- 
bation theory only in the case of weak interaction. They may differ 
from them very considerably for strong interaction. 

In the case of ionic crystals assumed to have a single vibrational fre- 
quency w/27 for long longitudinal waves, a static dielectric constant € 
and an optical refractive index «4, the energy H(0) of the ground state is 
then (cf. 2.29) found to be — W where (cf. I, 1.5) 


1 
W= 3(-—=) eu. s wall 9 207 
2i enw iae € 
Here w is an inverse length satisfying 
hole lim. Se Oe ed 


It should be noted (cf. I, 4.1) that w is connected with the frequency 
w ,/27 of residual waves by 


w=w,(ele,)?. «cts 8) ya lel le eee 


It will also be shown below that the ratio of effective mass m* to the 
electron mass ist 
m* W 
ie ee 
In the case of NaCl (for numerical values of parameters cf. Gatien 
leads to —-W=—0-17ev. This value is lower than both Markham and 
Seitz’s (1948) value of —0-13 ev obtained by Landau’s method, and 
Pekar’s (1949) value of —0-10 ev (obtained with an approximation which 
is good for large W/iw) and it thus represents a better approximation. Since 
hw=0-03 ev it follows from (1.4) that the effective mass is only about 
twice the electronic mass. It is thus confirmed that as was suggested in I 
no trapping takes place. 


(1.4) 


§2. DERIVATION OF THE RESULTS 
Consider the Hamiltonian used in I, eqn. (2.6), 


A= int Poll. (2.1) 
where Pp, is the momentum operator of the electron, 
LT eee HS. eons Mei (2.2) 


refers to the oscillators representing the lattice vibrations and H,,,is their 
interaction with the electron. The lattice vibrations can be expressed in 
terms of plane waves with wave number w, and H w refers to a single one 
of them. Clearly if |K) is the electronic wave function, 

Por |K)= ae 
2m 2m 


BG SII emer (3 


} From a preprint received recently I found that Tsung-Dao Lee and David 


Pines have independently derived the same results. They h 
published details—H. Fréhlich, See ee 
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where 7K is the electronic momentum. Also if | n,,) is an oscillator wave 
function 


Nie) Out, Tee UO | 1, yuk i (2.4) 
where the quantum number n,, is a positive integer or zero. Here we have 
omitted the zero point energy of the oscillator which can be done by an 
appropriate choice of the zero of energy. Clearly 

nh? K? 
(A—H,,,){| K) IT | n,,)}= Cae {]K) 7 |x} ie ue (2.5) 
w w w - 
which means that the product of the electronic and all oscillator wave 
functions are eigenfunctions of the operator H—H int Lhey form a 
complete set, obtained by using all combinations of all possible values for 
the quantum numbers n, and for K. 
It has been shown in I that the interaction term H,,, leads to transitions 
between these eigenfunctions provided the total wave number 


k=K+2 ny, w ger a are LP ew. G) 
is not altered. It is useful, therefore, to denote these eigenfunctions by 


|k;...m,...). Clearly in ordinary space (r) representation the 
normalized electronic wave function is 


hi -Ve exp (—1Kar)) ee se op. + (857) 
where V isthe volume. Thus using (2.6), 
[k;...my,-.- )=|K)J7|n,,)=V-# exp (ik.r) 7 exp (—in,w.r)|n,,) 


(2.8) 
are the eigenfunctions of H—H,,,.. Hence with (2.5), 
He 
Rem eae Ko Cte) | Ko 3 yes d= Sx + ++ Say my! + + Fy | 
X{((k—2 nyw)yP+2 nyu?} 2 . . (2.9) 


where use has been made of eqn. (1.2). Thus in the matrix representation 
based on the eigenfunctions (2.8) the operator H—H,,, has diagonal 
elements only. For the interaction operator H,,, on the other hand all 
diagonal elements vanish. The non-diagonal elements are different from 
zero only for transitions which leave k as well as all n,, but one unaltered, 
the excepted one differing by unity. Thus from I, eqns. (2.15) and (2.16) 


ee i eee tee Kg | Kon ily Nes oe = y Noe tL, 


CO [ty +1\t 
Ls Mey) Ct | Has | the 1)= 5 (A ) 5 (2.10) 


where using (1.2) 


C=(4nWhuf2m)t. 2. . . . . . (2.12) 


To find an approximate solution of our problem with the help of the 
variational method it is required that 


Sdb(k)|H[e(k))=0 2. 2 1 1. | (2,12) 
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where 5 represents a variation of the normalized trial wave function 
|¥(k)). The latter will be assumed of the form 


| u(k) = V-} exp (tk. r) JT|(k, w)) rae, hE 


where } ete ii; 144) 
0. + V-te(k, w, x) exp(—iw.r)| ly 
Kk, «))= aE) oa) 


Here the coefficients c(k, w, x) are assumed to depend on the magnitudes 
k and w and on the angle between k and w whose cosine is denoted by 2 ; 
they will be determined from (2.12). As required the wave functions (2.13) 
using (2.14) can be written as a linear superposition of eigenfunctions (2.8). 
They all have the same value of k and contain all combinations of the values 
Ny=0, and n,=1 of all w’s. We thus admit states which may contain a 
very great number of vibrational quanta provided they all belong to 
different wave vectors w. In contrast more restrictive wave functions were 
used in I where, as in perturbation theory, no states were used with more 
than a single quantum. 

Using orthogonality and normalization of the wave functions, eqns. 
(2.1), (2.13), (2.14), together with (2.3), (2.4), (1.2), (2.6), (2.7) and (2.10) 
yield 


Che [| Yl) = 5 fe +2 


(2.14) 


V-1c?(k, w, x)(u? + w?— 2kwe) 
1+ V-1c?(k, w, 2) 


es By Ve?(k, wy, %1)0°(k, We, %g)(Wy - We) \ 
Wi, Wa (ee V—1¢2(k, Wi; %))(1+- V¢*(k, We, Xp)) 
20 Vicky wa) 

= «wl + V-IC(k, w, @))° eR re eee 


Here the dash at the summation sign means that w,4w,. The last term 
arises from H;,,; the rest are due to the electronic terms apart from wu? 
which is due to H,,,. Clearly as V-> co, V-4c? can be neglected against 
unity (anticipating c to be independent of V). Furthermore, if polar 
coordinates in w-space are introduced with k as axis and ¢ as azimuth 
then only w, . w, depends on ¢, 


W1 . W2=W,W,(cos 8, cos 0.+sin 6 sin 0, cos (6, —¢,)). (2.16) 


Clearly 4 integration removes the cos (¢,:—¢2) term so that in (2.15), 
w,.W, can be replaced by 


WW, COS 0, COB Pg=W,X Wok, . . . . (2.17) 


For V — oo, eqn. (2.15) can thus be written as 


(hk) | | (ky) — EY) 


2m 


5 h? - 2Cc(k, 
+V 1 (w+ w?— 2khwa)c?(k, w, x) + ah . (2.18) 
where 


yaa Vie heaton le ee ee (2.19) 
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According to (2.12) the variation of (2.18) with respect to all coefficients 


c(k, w, x) must vanish. Hence 
(uw? + w*— 2kwa-+t 2ywx)e(k, w 2)=— 2m C 
wax Oe Oo 


which requires that 


(2.20) 


w+w—2(k—yjwrA0. . . 2... (2.21) 


Denoting by E(k) the expectation value of H when these values of 
c(k, w, x) have been inserted one finds from eqns. (2.18), (2.19) and (2.20) 
_Wike+y) 2m eek u? + w*?— 2khwa 
ey ee a rae yh 
2 eee) 
wer 2(k—y)\wal «am 
2m I 
ay ee, PY APS ee Sa ee Ee 
Om a OF Pig mee arap 
Introducing now for y the dimensionless variable 


Ga ie) MOC eae hoe a) os (Zone) 


(2.22) 


we find from (2.22) and (2.21) the following parametric representation 
of E(k), 


2m 1 


fh? hg 

2m\? Wx 
Aw hele 1772 a Ew eet aks DAS 
ee ( h ) red LE Sy , Ce 

To evaluate the sums they will be replaced by integrals, 
2ar 

a ——s 2 ee. 226 
V Basses cai dw de... (2.26) 


‘The w-integration can be extended to infinity as has been discussed in I. 
All integrations are elementary. With the use of (2.11) and (1.2) eqns. 
(2.24) and (2.25) then become 


B(k)=tw| (F-») |-w. eae (2.27) 


ue i\a 
and 
Ny coe es | aa satin 2.28 
yy —v= his OM (a sin ») . es ees ( 
Clearly if v=0 then from (2.28) also k=0. Hence 
E£(0)=—W. a eer 2.(229) 


Tf y+ 1 then from (2.28) k—> co and with (2.27) E(k) ~ 2hwk/u— o. 
'The case v>1 has to be excluded to satisfy condition (2.21). It cannot 
be expected that the choice of the wave function (2.13) and (2.14) permits 
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a successful treatment of energies E(k) which are above the ground state 
by more than iw. For states below this energy v<1. In this case (2.27) 
and (2.28) can be developed into a power series in 2, i.e. 


W 
E(k)=—W+tiw (250-0 (1+ aa) +...) . . (2.80) 


and 
k W 
5a (14 he) +... . s [et el recmmee (2.31) 
Eliminating v leads to 
hk? 
E(k)=—W+ om* (2.32) 


with a value of m* given by (1.4). Its validity requires 
E(k)—H(0)<hw. . . . . « « « (2.83) 


§3. GENERALIZATION FOR k=O 


Equation (2.29) shows that for k=0 the energy H(0) is the energy 
obtained from second order perturbation theory (although the wave 
functions in the two cases are equal only if W/hw<1). This result holds 
not only for the model used in §2 but also in a much more general case 
obtained by the replacements 


Clw > g(w), u—>u(w), eee tase (ETT. 


where g(w) and u(w) are arbitrary functions of w, provided no divergent 
integrals occur. This can best be seen by making the replacements 
(3.1) in eqns. (2.24) and (2.25). Then v=0 requires k=0 because the 
integrand in (2.25) becomes antisymmetric in w Equation (2.24) is 
then equal to the second order perturbation formula. This result can 
also be derived by a more elegant method which will be presented here 
because it will probably be of use in the discussion of excited states. 

Consider the matrix representation of the Hamiltonian H given in 
§2, and write H as the sum of two terms 


A=H HT, eo oe. ee oe 
such that H” is a diagonal matrix with elements 
u” i? 
oon Typ oe GEL ee tet ee = Fp (Zw Ww)? —2n,,w?} <aetoa) 
Then since the diagonal elements of H;,, vanish the diagonal elements 
of H are obtained from (2.9) by putting k=k’=0 and Nw=N,,. Use of 
(3.1), (3.2) and (3.3) then gives 


/ h? 
Coes twee [HI]... ye.) = 5 Sy (w®+-w?(w)), + eee Ord} 
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as diagonal elements of H’. The non-diagonal elements are the same as 
those of Hi, given by (2.10). Clearly then, the operator H’ can be 
written as 


H’'=ZH'(w), eee eee (315) 
where - 
q2 
(ih | H'(W) | tow = Mw 5— (wet Uw)), (3.6) 


are the diagonal elements of H’(w). All non-diagonal elements vanish 
except those next to the diagonal ; they are given by 


Cg | H'(W) | (ty + 1))=96w) (=) Spanner (a7) 


Thus H’(w) is identical with the Hamiltonian of an harmonic oscillator 
(apart from its zero point energy) with frequency h(w?-+u?(w))/2m whose 
origin, through a linear term, has been displaced by 


A= {16m>g"(w)/ Viit(w?--u2(w))?}*, ~~... (3.8) 


assuming its mass to be unity. If | %)'(w)) is its eigenfunction in the 
ground level, clearly 


H'(w) |! w))=— ppeatay IHW). 89) 
Using (3.5), therefore, 
PPG Viigo) == le |g) rea ses (3.1.0) 
where 
Pum eT AW) sn er eo > (3.11) 
and ‘ 
By — 7 ee el 


Thus the lowest eigenvalue of H’ is exactly equal to the expression for 
the lowest state of H calculated with the help of perturbation theory. 

To compare the eigenfunctions it will be noticed that the function 
| o'(w)) of the displaced oscillator can be developed in a series of eigen- 
functions | ,,) of the undisplaced oscillator. It is found that 


— 233 \"w 
| Yo! (w) =exp(—se8) EAE ngs = B18) 


where 
apenre I) 2 ae a _ (3.14) 
So BYR ew) 
It follows that terms with n,>2 can be neglected if Vo. The 


eigenfunction |,’ ) is then seen to be identical with | #(k) >, eqns. (2.13) and 
(2.14) if k=0 and if the coefficients ¢ are introduced from (2.20) considering 
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the replacements (3.1). Also, since for k=0 the c’s are independent of the 
direction of w, the diagonal elements (3.3) of H” vanish forn,—0 and n,=1. 


Clearly then 

(ho) A721 =O, ta atl cP en oaks 
and hence 

Ciba" | HL oy’ )=Ho oe ee ee 


It now follows from the variational principle (2.12) that E, is larger than 
the lowest eigenvalue of H. 

It is hoped that the present method will permit treatment of excited 
states by using solutions of H’ different from the ground state. Also 
use of H” as a perturbation should lead to an estimate of the accuracy of 
the above solutions. 
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THE influence of dislocations on the growth of crystals has been fully 
discussed by Burton, Cabrera and Frank (1951). It has been shown that 
in the case of normal growth (i.e. when the considered crystal face is 
completely surrounded by the phase out of which it grows) the exposed 
edge of the screw dislocation will wind itself up into a spirally terraced 
hill (vicinal faces), centred on the emerging dislocation. 

Salol (phenyl ester of salicylic acid) crystals, when grown from a 
solution, exhibit sometimes growth spirals of the Frank type (Amelinckx 
1952). In many cases, however, especially when grown slowly from 
heavy organic solvents, they show spirally terraced pits (hopper faces) 
not foreseen as a growth feature by Frank’s theory. On the basis of 
Frank’s theory spirally terraced pits could indeed only be generated by 
dissolution or evaporation, and not as a consequence of growth. 

It will be shown here how a dislocation mechanism can be devised 
to account for the occurrence of ‘ spiral growth pits ’. 

As is well known, hopper crystals are formed e.g. when crystals grow 
whilst floating on the surface of a solution or of a melt. In this case 
diffusion will make growth possible, but only on a relatively narrow strip 
along the crystal’s edge. Let us suppose that at a certain moment, 
when the supersaturation is low enough, a screw dislocation is created 
in such a floating crystal. Similarly, as for the Frank mechanism, growth 
will start at the exposed edge AB (fig. 1 (a)), but only a narrow strip 
parallel to the crystal’s edge will be developed, e.g. BC. When this 
strip approaches the crystal boundary at C another strip CD (fig. 1 (0) 
can develop which, when having reached D, will be able to go on growing 
towards E (fig. 1 (c)). At the end a growth layer will have surrounded 
completely the crystal face. After having reached B again, the crystal 
(at least its outer part) will have grown by a layer equal in thickness to 
the screw component of the Burgers vector of the dislocation (fig. 1 (d)). 
As in the meantime the crystal has grown laterally, this new layer will 
have its inner edge a little more outward than the previous one, and the 
result will be a spirally terraced pit (fig. 1 (e)). 

An example of a spiral pit on a salol crystal is shown in fig. 2 (Plate 12). 
Stepheights cannot be measured with extreme accuracy because multiple 
beam interferometry cannot be used in this case. Step-heights measured 
by means of two-beam fringes ranged from some 400 up to 18 000 A. 


* Communicated by Professor Dr. W. Dekeyser. 
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To decide whether a feature was a pit or a hill we made use of the 
asymmetry of the diffraction pattern at an edge when going off focus. 

Similar growth features have also been observed on citric acid 
crystallized from alcohol. 

The proposed mechanism and its experimental confirmation emphasize 
once more the important role played by dislocations in various types of 


growth. 


Fig. 1 


(e) 


Another instance where dislocations were though i 
role is growth from the melt. It has been aes ie pale Cause : 
and Frank (loc. cit.) that above a certain critical temperature fk the 
surface roughness s could probably account for the growth of cr eal 
so that neither two-dimensional nucleation nor dislocations ad ee 
necessary. Chalmers and Martius (1952) have proposed another mechani : 
for growth from the melt which is not based on dislocations <a 

On crystals prepared by slow cooling of molten salol we Bhsced that 
on almost every crystal growth spirals either of the Frank type or of th 
type discussed above were visible (figs. 3, 4, Plate 12). ie ate 


el 
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Two alternative conclusions are now possible. The first is that the 
surface roughness will play no role when growth-promoting dislocations 
are present. The second that 7’, is superior to the melting point in the 
actual case, a possibility which was already considered on theoretical 
grounds by Burton, Cabrera and Frank (loc. cit.). 

A striking feature for spirals from the melt is that the Burgers vectors 
are most of the time very large (up to a few thousand angstréms). 
Nevertheless no craters (Frank 1951) were observed, the dislocations being 
dissociated into ribbons of dislocations. For growth from solution there 
is evidence of craters in some rare cases, although the Burgers vectors 
were most of the time much smaller. 

It is hoped to publish elsewhere a detailed account of these 
observations and a more elaborate theory of the proposed mechanism. 
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Note added in proof—Thymol has now been found to be a second 
example-of a compound exhibiting the here described growth features 
(Votava, Amelinckx and Dekeyser, to be published in Naturw.). 
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ABSTRACT 


The formation of the unsaturated helium II film has been studied by 
measurements of superfluidity and adsorption isotherms on glass between 
1:35°K and 2:27°K. Whereas with decreasing percentage saturation the. 
onset of superflow is depressed to lower temperatures, the coverage shows. 
a discontinuity at the normal lambda-point of the bulk liquid. Formation 
of liquid clusters has been advanced as a possible explanation and the. 
existence of an anomalous surface tension effect has been suggested. 


A YEAR ago we reported on experiments in which the transfer rates of 
unsaturated helium IT films on a german silver surface were determined 
thermally (Bowers, Brewer and Mendelssohn 1951). These experi- 
ments have since been extended to film transfer over glass surfaces, and 
have been supplemented by measurements of the adsorption isotherms 
of helium from 0 to 95% saturation between 1:35°K and 2-27°R. 
Comparison of the flow phenomena on german silver and on glass, as 
well as the transfer of the saturated film over different surfaces (Mendel- 
ssohn and White 1950), indicates that the microscopic state of the 
surface very strongly influences the film, particularly at high percentage 
saturations. Since glass appears to be the only surface on which simple 
and reproducible phenomena of film flow have been observed, the adsorp- 
tion was studied on a glass surface. Excepting the earlier and rather 
qualitative adsorption observations of Kistemaker (1947), glass has not- 
been used in recent measurements of the adsorption isotherms. 

The results obtained on a glass surface of about 0-5 m2 area are shown. 
in fig. 1, where the adsorbed volume for different percentage saturations 
is plotted against the temperature. Our results differ from those of Long 
and Meyer (1949), and of Frederikse and Gorter (1950), obtained on 
Fe,O; and on steel, in that for constant percentage saturation the 
coverage rises with temperature until the lambda point (2-19° K) is reached, 
when there is a sharp drop in the volume adsorbed. 

The flow measurements of the unsaturated film on glass have confirmed 
our observation on german silver that for each percentage saturation, 
there exists a definite temperature, below 2-19°K, above which no film 


* Communicated by the Authors. 
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transfer takes place. This is in disagreement with the earlier report 
of Long and Meyer (1952 a) that even films of only 1} statistical layers 
show superflow up to the lambda point of the bulk liquid.* 

Combining the results of the flow experiments and of the adsorption 
experiments, a (dotted) line has been drawn in fig. 1 dividing the diagram 
Into a region where superflow takes place and one where no anomalous 
film transport occurs. It seems remarkable that (except for very low 
percentage saturations) the adsorption curves all break only at the 
lambda point of the bulk liquid, while the onset of superfluidity occurs 
at lower temperatures. The adsorption isotherm for 1-552° kK, shown in 


Fig. 1 


208} 


8} 


\bds. 
(c.at NTP) 


0:5 


17 Le 2-1 23 
Temperature °K. A pl. 


Volume of the gas adsorbed on 0°5 m? of glass for varying percentage saturation 
in dependence on temperature. The dotted line indicates the onset of 
superfluidity. 


fig. 2, has been measured with particular care in the region where the 
onset of superfluidity occurs, but no anomaly which is larger than the 
experimental scatter of the points was observed. 

A number of explanations suggest themselves for the curious fact 
that the films show a definite break in coverage at 2:19°K but become 
superfluid only at lower temperatures, depending on the percentage 


* However, these authors have quite recently, using our method, confirmed 
that the temperature of onset of superfluidity is depressed as the thickness 
of the film is decreased (Long and Meyer 1952 b). 
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saturation. First of all, one may think that the lambda phenomenon 
and superfluidity are not closely connected, but this would be out of 
keeping with almost all other observations on helium II. Secondly, it 
cannot be ruled out that the moving film differs in its properties from 
the stationary one, although, at least for the saturated film, the obser- 
vations of Henshaw and Jackson (1953) do not indicate this. Finally, 
one might regard the break at 2-:19°kK as due to the presence of bulk 
liquid. This cannot be due to ordinary capillary condensation, which — 
should be the same just below the lambda point as just above. However, 
it is possible that the formation of an unsaturated film covering a solid 
surface may be different in helium below the lambda point from that 
i with any other gas. There are indications that the de Broglie wavelength 


Fig. 2 


ONSET OF 


SUPERFLUIDITY 


" So aM 60 80 100 


% Saturation 
The adsorption isotherm for 1°552° x, 
Inset : The region between 68 and 76°/, saturation where superfluidity begins. 


d in helium IT is much greater than in helium I, As the pressure of helium 
over a solid surface is increased, it may therefore be energetically 
preferable to form clusters of liquid of the order of magnitude of d 
rather than to pile on further layers of adsorbed gas. On the basis of 
such a model it would be wrong, as has hitherto been done, to interpret 
throughout the measured coverages as a ‘number of layers’ of the 
unsaturated helium film. In fact the formation of the helium IT film 
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would proceed by the adsorption on the plane surface of a number of 
layers which is not greater just below the lambda point than just above, 
on which clusters of bulk liquid appear, which ultimately aggregate on 
saturation to the well-known film, with about 10~-® cm thickness. Using 
the relation proposed by one of us (Mendelssohn 1947), 


II 1 
T,=const. ( — z) : 


where / is the mean free path in helium I, we see that d must be at least 
of the order of magnitude of ~10-* cm in order to maintain the break 
of the curves in fig. 1 within 0-01° of the normal lambda temperature. 
The same formula can also be used to estimate the thickness of the 
adsorbed unsaturated film whose lambda point is depressed to tempera- 
tures indicated by the dotted line. For a percentage saturation of about 
55%, which has a lambda point of ~1-4°k, a thickness of ~10-7 cm 
is obtained. 

Another consequence of the suggested model of film formation should 
be an anomalous condensation in capillaries of the order of d below 
2-19°x, far in excess of that calculated from the measured (Allen and 
Misener 1938) surface tension. We might further expect that He, 
besides not showing superfluidity, should also not form thick films. 
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XXXIX. CORRESPONDENCE 


Electromagnetic Radiation from Radium E 


By D. G. E. Martin and G. Parry 
George Holt Physics Laboratory, Liverpool* 


[Received February 2, 1953] 


ABSTRACT 


A scintillation spectrometer has been used to study the electromagnetic 
radiation from a RaE source, to determine if a monochromatic line exists 
in the region of 85 kev. It appears that the intensity of such radiation 
is not greater than 1 quantum per 6000 beta disintegrations. The number 
of gamma quanta in the bremsstrahlung above 90 kev is estimated to be 
2-5 10-3 per disintegration. 


Accorping to Zavelski, Umarov and Matushevski (1949), the beta- 
spectrum of Rak is complex, 92% of the beta transitions going to an 
excited state of 21°Po with an energy of 85 kev above the ground state. 
They found a line in the beta-spectrum at 66 kev with an intensity of 
0-1% per disintegration, and they attributed this to the internal 
conversion of the de-excitation energy. The external gamma radiation 
was, however, weak. 

One would expect that this gamma radiation would be found in 
experiments on the inner bremsstrahlung of RaE. The spectrum of 
Madansky and Rasetti (1951) does show an anomaly at about 90 kev, 
but their curve for =P shows a similar anomaly in the same region. It 
therefore seemed worthwhile to investigate the RaE spectrum again. 

We have used an apparatus similar to that of Madansky and Rasetti. 
A source of about 200,¢e of Rak, prepared from material supplied by 
the Radiochemical Centre, was deposited on a thin Zapon film 
(100 ng/em?) and the gamma radiation from it measured by a cylindrical 
crystal of Nal(Tl) (83 em diameter x 1-8 em) used in conjunction with 
an E.M.I. type 5311 photomultiplier. The pulses from the photo- 
multiplier were fed to a scaling circuit through a linear amplifier and a 
single channel analyser, and energy calibration of the pulse heights was 
obtained from the gamma-ray lines of RaD and ®Co. The primary 
beta radiation of Rak was absorbed in perspex, and the gamma-rays 
were allowed through a hole of 1-25 cm diameter in a lead screen 4 cm 
thick. 


* Communicated by Professor H. W. B. Skinner, F.R.S. 
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Figure 1 shows the observed distribution from RaD and fig. 2 that 
due to the electromagnetic radiation from the RaE source corrected 
for absorption. From this data, there may be a slight sign of a mono- 
chromatic radiation in the region of 85 kev but its intensity would not 
be more than 1 quantum per 6000 disintegrations. This result appears 
to be in agreement with that of N ovey (1952) 


Fig. 1 


2Omm Omg 60) 6.90 
KEV 


Observed distribution due to radiation from RaD. 


Fig. 2 


Counting rate, arbitrary units 


KEV 


Observed distribution, corrected for absorption, due to radiation from Rak. 


There may be a non-uniformity in the curve at about 45 kev, but 
this could be entirely due to statistical error. To investigate it further 
would require an apparatus having better resolution, and a lower 
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absorption correction in this region. Certainly, less than 1 in 500 of the 
counts at 46 kev can be attributed to RaD contamination. 

It seems probable that practically the whole of the observed 
distribution is due to inner bremsstrahlung, and an estimate of the 
number of quanta above 90 kev is 2-5 10~* per beta disintegration. 
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The Virial Theorem in Free Electron and Generalized 
Thomas—Fermi Models 


By N. H. Marcu 
Department of Physics, The University, Sheffield* 


[Received February 3, 1953] 


Ty a recent paper (March 1952, referred to as I) the writer has discussed 
the virial theorem in the Thomas—Fermi (T.F.) and Thomas—Fermi—Dirac 
(T.F.D.) theories. These models refer to absolute zero, and have recently 
been generalized to apply when 740. The purpose of this note is to 
record briefly some results arising from an investigation of the form of 
the virial theorem in these generalized models. 

Marshak and Bethe (1940) have extended the T.F. theory to the case 
of low temperatures, whilst a full account of the temperature has been 
taken by Sakai (1942) and Feynman, Metropolis and Teller (1949, 
referred to as F.M.T.). 

In the case of a single centre F.M.T. have obtained the result 


27. OU = 3870, > 2) (1) 
where 7 is the kinetic energy, U the potential energy, v the volume and 


p the pressure. In the many-nuclei case, not hitherto discussed when 
T' £0, the methods of I give the result 


4a eV 
TJ = 9, 2 j ¢ 
27 +U = 313m (2mkT')>/ [tar (7 -1) r:. climes cL ~e,(he)soe) 
the notation being that used by F.M.T. and in I. 
One remark aside from the virial theorem is worth making here. 
It is that the entropy S in the T.F. model can be written down immediately 
from the free electron result (see for example, Wohlfarth 1950) as 


k oo 
S=— | arl (finf+(I—f)In (If) dk, . . (3) 
where ifs [ exp (ermine + n)+ Uli 
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Correspondence 347 


This procedure is much more direct than that of Brachman (1951) who 
obtained S by integrating the Gibbs—-Helmholtz equation. Brachman’s 
eqn. (7 b) follows easily from our expression.* 

The generalization of the T.F.D. model presents more difficulties. 
Yokota (1949) has given a discussion in which he takes the distribution 
function f(&, 7’) to be the Fermi—Dirac function. However, this is not 
correct (as Yokota in fact noted), and improved approximations to f have 
been obtained recently by Wohlfarth (1950) and Lidiard (1951). It is 
somewhat surprising to find that with both these approximations for f the 
virial theorem in the free electron model is exactly satisfied and may be 


written, BAe a pUer ee Ab =e, ic are! (4) 


where A is the exchange energy and p the pressure obtained from the 
free energy F' by means of the relation 


omer 


If, on the other hand, following Yokota, one takes f to be the Fermi—Dirac 
function, then the theorem is no longer satisfied. 

We can improve the T.F.D. model of Yokota if we use Lidiard’s 
approximation for f. The equations obtained by minimizing the free 
energy are somewhat complicated and we do not propose to give them 
here.f But again in the case of a single centre one obtains after some 


manipulation 27 Ue Ase spv oe. |, (6) 


That the virial theorem should be satisfied when approximate 
distribution functions are used seems to be closely related to the fact 
that the free energy has been minimized with respect to the parameters 
in the forms chosen for f, but up to the moment we have achieved no 


general result. 
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Delayed Disintegration of a Heavy Nuclear Fragment : A ee 


By M. Danysz and J. PNInWSKI 
Institute of Experimental Physics, University of Warsaw 


[Received December 1, 1952] 


A REMARKABLE coincidence of two events recorded in a photographic 
emulsion has recently been observed in this laboratory. It occurred in a 
G5 emulsion, 600 » thick, which had been exposed to cosmic radiation 
at an altitude of 85 000 feet,t and consists of two stars marked A and B 
in the photo-micrograph reproduced in Plate 13. The centre of the star B 
coincides with the end of the track of a heavy fragment ejected from the 
star A. If this coincidence is not accidental, it must be considered as an 
example: of the delayed disintegration of a heavy fragment. The pro- 
bability of a fortuitous coincidence is very small, and it therefore seemed 
appropriate to analyse. the events more closely. It is clear, of course, 
that any novel conclusions drawn from a single observation should be 
treated with proper reserve. 

The disintegration referred to as star A represents the result of a high- 
energy interaction, presumably between a singly-charged particle and a 
silver or bromine nucleus. In the nomenclature of Brown et al. (1949), 
the star A may be described as of type 21+18p. The track p—see 
Plate 13 —is the only track in the upper hemisphere with minimum grain 
density, of which the direction coincides approximately with the ‘ axis’ 
of the cone of the emitted shower particles. It therefore probably 
represents the incident particle which produced the disintegration A. 
Scattering measurements do not allow an accurate evaluation of its 
energy, but they are not inconsistent with the above assumption. Using 
statistical data of Camerini et al. (1951), and of Daniel et al. (1952), we 
estimate the energy of the primary particle to be about 30 Bev. Energy 
measurements on some of the longer tracks produced by the shower 
particles are also in agreement with this estimate. 

From the analysis of the ‘ black’ and ‘ grey ’ tracks of star A, we could 
identify nine «-particles, one heavy fragment and eleven singly-charged 
particles. As no slow mesons were detected, the latter were probably due 
to protons, deuterons or tritons. The comparison of the shape of the 
track ‘f’ with those of other multiply-charged particles stopping in the 
emulsion allows a rough estimate of the charge of fragment f; the most 
probable value is about 5. 

The initial kinetic energy of the fragment ‘f’, estimated from its 
range (90 #) and charge (Perkins 1950), is of the order of 60 Mev. The 


fragment appears to stop in the emulsion exactly in the centre of star B. 
ee a et eee 

*A preliminary report of this observation was given by Professor A. Soltan 
at the meeting of the Polish Academy of Science held on 20th October 1952. 

i Communicated by Professor C. ¥. Powell, F.B.S. 

t We are indebted to Professor C, F. Powell for the opportunity of exposing 
the plates in a balloon flight and processing them at Bristol. 
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Tts final energy, as deduced from the characteristics of the end of its track 
1s very small and certainly much legs than that necessary to produce the 
star B. The time of flight of the fragment was of the order of 3x 10-22 sec. 

The results of the analysis of the star B, which is a four-prong star of 
low energy are given in the following table : 


Track Range | Grain density Identity Energy 


i: 9 pw Black p, d, t, or « 0-7 Mev 
assumed proton 


123 uw Black p, d, t, or « 16 Mev 
assumed «a-particle 


— 3°7+0-2 py ds t0r 7 82 Mev 
assumed proton 


2 ps . Black Recoil — 


Taking into account the binding energy of the emitted particles, we 
find for the total energy released in event B, a value equal to or greater 
than 120 Mev. We may, however, expect that some neutrons should also 
be emitted. If so, the total internal energy of the fragment should be 
increased to ~140 Mev. If track 3 was produced by a 7-meson, the total 
energy, including the rest-mass of the particle, is of the order of 180 Mev, 
not taking into account possible neutron emission. 

The probability of a chance coincidence between the end of the track f 
and the centre of star B, in the conditions of the experiment, was calcu- 
lated from, (a) the average star density in the exposed plate and the local 
increase in the radiation density, due to the proximity of star A; and, 
(6) the known optical resolution of the microscope. This probability is 
less than 1 in 107. 

CONCLUSION 


Assuming that the event is not due to a chance coincidence, we are left 
with various alternative possibilities. It might be attributed either to an 
interaction between a heavy fragment and a nucleus of the emulsion, or to 
the spontaneous decay of the heavy fragment (Schopper 1947, Lovera 
1947, Hodgson and Perkins 1949). The first interpretation fails because 
of the small, if not zero, final kinetic energy of the fragment. For the 
second interpretation to be valid, the fragment must have been emitted 
with a high internal energy, at least 120 Mev and probably more. Further, 
it must have remained stable, against both y-transitions and the emission 
of particles, during a time greater than 3x10~™sec. These consider- 
ations make it difficult to interpret the event in terms of a highly excited 
state of the nucleus. . 

It might be supposed alternatively, that the explosion was due to a 
m-meson capture at B, the meson being picked up in a Coulomb orbit 
round the heavy fragment as the latter left the disintegration at A. 
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B would then be regarded as a kind of ‘ delayed ’ o star. The weight to 
be given to this assumption depends on estimates of the probability of 
the heavy fragment picking up the meson in the disintegration A—if 
such a process is indeed possible— and of the time likely to elapse between 
the instant of capture of the meson into the orbit and its interaction with 
the nucleus. This time interval is generally considered to be of the 
order of 10~1? sec or less. 

An alternative explanation of the event may be sought in terms of the 
heavy neutral V,°-particle, or of similar charged particles, which may be 
considered as a nucleons in excited states, with mean lifetime greater than 
10-19 sec. It is possible that such particles exist not only as free particles, 
but also in bound states within nuclei. If the fragment were formed with 
such an excited particle among its nucleons, this could perhaps account 
for the delayed disintegration as well as for the observed release of energy. 
The kinetic energy Q, released in the decay V,;°>p*+7a would be 
augmented by the rest-energy of the created 7-meson, if the latter were 
absorbed in the same nucleus. 
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Delayed Disintegration of a Heavy Nuclear Fragment: II 


By D. A. Trpman, G. Davis, A. J. Herz and R. M. TENNENT 
Imperial College, London, 8.W.7 * 


[Received December 15, 1952] 


Soon after Dr. Menon of Bristol University had informed us about the 
event found by Danysz and Pniewski (see the preceding note), a similar 
observation was made in this laboratory. Photo-micrographs of the 
event are shown in Plate 14. 

The event occurred in a G5 emulsion, 400 » thick, exposed in a 
balloon flight at geomagnetic latitude 40° N. The balloon floated for 
eight hours at about 96000 feet. The parent disintegration (star A) 
is of type 16+ 0p, excluding the heavy fragment f. The track of the 


* Communicated by Dr. E. P. George. 
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ejected fragment has a length of 219 p, and its end coincides with the 
centre of a three-track star. The following table shows the results of 
measurements on the tracks. The references given refer to the method 
of analysis. 


(i) Fragment f. 


Range si oh a 28 a a> 219 fb 
Charge determined by comparison with tracks of 
identified particles—see Plate 15 % -. 2or3e 
Number of delta-rays of more than 15 kev .. None observed 
(ii) Star A. 

Energy (excluding fragment f) (Brown et al. 1949) -. 1620 Mev 
Kinetic energy of fragment f, if a helium nucleus a4 20 Mev 
Total visible energy released in star A .. tes .. 1640 Mev 
Energy of primary (Camerini et al. 1951) ue .. 6000 Mev 


(iii) Star B. 
Energies of particles 
Kinetic energy in Mev 


Particle Range Grain density fe retonmnnt opartiels 
Bl 24 w black 1-6 6-4 
B2 NOW black 4-0 16-0 
B3 leaves emulsion 109 gr/100 uw 4843 
(5:1 Simin) 
Total visible energy 5443 70+3 


Scattermg measurements on track B3 gave &199=(0-4-+0-15)°, and 
this suggests that the track is that of a proton rather than that of a 
a-meson. If the particle were a 7-meson its energy would be 6 Mey. 

It is of interest to make a rough estimate of the probability that the 
star B occurred by chance at the end of the range of the ejected particle 
from star A, the two disintegrations being unrelated. It seemed best to 
estimate a safe upper limit for such a chance juxtaposition. This involves 
estimating the probability that the centre of a star of 3, 4 or 5 tracks 
shall occur inside the cubic micron of emulsion at the end of the range of 
a heavy particle ejected from a star of six or more tracks. We assume 
that every star of six or more tracks contains one heavy fragment ending 
in the emulsion (Perkins 1950). 

Lord (1951) gives the rate of production of stars of 3, 4 and 5 tracks 
at great altitudes as 1000 cm~%d-!. The average duration of balloon 
flights is not more than about eight hours. It follows that the probability 
of observing a small star in a specified cubic micron of emulsion is 
SB 3< 10-10. 

The rate of production of stars with six or more tracks is given by 
Lord (1951) as about 1600 cm-8d-!. Assuming that 1000 cm? of 
emulsion have so far been examined carefully in laboratories all over the 
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world, we find that 5-3x10*® large stars, each containing a heavy 
fragment, should have been observed. The probability of observing 
one chance coincidence would thus be 10-4. The probability of finding 
two such chance events in the 1000 cm? of emulsion would be 10°°. 

It appears certain that the charge of the ejected heavy particle from 
star A is not greater than 3e (see Plate 15). It is rather unlikely that. 
the event can be attributed to the capture of a negative 7-meson by the 
heavy fragment (Danysz and Pniewski, preceding letter), for the visible 
energy in star B is somewhat high for such an event (Menon, Muirhead 
and Rochat 1950). The observations are, however, not inconsistent with 
the view that one of the nucleons of the fragment was a V,°-particle. 
The observed release of energy among the charged particles can- be 
reconciled with this assumption and the measured mass of the V;°-particle. 


f 
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A Neutral V-Particle from a Low-Energy Nuclear Interaction 


By J. P. Astsury, J. 8. Bucnanan, G. D. Jamus, D. D. Minuar,, 
J. A. Newta, D. I. Pacs, A. B. Sanrar and A. Rytz 
The Physical Laboratories, The University, Manchester 13 


[Received February 3, 1953] 


THE study of low-energy interactions in which V-particles are produced 
has a double importance. In the first place, it can be expected to provide 
fundamental information about the production process—whether, for 
example, the V,°-particle should be regarded as an excited nucleon or 
as a produced particle. Secondly, in the design of a large particle 
accelerator, the likely yield of V-particles and heavy mesons must be a. 
major consideration. This yield can only be estimated at present from 
observations made on cosmic-ray nuclear interactions. 

In a recent paper (Astbury et al. 1952) we have drawn attention to: 
the fact that charged V-particles can be produced in nuclear interactions. 
where the total energy involved is ~5x10%ev. In this letter we 
describe an event which shows that at least one type of neutral V-particle 
can be produced in interactions of similar energy. 

During the course of the cloud chamber experiment previously reported 
some 25000 photographs have been taken. Among these there are 
fifty showing characteristic nuclear ‘ stars ’ in the gas of the cloud chamber. 
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Thirty of these stars correspond to low-energy disintegrations in which 
no fast charged particles are produced. In the remaining twenty events 
there is evidence for the production of z-mesons or of fast protons. 

The photograph of one of these last events (see Plate 16) shows the decay 
of a neutral V-particle at a distance of 3-5em from the nuclear inter- 
action. Figure 1 shows a scale drawing of the event in which all the 
relevant tracks are numbered; the measurements made on them are 
given in the table. From the geometry and dynamics of the decay it is 
clear that the V-particle came from the star in the gas. 


Fig. 1 


Diagram of the interaction showing the primary particle (1) and the eight 
charged products (2-9). X and Y are the decay products of the V-particle 
produced in the interaction. 


If the star is related to the main shower occurring in the lead above 
the cloud chamber either particle 1 or a nearby neutral particle must be 
the primary causing the interaction. We consider it most improbable 
that particle 1 should be a secondary particle from the star since only one 
fast charged particle is emitted in the forward direction and all three 
identified light mesons move backwards in the laboratory system. 
Were particle 1 a secondary, this asymmetry would be even more pro- 
nounced. For this reason, our interpretation of the event assumes 
particle 1 to be the primary. 
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Excluding particle 7, the energy given to the charged secondary 
particles from the star is ~ 10° ev; this figure includes the rest-energy 
of three 7-mesons. From measurements made on stars in photographic 
emulsion (Camerini, Lock and Perkins 1952) it is known that the average 
energy of shower particles in low multiplicity stars is less than 10° ev. 
It is reasonable to take this as an estimate of the energy of particle 7 
although the photographic plate work shows that there are large fluctua- 
tions about the average value. 


Table of Measurements 


Track Sign and Tonization Nature of Kinetic 


* 
number m (108 a density (I/I) particle energy (Mev) 


¢ >5 a 

0:6 7-meson 

+2:0 y 7-meson 

aoe proton 

+3°6 proton 

Ae proton 

pa 2 

+2-0 proton 

+3:2 7-meson 

—1-0 : light meson See text 
>3 2 See text 


1 
2 
3 
4 
5 
6 
7 
8 
9 
x 
a7 


* No errors are given for the measurements of momenta since their effect is. 
negligible in computing the energy of the interaction. 


The total energy of the neutral mesons and neutrons emitted from the 
star is unlikely to be very much larger than that of the protons and 
charged mesons. What evidence there is suggests that the number of 
7-mesons from showers of low multiplicity is about the same as the 
number of charged mesons (Camerini et al. 1951). The numbers of 
neutrons and protons must, of course, be very nearly the same. 

Owing to the shortness of the tracks X and Y, the V-decay cannot 
be fully analysed. The plane of the fork can only be determined to. 
within 5° and, to this accuracy, the star is coplanar with the decay 
tracks. The angle between the tracks is 55°+5° and the measurements. 
are all equally consistent with the decay of a V,°- or a V,°-particle 
(Armenteros et al. 1950). Assuming a V,°-decay, the eineee energy 
oe pace and time of flight of the particle were about 390 mev, 1 500 Mev 
ae cae ; rae eit figures for a V,°-particle are 330 Mev, 

From the preceding argument we consider that the best estimate of 
the energy of the interaction is ~6 x 10° ev and that the value is uniikel a 
to be much greater than 10? ey. f : 
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It is interesting to consider how this event would have appeared had 
it occurred in a lead plate in a cloud chamber or in photographic emulsion. 
Only the secondary particles 3, 7 and 9 have ranges in lead greater than 
one centimetre and, of these, 3 has a nearly horizontal trajectory and 
would probably not emerge from a horizontal lead plate 1 cm thick. 
Thus the only secondary particles seen would be 7 and 9 and, of course, 
the V-particle.. 

If it had occurred in photographic emulsion the event would have been 
classified as a star with three shower particles (cither 3, 7, and 9 or 1, 3, 
and 9 depending upon which of the particles 1 and 7 were taken as 
primary) and the energy roughly estimated at 6 to 10x16%ev. The 
V-decay would not have been observed. 


The photograph discussed here was obtained at the Hochalpine 
Forschungsstation, Jungfraujoch, Switzerland. We are indebted to the 
administration of the station and to Mr. H. Wiederkehr for enabling us 
to work there. 

We are grateful to Mr. R. W. Wilkins and Mr. A. H. Chapman for 
assistance in running the apparatus and to Miss P. Miles for making many 
of the measurements. 
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XL. Notices of New Books and Periodicals received 


Light. By R. W. Drrconpurn. (Blackie & Son Ltd.) [Pp. 680.] Price.45s. 


Proressor DrrcHBURN’S book is very well produced and is a pleasure to read ; 
it is a valuable addition to the range of post-war English textbooks. The 
aim stated is ‘‘ to describe the wave theory in such a way that the quantum 
theory may appear as a natural development rather than as an alternative 
theory”, and the treatment throughout adheres closely to this ideal. The 
book covers the requirements of a full Honours Degree Course and has been 
arranged to be read through twice. The material for omission on a first 
reading is clearly marked, but it may perhaps be doubted whether a student 
starting with only an ‘ Intermediate ’ knowledge of Mathematics could under- 
stand all the earlier analysis. 

The material discussed has been very thoughtfully selected, the sequence 
in which topics are arranged presents a satisfactorily coherent picture, and, in 
general, the book can be strongly recommended as a thoroughly suitable basis 
for an Honours Course. S.H. P. 


Theoretical Nuclear Physics. By J. M. Buatr and V. F. WEIssKopr. 
(Wiley & Sons.) [Pp. 864.] Price $12.50. 


THE enormous development of nuclear physics during the last fifteen years 
has made obsolete most of the books and reviews written during the 1930's. 
Several new phenomena have come to light, and many others are found to 
require fresh theoretical treatment. This book deals mainly with the physics 
of nuclear processes involving energies less than about 50 Mey, a notable 
exception being the analysis of n—p and p-p scattering up to energies of 350 Mev. 
It is essentially, and in the best sense, an elementary book, dealing more with 
the broad physical principles of nuclear physics than with detailed mathematical 
solutions ; nevertheless, several sections are devoted to more mathematical 
aspects of the theory. There is full discussion of the experimental data, and 
some of the interpretation has not been published elsewhere. 

It is easy to find fault with a work of such magnitude written while the 
subject matter is still growing and changing. For example, no account is 
given of the fruitful Butler theory of (d, p)-reactions, while all mention of the 
important interference terms in f-decay is relegated to a footnote. There is 
a general tendency towards repetition, especially through the use of over- 
simplified models to pave the way for a more exact and complicated treatment 
later on. Thus, the theory of f-decay is first formulated non-relativistically : 
this might be helpful if it enabled me to avoid a matrix formulation of the 
problem, but naturally it does not. Again, a one-particle model is used for a 
first discussion of the dispersion theory of nuclear reactions in spite of the fact 
that much the same results can be obtained with greater validity and hardly 
more complication from a Kapur-—Peierls treatment. 

It is a pity that very little recent work published outside the Physical Review 
has been included: the authors’ apology in the preface does not excuse this 
omission in a work which promises to be a standard text. And it is strange 
that in an expensive book of over 800 pages, ten per cent of the space is devoted 
to lists of symbols. Nevertheless, the book gives an excellent comprehensive 


picture of low energy nuclear physics, and no nuclear physicist today can 
afford to do without it. SoHo 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


J. P. ASTBURY et al. 


Phil. Mag. Ser. 7, Vol. 44, Pl. 10. 


a 


LALO AR RUNES. agate xing j 
ee 


(a) Event DP22. A negative heavy particle comes from a low-energy 
interaction in the lower lead plate. (See table 1.) | 


(b) Event JM41. A negative heavy particle from a nuclear interaction above i 
the cloud chamber stops in the lead plate. (See table 1.) 
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(a) Event KD/9. A negative heavy particle with a time of flight in the 
cloud chamber of 3x 10-® sec. (See table 1.) 
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(b) Event CC17. A slow particle whose direction of travel cannot he 
determined definitely. If it is moving downwards its charge is negative. 
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(a) (0) (c) 
Track (a) Lithium nucleus ( 6) 
(6) Heavy fragment f. 
(¢) «-particle. 
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